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1 Introduction 



These notes summarize the lectures presented at the 2005 CERN-CLAF school in Malargiie, Ar- 
gentina and the 2009 CERN-CLAF school in Medellm, Colombia. The audience in both occasions 
was composed to a large extent by students in experimental High Energy Physics with an important 
minority of theorists. In nearly ten hours it is quite difficult to give a reasonable introduction to a 
subject as vast as Quantum Field Theory. For this reason the lectures were intended to provide a re- 
view of those parts of the subject to be used later by other lecturers. Although a cursory acquaitance 
with th subject of Quantum Field Theory is helpful, the only requirement to follow the lectures it is 
a working knowledge of Quantum Mechanics and Special Relativity. 

The guiding principle in choosing the topics presented (apart to serve as introductions to later 
courses) was to present some basic aspects of the theory that present conceptual subtleties. Those 
topics one often is uncomfortable with after a first introduction to the subject. Among them we have 
selected: 

- The need to introduce quantum fields, with the great complexity this implies. 

- Quantization of gauge theories and the role of topology in quantum phenomena. We have in- 
cluded a brief study of the Aharonov-Bohm effect and Dirac's explanation of the quantization 
of the electric charge in terms of magnetic monopoles. 

- Quantum aspects of global and gauge symmetries and their breaking. 

- Anomalies. 

- The physical idea behind the process of renormalization of quantum field theories. 

- Some more specialized topics, like the creation of particle by classical fields and the very 
basics of supersymmetry. 

These notes have been written following closely the original presentation, with numerous 
clarifications. Sometimes the treatment given to some subjects has been extended, in particular the 
discussion of the Casimir effect and particle creation by classical backgrounds. Since no group 
theory was assumed, we have included an Appendix with a review of the basics concepts. 

By lack of space and purpose, few proofs have been included. Instead, very often we illustrate 
a concept or property by describing a physical situation where it arises. Full details and proofs 
can be found in the many textbooks in the subject, and in particular in the ones provided in the 
bibliography [1-10]. Specially modern presentations, very much in the spirit of these lectures, can 
be found in references [4, 5, 9, 10]. We should nevertheless warn the reader that we have been a bit 
cavalier about references. Our aim has been to provide mostly a (not exhaustive) list of reference 
for further reading. We apologize to those authors who feel misrepresented. 

Acknowlegments. It is a great pleasure to thank the organizers of the CERN-CLAF 2005 and 
2009 Schools, and in particular Teresa Dova, Marta Losada and Enrico Nardi, for the opportunity to 
present this material, and for the wonderful atmosphere they created during the schools. The work 
of M.A.V.-M. has been partially supported by Spanish Science Ministry Grants FPA2002-02037, 
FPA2005-04823 and BFM2003-02121. 
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1.1 A note about notation 

Before starting it is convenient to review the notation used. Through these notes we will be using 
the metric rj^ = diag (1,-1,— 1,-1). Derivatives with respect to the four- vector x M = (ct, x) will 
be denoted by the shorthand 

d (Id ~\ 

As usual space-time indices will be labelled by Greek letters (fx, u, . . . = 0, 1, 2, 3) while Latin 
indices will be used for spatial directions (i, j, ... = 1, 2, 3). In many expressions we will use the 
notation = (1, <r l ) where a 1 are the Pauli matrices 




(1.2) 



Sometimes we use of the Feynman's slash notation $ = ^a^. Finally, unless stated otherwise, we 
work in natural units h = c = 1. 



2 Why do we need Quantum Field Theory after all? 

In spite of the impressive success of Quantum Mechanics in describing atomic physics, it was 
immediately clear after its formulation that its relativistic extension was not free of difficulties. 
These problems were clear already to Schrodinger, whose first guess for a wave equation of a free 
relativistic particle was the Klein-Gordon equation 

(^- v2 + m2 )^( t ' f ) = °- ( 2J ) 

This equation follows directly from the relativistic "mass-shell" identity E 2 = p 2 + m 2 using the 
correspondence principle 

p -iV. (2.2) 
Plane wave solutions to the wave equation (2.1) are readily obtained 

ip(t, x) = e~ ip ^ = e - iEt+ip - s with E = ±cjp = ± v / p 2 + m 2 . (2.3) 

In order to have a complete basis of functions, one must include plane wave with both E > and 
E < 0. This implies that given the conserved current 

j, = ^(rd^-d,r^), (2.4) 

its time-component is j° = E and therefore does not define a positive-definite probability density. 
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Fig. 1 : Spectrum of the Klein-Gordon wave equation 



A complete, properly normalized, continuous basis of solutions of the Klein-Gordon equation 
(2.1) labelled by the momentum p can be defined as 



f P (t,x) = 
f- p (t,x) = 



1 



-iuipt+ip-x 



(2tt) 2 v /2^ 
1 



(2tt) 2 v ^ 
Given the inner product 

(ipi\*h) = * J (Fxfadofo - dottfo) 

the states (2.5) form an orthonormal basis 

(f P \fp>) = 5(p-p% 
(f-plf-v?) = S(p-p'), 
(f P \f-P>) = o. 



(2.5) 



(2.6) 
(2.7) 



The wave functions f p (t,x) describes states with momentum p and energy given by u p = 
\lp 2 + m 2 . On the other hand, the states | /_ p ) not only have a negative scalar product but they 
actually correspond to negative energy states 



id f- p (t,x) = -\Jp 2 + m 2 f- p (t,x). 



(2.8) 
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Therefore the energy spectrum of the theory satisfies \E\ > m and is unbounded from below (see 
Fig. 1). Although in a case of a free theory the absence of a ground state is not necessarily a fatal 
problem, once the theory is coupled to the electromagnetic field this is the source of all kinds of 
disasters, since nothing can prevent the decay of any state by emission of electromagnetic radiation. 

The problem of the instability of the "first-quantized" relativistic wave equation can be heuris- 
tically tackled in the case of spin-| particles, described by the Dirac equation 

+ -n?jij(t,x) = 0, (2.9) 

where a and j3 are 4 x 4 matrices 

with a 1 the Pauli matrices, and the wave function tp(t, x) has four components. The wave equation 
(2.9) can be thought of as a kind of "square root" of the Klein-Gordon equation (2.1), since the latter 
can be obtained as 

+ a ■ V - + a • V - ^(f , £ ) = ^ - V 2 + m 2 ) ^(t, £). (2.1 1) 

An analysis of Eq. (2.9) along the lines of the one presented above for the Klein-Gordon 
equation leads again to the existence of negative energy states and a spectrum unbounded from 
below as in Fig. 1. Dirac, however, solved the instability problem by pointing out that now the 
particles are fermions and therefore they are subject to Pauli's exclusion principle. Hence, each 
state in the spectrum can be occupied by at most one particle, so the states with E = m can be made 
stable if we assume that all the negative energy states are filled. 

If Dirac 's idea restores the stability of the spectrum by introducing a stable vacuum where all 
negative energy states are occupied, the so-called Dirac sea, it also leads directly to the conclusion 
that a single-particle interpretation of the Dirac equation is not possible. Indeed, a photon with 
enough energy (E > 2m) can excite one of the electrons filling the negative energy states, leaving 
behind a "hole" in the Dirac see (see Fig. 2). This hole behaves as a particle with equal mass 
and opposite charge that is interpreted as a positron, so there is no escape to the conclusion that 
interactions will produce pairs particle-antiparticle out of the vacuum. 

In spite of the success of the heuristic interpretation of negative energy states in the Dirac 
equation this is not the end of the story. In 1929 Oskar Klein stumbled into an apparent paradox 
when trying to describe the scattering of a relativistic electron by a square potential using Dirac's 
wave equation [11] (for pedagogical reviews see [12, 13]). In order to capture the essence of the 
problem without entering into unnecessary complication we will study Klein's paradox in the con- 
text of the Klein-Gordon equation. 

Let us consider a square potential with height V > of the type showed in Fig. 3. A solution 
to the wave equation in regions I and II is given by 

il)i(t,x) = e - lEt+mx + Re~ iEt - ipiX , 

ipn(t,x) = Te~ iEt+P2X , (2.12) 
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Fig. 2: Creation of a particle-antiparticle pair in the Dirac see picture 
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Fig. 3: Illustration of the Klein paradox. 
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where the mass- shell condition implies that 

pi = \/E 2 - m 2 , p 2 = y/(E - V ) 2 - m 2 . (2.13) 

The constants R and T are computed by matching the two solutions across the boundary x — 0. 
The conditions ^/(^ 0) = tpn(t, 0) and d x ipi(t, 0) = d x ?pjj(t, 0) imply that 

T=^^, i?=^^. (2.14) 

Pi +P2 Pl+ P2 

At first sight one would expect a behavior similar to the one encountered in the nonrelativistic 
case. If the kinetic energy is bigger than V both a transmitted and reflected wave are expected, 
whereas when the kinetic energy is smaller than V one only expect to find a reflected wave, the 
transmitted wave being exponentially damped within a distance of a Compton wavelength inside 
the barrier. 

Indeed this is what happens if E — m > V . In this case both p 1 and p 2 are real and we have a 
partly reflected, and a partly transmitted wave. In the same way, if E—m < V and E—m < V —2m 
then p 2 is imaginary and there is total reflection. 

However, in the case when V > 2m and the energy is in the range V — 2m < E — m < V a 
completely different situation arises. In this case one finds that both p 1 and p 2 are real and therefore 
the incoming wave function is partially reflected and partially transmitted across the barrier. This is 
a shocking result, since it implies that there is a nonvanishing probability of finding the particle at 
any point across the barrier with negative kinetic energy (E — m — V < 0)1 This weird result is 
known as Klein's paradox. 

As with the negative energy states, the Klein paradox results from our insistence in giving a 
single-particle interpretation to the relativistic wave function. Actually, a multiparticle analysis of 
the paradox [12] shows that what happens when E — m > V — 2m is that the reflection of the 
incoming particle by the barrier is accompanied by the creation of pairs particle-antiparticle out of 
the energy of the barrier (notice that for this to happen it is required that V > 2m, the threshold for 
the creation of a particle-antiparticle pair). 

Actually, this particle creation can be understood by noticing that the sudden potential step in 
Fig. 3 localizes the incoming particle with mass m in distances smaller than its Compton wavelength 
A = — . This can be seen by replacing the square potential by another one where the potential varies 
smoothly from to Vq > 2m in distances scales larger than 1/m. This case was worked out by 
Sauter shortly after Klein pointed out the paradox [14]. He considered a situation where the regions 
with V = and V = V axe connected by a region of length d with a linear potential V(x) = ^p. 
When d > ^ he found that the transmission coefficient is exponentially small 1 . 

The creation of particles is impossible to avoid whenever one tries to locate a particle of mass 
m within its Compton wavelength. Indeed, from Heisenberg uncertainty relation we find that if 
Ax ~ ^, the fluctuations in the momentum will be of order Ap ~ m and fluctuations in the energy 
of order 

AE ~ m (2.15) 



'in section (9.1) we will see how, in the case of the Dirac field, this exponential behavior can be associated with the 
creation of electron-positron pairs due to a constant electric field (Schwinger effect). 
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Fig. 4: Two regions R\, R 2 that are causally disconnected. 



can be expected. Therefore, in a relativistic theory, the fluctuations of the energy are enough to 
allow the creation of particles out of the vacuum. In the case of a spin-| particle, the Dirac sea 
picture shows clearly how, when the energy fluctuations are of order m, electrons from the Dirac 
sea can be excited to positive energy states, thus creating electron-positron pairs. 

It is possible to see how the multiparticle interpretation is forced upon us by relativistic invari- 
ance. In non-relativistic Quantum Mechanics observables are represented by self-adjoint operator 
that in the Heisenberg picture depend on time. Therefore measurements are localized in time but 
are global in space. The situation is radically different in the relativistic case. Because no signal 
can propagate faster than the speed of light, measurements have to be localized both in time and 
space. Causality demands then that two measurements carried out in causally-disconnected regions 
of space-time cannot interfere with each other. In mathematical terms this means that if Rl and 
Or 2 are the observables associated with two measurements localized in two causally-disconnected 
regions R±, R 2 (see Fig. 4), they satisfy 

[O Rl ,O R2 } = 0, if (xi — X2) 2 < 0, for all x\ G R±, x 2 G R 2 . (2.16) 

Hence, in a relativistic theory, the basic operators in the Heisenberg picture must depend on 
the space-time position x* 1 . Unlike the case in non-relativistic quantum mechanics, here the position 
x is not an observable, but just a label, similarly to the case of time in ordinary quantum mechanics. 
Causality is then imposed microscopically by requiring 

[O(x),O(y)] = 0, if(x-y) 2 <0. (2.17) 

A smeared operator O r over a space-time region R can then be defined as 

O r = [ d A xO{x)f R (x) (2.18) 
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where Jr(x) is the characteristic function associated with R, 



/«(«) = { J X X %1- (2-19) 

Eq. (2.16) follows now from the microcausality condition (2.17). 

Therefore, relativistic invariance forces the introduction of quantum fields. It is only when 
we insist in keeping a single-particle interpretation that we crash against causality violations. To 
illustrate the point, let us consider a single particle wave function ip(t, x) that initially is localized 
in the position x = 

ip(0,x) = 5(x). (2.20) 



Evolving this wave function using the Hamiltonian H = V— V 2 + m 2 we find that the wave func- 
tion can be written as 

i/j(t,x) = e~ it ^ WT ^5(x) = [ ^L e ik-s-uV^+^_ (2 .21) 



., (2vr) 3 

Integrating over the angular variables, the wave function can be recast in the form 

1 f°° 

^(t,x) = —— kdke ik ^ e ~^ k2 + m \ (2.22) 

Z7T \X\ J_ 00 

The resulting integral can be evaluated using the complex integration contour C shown in Fig. 5. 
The result is that, for any t > 0, one finds that ip(t, x) ^ for any x. If we insist in interpreting the 
wave function ip(t, x) as the probability density of finding the particle at the location x in the time t 
we find that the probability leaks out of the light cone, thus violating causality. 



3 From classical to quantum fields 

We have learned how the consistency of quantum mechanics with special relativity forces us to 
abandon the single-particle interpretation of the wave function. Instead we have to consider quantum 
fields whose elementary excitations are associated with particle states, as we will see below. 

In any scattering experiment, the only information available to us is the set of quantum number 
associated with the set of free particles in the initial and final states. Ignoring for the moment other 
quantum numbers like spin and flavor, one-particle states are labelled by the three-momentum p and 
span the single-particle Hilbert space "Hi 

\p)eHu (p\p') = 8(p-p')- (3-D 
The states {\p}} form a basis of Hi and therefore satisfy the closure relation 

J d 3 p\p)(p\ = l (3.2) 
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Fig. 5: Complex contour C for the computation of the integral in Eq. (2.22). 

The group of spatial rotations acts unitarily on the states \p). This means that for every rotation 
R E SO(3) there is a unitary operator U(R) such that 

U(R)\p) = \Rp) (3.3) 

where Rp represents the action of the rotation on the vector k, (Rp) 2 = R 1 ,jkK Using a spectral 
decomposition, the momentum operator P l can be written as 




(3.4) 



With the help of Eq. (3.3) it is straightforward to check that the momentum operator transforms as 
a vector under rotations: 

U(R)- 1 ptUiR) = J d 3 p \R~ l p) p i (R^pl = R^pl, (3.5) 

where we have used that the integration measure is invariant under SO(3). 

Since, as we argued above, we are forced to deal with multiparticle states, it is convenient to 
introduce creation-annihilation operators associated with a single-particle state of momentum p 

Hf), a\p')} = 5(p- p'), [a(p), a(p')} = \a\p), a\p')} = 0, (3.6) 

such that the state \p) is created out of the Fock space vacuum |0) (normalized such that (0 1 0) = 1) 
by the action of a creation operator a* (p) 

\p) = <f(p)\0), a(p)|0) = Vp. (3.7) 
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Covariance under spatial rotations is all we need if we are interested in a nonrelativistic theory. 
However in a relativistic quantum field theory we must preserve more that SO(3), actually we 
need the expressions to be covariant under the full Poincare group ISO(l, 3) consisting in spatial 
rotations, boosts and space-time translations. Therefore, in order to build the Fock space of the 
theory we need two key ingredients: first an invariant normalization for the states, since we want a 
normalized state in one reference frame to be normalized in any other inertial frame. And secondly 
a relativistic invariant integration measure in momentum space, so the spectral decomposition of 
operators is covariant under the full Poincare group. 

Let us begin with the invariant measure. Given an invariant function f(p) of the four-momen- 
tum p^ of a particle of mass m with positive energy p° > 0, there is an integration measure which 
is invariant under proper Lorentz transformations 2 

/ ^( 2 ^(P 2 -™W)/(P), (3-8) 

where 9{x) represent the Heaviside step function. The integration over p° can be easily done using 
the 5-function identity 

*[/(*)]= E wLx\ 6{x - Xi ^ (3 - 9) 

Zi=zeros of / 

which in our case implies that 

S(p 2 - m 2 ) = _L 5 (p° - + m 2 j + ±_ S (p°+ + m 2 ) . (3.10) 

The second term in the previous expression correspond to states with negative energy and therefore 
does not contribute to the integral. We can write then 

Hence, the relativistic invariant measure is given by 

S§fk wi,h < 3 - 12 ' 

Once we have an invariant measure the next step is to find an invariant normalization for the 
states. We work with a basis {\p)} of eigenstates of the four-momentum operator 

P°\p)=cu p \p), P 1 \P)=P 1 \P)- (3-13) 

Since the states \p) are eigenstates of the three-momentum operator we can express them in terms 
of the non-relativistic states \p) that we introduced in Eq. (3.1) 

\p) = N(p)\p) (3.14) 



2 The factors of 2n are introduced for later convenience. 
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with N{p) a normalization to be determined now. The states {\p)} form a complete basis, so they 
should satisfy the Lorentz invariant closure relation 

/ (0* {2n)5{p2 - m2) 9{p0) |P> {pl = 1 (3 - 15) 

At the same time, this closure relation can be expressed, using Eq. (3.14), in terms of the nonrela- 
tivistic basis of states {\p}} as 

/ (0* {27T)6{p2 ~ m2) 6{p0) lp) {pl = I j^ij~J N{p)l2 lp) ^ (3 - 16) 

Using now Eq. (3.4) for the nonrelativistic states, expression (3.15) follows provided 

\N(p)\ 2 = (2n) 3 (2co p ). (3.17) 

Taking the overall phase in Eq. (3.14) so that N(p) is real, we define the Lorentz invariant states \p) 
as 

|p) = (2 7 r)i v ^|p), (3.18) 
and given the normalization of \p) we find the normalization of the relativistic states to be 

(p\p') = (2ixf(2uj p )5(p-p'). (3.19) 

Although not obvious at first sight, the previous normalization is Lorentz invariant. Although 
it is not difficult to show this in general, here we consider the simpler case of 1+1 dimensions where 
the two components (p ,]? 1 ) of the on-shell momentum can be parametrized in terms of a single 
hyperbolic angle A as 

p° = mcoshA, p 1 = rasinhA. (3.20) 

Now, the combination 2u p 5(p 1 — p 1 ') can be written as 

2u p S(p 1 - p 1 ') = 2mcoshA<5(msinhA -msinhA') = 25 (X - A'), (3.21) 

where we have made use of the property (3.9) of the 5-function. Lorentz transformations in 1 + 1 
dimensions are labelled by a parameter £ e K and act on the momentum by shifting the hyperbolic 
angle A — > X + £. However, Eq. (3.21) is invariant under a common shift of A and A', so the whole 
expression is obviously invariant under Lorentz transformations. 

To summarize what we did so far, we have succeed in constructing a Lorentz covariant basis 
of states for the one-particle Hilbert space Hi. The generators of the Poincare group act on the 
states \p) of the basis as 

^P»=P», W(A)|p) = lA^p") = \Ap) with AeSO(l,3). (3.22) 
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This is compatible with the Lorentz invariance of the normalization that we have checked above 

W) = (v\U(k)- l U^)\p) = (Ap|Ap')- (3-23) 
On Hi the operator P M admits the following spectral representation 

^ = /S^"W' <3 - 24) 

Using (3.23) and the fact that the measure is invariant under Lorentz transformation, one can easily 
show that P M transform covariantly under SO(l, 3) 

W(A)- 1 P^(A) = J ^-L-IA- 1 ?)^ (A" VI = A^P". (3.25) 

A set of covariant creation- annihilation operators can be constructed now in terms of the 
operators a(p), at(p) introduced above 

a(p) = (27r)i v ^a(p), a\p) = {2tt)I y^a^p) (3.26) 

with the Lorentz invariant commutation relations 

[a(p),a(p')} = [a\p),a\p')]=0. (3.27) 

Particle states are created by acting with any number of creation operators a(p) on the Poincare 
invariant vacuum state |0) satisfying 



(0|0) = 1, P"|0) = 0, W(A)|0) = |0), VAeSO(l,3). (3.28) 
A general one-particle state |/) e Hi can be then written as 

d 3 p 1 

j p 

while a n-particle state |/) G Hf n can be expressed as 



l/> = /n^^/(Pi.---.P> , (Pi)^^' , (P»)|0>^ (3.30) 

That this states are Lorentz invariant can be checked by noticing that from the definition of the 
creation- annihilation operators follows the transformation 

U(A)a(p)U(Ay = a(Ap) (3.31) 

and the corresponding one for creation operators. 

As we have argued above, the very fact that measurements have to be localized implies the 
necessity of introducing quantum fields. Here we will consider the simplest case of a scalar quantum 
field (j)(x) satisfying the following properties: 
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- Hermiticity. 

<fi(x) = <t>{x). (3.32) 

- Microcausality. Since measurements cannot interfere with each other when performed in 
causally disconnected points of space-time, the commutator of two fields have to vanish out- 
side the relative ligth-cone 

[0(z),0(y)] = O, (x-y) 2 <0. (3.33) 

- Translation invariance. 

e ip - a <P(x)e~ tp - a = <j>{x - a). (3.34) 

- Lorentz invariance. 

W(A)V(a:)W(A) = 0(A~V). (3.35) 

- Linearity. To simplify matters we will also assume that 4>(x) is linear in the creation- 
annihilation operators a(p), a\p) 

H x ) = J (2^3 2uj x ) a (P) + 9(P, x)a\p)] . (3.36) 

Since <f)(x) should be hermitian we are forced to take f(p,x)* = g(p,x). Moreover, 0(x) 
satisfies the equations of motion of a free scalar field, (d^d^ + m 2 )(p(x) = 0, only if f(p, x) 
is a complete basis of solutions of the Klein-Gordon equation. These considerations leads to 
the expansion 

^ = J (1^277 [ e ~ Wpt+ ^ a (P) + e^-^a^p)} . (3.37) 

Given the expansion of the scalar field in terms of the creation- annihilation operators it can be 
checked that 4>(x) and d t (j)(x) satisfy the equal-time canonical commutation relations 

[<f>(t,x),d t <f>(t,$)}=i5(x-$) (3.38) 

The general commutator [4>(x), 4>{y)} can be also computed to be 

[<f)(x), <f)(x')} = iA(x - x'). (3.39) 

The function A(x — y) is given by 

iA(x-v) = -Im I L e -<wp(*-* , )+*?-(*-*') 

= J^(2n)5(p 2 -rn 2 )e(p°)e-^'\ (3-40) 
15 



where e(x) is defined as 



e(x) = 9(x) - 9(-x) = 



1 x > 
-1 x < 



(3.41) 



Using the last expression in Eq. (3.40) it is easy to show that iA(x — x') vanishes when x 
and x' are space-like separated. Indeed, if (x — x') 2 < there is always a reference frame in which 
both events are simultaneous, and since iA(x — x') is Lorentz invariant we can compute it in this 
reference frame. In this case t — t' and the exponential in the second line of (3.40) does not depend 
on p°. Therefore, the integration over k° gives 



So we have concluded that iA(x — x') — if (x — x') 2 < 0, as required by microcausality. Notice 
that the situation is completely different when (x — x') 2 > 0, since in this case the exponential 
depends on p° and the integration over this component of the momentum does not vanish. 

3.1 Canonical quantization 

So far we have contented ourselves with requiring a number of properties to the quantum scalar field: 
existence of asymptotic states, locality, microcausality and relativistic invariance. With these only 
ingredients we have managed to go quite far. The previous can also be obtained using canonical 
quantization. One starts with a classical free scalar field theory in Hamiltonian formalism and 
obtains the quantum theory by replacing Poisson brackets by commutators. Since this quantization 
procedure is based on the use of the canonical formalism, which gives time a privileged role, it 
is important to check at the end of the calculation that the resulting quantum theory is Lorentz 
invariant. In the following we will briefly overview the canonical quantization of the Klein-Gordon 
scalar field. 

The starting point is the action functional S[<f>(x)] which, in the case of a free real scalar field 
of mass m is given by 




(3.42) 




d A x (d^dy - m 2 (j) 2 ) . 



(3.43) 



The equations of motion are obtained, as usual, from the Euler-Lagrange equations 



dC ] dC 




(3.44) 



The momentum canonically conjugated to the field <j)(x) is given by 

. _ dC d(j) 



(3.45) 



d(d (j)) dt ' 
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In the Hamiltonian formalism the physical system is described not in terms of the generalized coor- 
dinates and their time derivatives but in terms of the generalized coordinates and their canonically 
conjugated momenta. This is achieved by a Legendre transformation after which the dynamics of 
the system is determined by the Hamiltonian function 



H 



d 3 x [n^-C 



I / d 3 x 



7T 2 + (V0) 2 + m 5 



(3.46) 



The equations of motion can be written in terms of the Poisson rackets. Given two functional 

A[(f), 7r], B[6, it] of the canonical variables 



4<M 



d 3 xA(6,n), 



B 



[<M = J 



d 3 xB(6,n). 



Their Poisson bracket is defined by 

{A, B} = J d 3 x 
where 4i denotes the functional derivative defined as 



5 A 5B 5 A 5B 
56 5n 5tt 56 



(3.47) 



(3.48) 



5A dA 



d„ 



dA 



56 86 11 [d(d^6)_ 
Then, the canonically conjugated fields satisfy the following equal time Poisson brackets 



(3.49) 



{6(t,x),6(t,x' )} 
{6(t,x),n(t,x')} 



{n(t,x),n(t,x')} = 0, 
5(x — x'). 



(3.50) 



Canonical quantization proceeds now by replacing classical fields with operators and Poisson 
brackets with commutators according to the rule 



*{•>•} 



(3.51) 



In the case of the scalar field, a general solution of the field equations (3.44) can be obtained by 
working with the Fourier transform 



(d^ + m 2 )6(x) = 
whose general solution can be written as 3 

d 4 p 



(-p 2 + m 2 )6(p) = 0, 



(3.52) 



6{x) = 



(2n)5(p 2 - m 2 )6(p°) [a{p)e~ ip - x + a(p)*e ip - x ] 
d 3 p 1 



(27T) 4 



(2vr) 3 2u p 



[a(p)e- iUpt+ P- s + a(p)*e luJpt -^] 



(3.53) 



3 In momentum space, the general solution to this equation is cj>(p) = f(p)5(p 2 — to 2 ), with f(p) a completely 
general function of The solution in position space is obtained by inverse Fourier transform. 
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and we have required <f>(x) to be real. The conjugate momentum is 

= ~J [a(p)e~^ t+ ^ + atfye***-**] . (3.54) 

Now (j)(x) and n(x) are promoted to operators by replacing the functions a(p), a(p)* by the 
corresponding operators 

a(p) — > a(p), a(p)* — > a j (p). (3.55) 

Moreover, demanding [<j>(t, x), ir(t, x')] = i5(x — x') forces the operators a(p), a(p)^ to have 
the commutation relations found in Eq. (3.27). Therefore they are identified as a set of creation- 
annihilation operators creating states with well-defined momentum p out of the vacuum |0). In the 
canonical quantization formalism the concept of particle appears as a result of the quantization of a 
classical field. 

Knowing the expressions of <\> and n in terms of the creation- annihilation operators we can 
proceed to evaluate the Hamiltonian operator. After a simple calculation one arrives to the expres- 
sion 



H = d 3 p 



u p a\p)a{p) + \u P 5(0) 



(3.56) 



The first term has a simple physical interpretation since a*(p)a(p) is the number operator of par- 
ticles with momentum p. The second divergent term can be eliminated if we defined the normal- 
ordered Hamiltonian :H: with the vacuum energy subtracted 

:H:= H - (0\H\0) = J <Ppu p a\p)a(p) (3.57) 

It is interesting to try to make sense of the divergent term in Eq. (3.56). This term have two 
sources of divergence. One is associated with the delta function evaluated at zero coming from the 
fact that we are working in a infinite volume. It can be regularized for large but finite volume by 
replacing 5(0) ~ V. Hence, it is of infrared origin. The second one comes from the integration of 
lu p at large values of the momentum and it is then an ultraviolet divergence. The infrared divergence 
can be regularized by considering the scalar field to be living in a box of finite volume V. In this 
case the vacuum energy is 

E vac = (0\H\0) = Y,\u p . (3.58) 

p 

Written in this way the interpretation of the vacuum energy is straightforward. A free scalar quantum 
field can be seen as a infinite collection of harmonic oscillators per unit volume, each one labelled 
by p. Even if those oscillators are not excited, they contribute to the vacuum energy with their zero- 
point energy, given by \<jJ v - This vacuum contribution to the energy add up to infinity even if we 
work at finite volume, since even then there are modes with arbitrary high momentum contributing 
to the sum, p« = ! r 2E , with Li the sides of the box of volume V and an integer. Hence, this 
divergence is of ultraviolet origin. 
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Region I 



Region II 



Region III 



Conducting plates 

Fig. 6: Illustration of the Casimir effect. In regions I and II the spetrum of modes of the momentum p± is 
continuous, while in the space between the plates (region II) it is quantized in units of ^. 



3.2 The Casimir effect 

The presence of a vacuum energy is not characteristic of the scalar field. It is also present in other 
cases, in particular in quantum electrodynamics. Although one might be tempted to discarding this 
infinite contribution to the energy of the vacuum as unphysical, it has observable consequences. In 
1948 Hendrik Casimir pointed out [15] that although a formally divergent vacuum energy would 
not be observable, any variation in this energy would be (see [16] for comprehensive reviews). 

To show this he devised the following experiment. Consider a couple of infinite, perfectly 
conducting plates placed parallel to each other at a distance d (see Fig. 6). Because the conducting 
plates fix the boundary condition of the vacuum modes of the electromagnetic field these are discrete 
in between the plates (region II), while outside there is a continuous spectrum of modes (regions 
I and III). In order to calculate the force between the plates we can take the vacuum energy of 
the electromagnetic field as given by the contribution of two scalar fields corresponding to the two 
polarizations of the photon. Therefore we can use the formulas derived above. 

A naive calculation of the vacuum energy in this system gives a divergent result. This infinity 
can be removed, however, by substracting the vacuum energy corresponding to the situation where 
the plates are removed 



E(d) ies = E(d) vac - E(oo) 



(3.59) 



This substraction cancels the contribution of the modes outside the plates. Because of the bound- 
ary conditions imposed by the plates the momentum of the modes perpendicular to the plates are 
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quantized according to p± = B f, with n a non-negative integer. If we consider that the size of the 
plates is much larger than their separation d we can take the momenta parallel to the plates p\\ as 
continuous. For n > we have two polarizations for each vacuum mode of the electromagnetic 

field, each contributing like \ y^p| + v\ t0 the vacuum energy. On the other hand, when p± = the 

corresponding modes of the field are effectively (2+l)-dimensional and therefore there is only one 
polarization. Keeping this in mind, we can write 



2Sd I tI^IpI (3.60) 



(2tt) 3 2' 

where S is the area of the plates. The factors of 2 take into account the two propagating degrees 
of freedom of the electromagnetic field, as discussed above. In order to ensure the convergence of 
integrals and infinite sums we can introduce an exponential damping factor 4 

where A is an ultraviolet cutoff. It is now straightforward to see that if we define the function 

m = I [yiy^^ff^f = -L dze ~4^ (3 .62) 

the regularized vacuum energy can be written as 



I 00 PCX) 

-F(0) + 5>W-/ dxF( 

n=l ^° 



X) 



E(d) rcg — S 

This expression can be evaluated using the Euler-MacLaurin formula [18] 



(3.63) 



00 

n=l 



J~dxF(x) = -I[F(0) + F(oo)] + l[F'(oo)-nO)] 

- ^ [^"(00) - F"'(0)] + . . . (3.64) 

Since for our function F{oo) = F'(oo) = F"'(oo) = and F'(0) = 0, the value of E(d) reg is 
determined by F"'(0). Computing this term and removing the ultraviolet cutoff, A — > 00 we find 
the result 

E(d), v = 4f-(0) = (3.65) 



4 Actually, one could introduce any cutoff function f{p\ going to zero f ast enough as p±, — > 00. The result 
is independent of the particular function used in the calculation. 
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Then, the force per unit area between the plates is given by 



iW = ~240 (3 - 66) 

The minus sign shows that the force between the plates is attractive. This is the so-called Casimir 
effect. It was experimentally measured in 1958 by Sparnaay [17] and since then the Casimir effect 
has been checked with better and better precission in a variety of situations [16]. 

4 Theories and Lagrangians 

Up to this point we have used a scalar field to illustrate our discussion of the quantization procedure. 
However, nature is richer than that and it is necessary to consider other fields with more complicated 
behavior under Lorentz transformations. Before considering other fields we pause and study the 
properties of the Lorentz group. 

4.1 Representations of the Lorentz group 

In four dimensions the Lorentz group has six generators. Three of them correspond to the generators 
of the group of rotations in three dimensions SO(3). In terms of the generators Jj of the group a 
finite rotation of angle if with respect to an axis determined by a unitary vector e can be written as 



R(e, ip) = 




(4.1) 



The other three generators of the Lorentz group are associated with boosts Mj along the three spatial 
directions. A boost with rapidity A along a direction u is given by 



B(u, A) = 




(4.2) 



These six generators satisfy the algebra 

[Ji,M k ] = ie ljk M k) (4.3) 
[Mi,Mj} = -ie ijk J k , 

The first line corresponds to the commutation relations of SO(3), while the second one implies that 
the generators of the boosts transform like a vector under rotations. 

At first sight, to find representations of the algebra (4.3) might seem difficult. The problem is 
greatly simplified if we consider the following combination of the generators 

J± = ±(J k ±iM k ). (4.4) 
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Representation 


Type of field 


(0,0) 


Scalar 


(1,0) 


Right-handed spinor 


(0,1) 


Left-handed spinor 


f 1 1 ) 

V2' 2' 


Vector 


(1,0) 


Selfdual antisymmetric 2-tensor 


(0,1) 


Anti-selfdual antisymmetric 2-tensor 



Table 1 : Representations of the Lorentz group 



Using (4.3) it is easy to prove that the new generators Jf satisfy the algebra 

[J^iJf] = ieijk-Jk : , 

[J+,J-] = 0. (4.5) 

Then the Lorentz algebra (4.3) is actually equivalent to two copies of the algebra of SU(2) « SO (3). 
Therefore the irreducible representations of the Lorentz group can be obtained from the well-known 
representations of SU(2). Since the latter ones are labelled by the spin s = k + |, k (with k € N), 
any representation of the Lorentz algebra can be identified by specifying (s + , s_), the spins of the 
representations of the two copies of SU(2) that made up the algebra (4.3). 

To get familiar with this way of labelling the representations of the Lorentz group we study 
some particular examples. Let us start with the simplest one (s + , s_) = (0, 0). This state is a singlet 
under jf and therefore also under rotations and boosts. Therefore we have a scalar. 

The next interesting cases are (|, 0) and (0, |). They correspond respectively to a right- 
handed and a left-handed Weyl spinor. Their properties will be studied in more detail below. In 
the case of (1,1), since from Eq. (4.4) we see that J, = Jf + J~ the rules of addition of angular 
momentum tell us that there are two states, one of them transforming as a vector and another one as 
a scalar under three-dimensional rotations. Actually, a more detailed analysis shows that the singlet 
state corresponds to the time component of a vector and the states combine to form a vector under 
the Lorentz group. 

There are also more "exotic" representations. For example we can consider the (1, 0) and 
(0, 1) representations corresponding respectively to a selfdual and an anti-selfdual rank-two anti- 
symmetric tensor. In Table 1 we summarize the previous discussion. 

To conclude our discussion of the representations of the Lorentz group we notice that under a 
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parity transformation the generators of SO(l,3) transform as 

P-.Ji^Ji, P:M % ^ -Mi (4.6) 

this means that P : jf — > and therefore a representation (si, s 2 ) is transformed into (s 2 , Si). 
This means that, for example, a vector (|, |) is invariant under parity, whereas a left-handed Weyl 
spinor (|, 0) transforms into a right-handed one (0, |) and vice versa. 

4.2 Spinors 

Weyl spinors. Let us go back to the two spinor representations of the Lorentz group, namely (§, 0) 
and (0, |). These representations can be explicitly constructed using the Pauli matrices as 

J+ = \a\ J-=0 for (1,0), 

We denote by u± a complex two-component object that transforms in the representation s± = \ of 
J l ± . If we define a± = (1, ±<r l ) we can construct the following vector quantities 

v} + a+u +) uLalu-. (4.8) 

Notice that since ( jfY = the hermitian conjugated fields v)± are in the (0, |) and (|, 0) respec- 
tively. 

To construct a free Lagrangian for the fields u± we have to look for quadratic combinations 
of the fields that are Lorentz scalars. If we also demand invariance under global phase rotations 

u± — > e ie u± (4.9) 

we are left with just one possibility up to a sign 

£ Wey i = i«± (d t ±3- V) u ± = iul<j^u±. (4.10) 

This is the Weyl Lagrangian. In order to grasp the physical meaning of the spinors u± we write the 
equations of motion 

(db±a-v)u± = 0. (4-11) 

Multiplying this equation on the left by (^d =F a ■ V j and applying the algebraic properties of the 
Pauli matrices we conclude that u± satisfies the massless Klein-Gordon equation 

d^u ± = 0, (4.12) 

whose solutions are: 

u±(x)=u±(k)e- ik ' x , with fc° = |jfe|. (4.13) 
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Plugging these solutions back into the equations of motion (4.11) we find 



(\k\^k.ff y )u ± = 0, (4.14) 

which implies 

a ■ k 

u+ : = 1, 

\k\ 

— * 

U- : = -1. (4.15) 

| A; | 

Since the spin operator is defined as s = \a, the previous expressions give the chirality of the states 
with wave function u±, i.e. the projection of spin along the momentum of the particle. Therefore 
we conclude that u + is a Weyl spinor of positive helicity A = |, while w_ has negative helicity 
A = — \. This agrees with our assertion that the representation (|, 0) corresponds to a right-handed 
Weyl fermion (positive chirality) whereas (0, \) is a left-handed Weyl fermion (negative chirality). 
For example, in the Standard Model neutrinos are left-handed Weyl spinors and therefore transform 
in the representation (0, |) of the Lorentz group. 

Nevertheless, it is possible that we were too restrictive in constructing the Weyl Lagrangian 
(4.10). There we constructed the invariants from the vector bilinears (4.8) corresponding to the 
product representations 

(|, |) = (|,0)® (0,|) and (|, |) = (0,|)® (|,0). (4.16) 

In particular our insistence in demanding the Lagrangian to be invariant under the global symmetry 
u± — > e t9 u± rules out the scalar term that appears in the product representations 

(|, 0)®(|,0) = (1, 0) © (0, 0), (0, §) ® (0, §) = (0, 1) © (0, 0). (4.17) 

The singlet representations corresponds to the antisymmetric combinations 

e ab u a ± u b ± , (4.18) 

where e ab is the antisymmetric symbol e 12 = — e 2 i = 1. 

At first sight it might seem that the term (4.18) vanishes identically because of the antisym- 
metry of the e-symbol. However we should keep in mind that the spin-statistic theorem (more on 
this later) demands that fields with half-integer spin have to satisfy the Fermi-Dirac statistics and 
therefore satisfy anticommutation relations, whereas fields of integer spin follow the statistic of 
Bose-Einstein and, as a consequence, quantization replaces Poisson brackets by commutators. This 
implies that the components of the Weyl fermions u± are anticommuting Grassmann fields 

u a ± u b ± + u b ± u a ± = 0. (4.19) 

It is important to realize that, strictly speaking, fermions (i.e., objects that satisfy the Fermi-Dirac 
statistics) do not exist classically. The reason is that they satisfy the Pauli exclusion principle and 
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therefore each quantum state can be occupied, at most, by one fermion. Therefore the naive defini- 
tion of the classical limit as a limit of large occupation numbers cannot be applied. Fermion field 
do not really make sense classically. 

Since the combination (4.18) does not vanish and we can construct a new Lagrangian 

^weyi = i«±<T±9 M «± + ^e«6<«± + h.c. (4.20) 

This mass term, called of Majorana type, is allowed if we do not worry about breaking the global 
U(l) symmetry u± — > e %6 u±. This is not the case, for example, of charged chiral fermions, since the 
Majorana mass violates the conservation of electric charge or any other gauge U(l) charge. In the 
Standard Model, however, there is no such a problem if we introduce Majorana masses for right- 
handed neutrinos, since they are singlet under all standard model gauge groups. Such a term will 
break, however, the global U(l) lepton number charge because the operator e a bi'%i' b R changes the 
lepton number by two units 

Dirac spinors. We have seen that parity interchanges the representations (|, 0) and (0, |), 
i.e. it changes right-handed with left-handed fermions 

P : u± — > u T . (4.21) 

An obvious way to build a parity invariant theory is to introduce a pair or Weyl fermions u + and u + . 
Actually, these two fields can be combined in a single four-component spinor 

*=(::) < 4 - 22 > 

transforming in the reducible representation (|, 0) © (0, |). 

Since now we have both u + and u_ simultaneously at our disposal the equations of motion 
for u±, ia±d ll u± = can be modified, while keeping them linear, to 



°+ jj, W = m( J J )r. (4.23) 



ta^d^U- = mu + 

These equations of motion can be derived from the Lagrangian density 

£ D irac = # f ( ^ i ) - ( ? 1) V>. (4.24) 



at ) ^ r '" r V 1 
To simplify the notation it is useful to define the Dirac 7-matrices as 



and the Dirac conjugate spinor ip 



i> = = ^ I J J V (4.26) 
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Now the Lagrangian (4.24) can be written in the more compact form 

£Dirac = ^7^-™)^ ( 4 - 27 ) 

The associated equations of motion give the Dirac equation (2.9) with the identifications 

7° = p, 7* = ia\ (4.28) 

In addition, the 7-matrices defined in (4.25) satisfy the Clifford algebra 

{/,f} = 2f. (4-29) 

In D dimensions this algebra admits representations of dimension 2^. When D is even the Dirac 
fermions ip transform in a reducible representation of the Lorentz group. In the case of interest, 
D = 4 this is easy to prove by defining the matrix 

7 5 = -iVVrV = ( I (4.30) 

We see that 7 5 anticommutes with all other 7-matrices. This implies that 

[y>^] = o, with a^ = -^[Y,Y]- (4.31) 

Because of Schur's lemma (see Appendix) this implies that the representation of the Lorentz group 
provided by is reducible into subspaces spanned by the eigenvectors of 7 s with the same eigen- 
value. If we define the projectors P± — |(1 ± 7 s ) these subspaces correspond to 

^ "-*={!!-)• (432) 

which are precisely the Weyl spinors introduced before. 

Our next task is to quantize the Dirac Lagrangian. This will be done along the lines used for 
the Klein-Gordon field, starting with a general solution to the Dirac equation and introducing the 
corresponding set of creation-annihilation operators. Therefore we start by looking for a complete 
basis of solutions to the Dirac equation. In the case of the scalar field the elements of the basis were 
labelled by their four- momentum k^. Now, however, we have more degrees of freedom since we 
are dealing with a spinor which means that we have to add extra labels. Looking back at Eq. (4.15) 
we can define the helicity operator for a Dirac spinor as 

A=i?.i^?V (4.33) 



k\\0 1 

Hence, each element of the basis of functions is labelled by its four-momentum k^ and the corre- 
sponding eigenvalue s of the helicity operator. For positive energy solutions we then propose the 
ansatz 

u{k,s)e- ik - x , s = ± l - 1 (4.34) 
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where u a (k, s) (a = 1, . . . , 4) is a four-component spinor. Substituting in the Dirac equation we 
obtain 

Qk-m)u(k,s) = 0. (4.35) 
In the same way, for negative energy solutions we have 

v(k,s)e ik - x , s = ± 1 -, (4.36) 

where v(k, s) has to satisfy 

($ + m)v(k,s)=0. (4.37) 

Multiplying Eqs. (4.35) and (4.37) on the left respectively by =f m) we find that the momentum 
is on the mass shell, k 2 = m 2 . Because of this, the wave function for both positive- and negative- 

— * — * — * 

energy solutions can be labeled as well using the three-momentum k of the particle, u(k, s), v(k, s). 
A detailed analysis shows that the functions u(k, s), v(k, s) satisfy the properties 

u(k, s)u(k, s) = 2m, v(k, s)v(k, s) = —2m, 

u(k, s)yu(jfe, s) = 2A/\ v(k, s)^v(k, s) = 2P, (4.38) 

^2 u a(k, s)up{k, s) = + m) a p, ^2 u «(*> s )vp$, s ) = - m W' 

*=±s s =±h 

with k° = oo k = Vk 2 + rn 2 . Then, a general solution to the Dirac equation including creation and 
annihilation operators can be written as: 

= y T^p^ E [u(k, 8 )b(k,8)e^™* + v(k,8)$(k, 8^*^*1 . (4.39) 



The operators b\(k, s), b a (k) respectively create and annihilate a spin-| particle (for example, 
an electron) out of the vacuum with momentum k and helicity s. Because we are dealing with 
half-integer spin fields, the spin-statistics theorem forces canonical anticommutation relations for ip 
which means that the creation-annihilation operators satisfy the algebra 5 

{b a (k,s),b^(k',s')} = S(k - k ')5 a/3 5 ss >, 

{b a (k,s),b p (k',s')} = {&t(M), a')} = 0. (4.40) 

In the case of d a (k, s), d\(k, s) we have a set of creation-annihilation operators for the corre- 
sponding antiparticles (for example positrons). This is clear if we notice that d\{k, s) can be seen 
as the annihilation operator of a negative energy state of the Dirac equation with wave function 

5 To simplify notation, and since there is no risk of confusion, we drop from now on the hat to indicate operators. 
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v a (k, s). As we saw, in the Dirac sea picture this corresponds to the creation of an antiparticle out 
of the vacuum (see Fig. 2). The creation-annihilation operators for antiparticles also satisfy the 
fermionic algebra 



{d a (k,s),(Vp(k',s')} = 6(k - k ')5 a/3 5 ss/ , 
{d a (k,s),dp(k',s')} = {dl(k,s),dl(k',s')} 



0. 



(4.41) 



— * , — * — * — * 

All other anticommutators between b a (k, s), b* a (k, s) and d a (k, s), d* a (k, s) vanish. 
The Hamiltonian operator for the Dirac field is 

H = ^2 / ^k u}kbl(k,s)b a (k,s) - u k d a (k,s)dl l (k,s) 



(4.42) 



At this point we realize again of the necessity of quantizing the theory using anticommutators in- 
stead of commutators. Had we use canonical commutation relations, the second term inside the 
integral in (4.42) would give the number operator dl(k, s)d a (k, s) with a minus sign in front. As a 
consequence the Hamiltonian would be unbounded from below and we would be facing again the 
instability of the theory already noticed in the context of relativistic quantum mechanics. However, 
because of the anticommutation relations (4.41), the Hamiltonian (4.42) takes the form 



H 



E 



d A k 



Lo k bl(k, s)b a (k, s) + uj k d j a (k, s)d a (k, s) - u k 5(0) 



(4.43) 



As with the scalar field, we find a divergent vacuum energy contribution due to the zero-point energy 
of the infinite number of harmonic oscillators. Unlike the Klein-Gordon field, the vacuum energy 
is negative. In section 9.2 we will see that in certain type of theories called supersymmetric, where 
the number of bosonic and fermionic degrees of freedom is the same, there is a cancellation of the 
vacuum energy. The divergent contribution can be removed by the normal order prescription 



:H:= ^ / rf3 ^ Ukbl(k,s)b a (k,s) +u k d 1; a (k,s)d a (k,s) 



(4.44) 



Finally, let us mention that using the Dirac equation it is easy to prove that there is a conserved 
four-current given by 



d,f = 0. 



(4.45) 



As we will explain further in sec. 6 this current is associated to the invariance of the Dirac La- 
grangian under the global phase shift ip — > e' t9 ip. In electrodynamics the associated conserved 
charge 



Q 



d?xf 



(4.46) 



is identified with the electric charge. 
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4.3 Gauge fields 

— * — * 

In classical electrodynamics the basic quantities are the electric and magnetic fields E, B. These 
can be expressed in terms of the scalar and vector potential (</?, A) 



B = V x A. (4.47) 

From these equations it follows that there is an ambiguity in the definition of the potentials given by 
the gauge transformations 

d - 

ip(t,x) ^ ip(t,x) + —e(t,x), A(t,x) ->■ A(t,x) + Ve(t,£). (4.48) 

Classically (</?, A) are seen as only a convenient way to solve the Maxwell equations, but without 
physical relevance. 

The equations of electrodynamics can be recast in a manifestly Lorentz invariant form using 
the four-vector gauge potential A^ = (</?, A) and the antisymmetric rank-two tensor: F^ u = d ll A v — 
d v A pi . Maxwell's equations become 

d,F^ = f, 
e^d u F av = 0, (4.49) 

where the four-current = (p,j) contains the charge density and the electric current. The field 
strength tensor and the Maxwell equations are invariant under gauge transformations (4.48), 
which in covariant form read 



A/x > A^ + d^e. 

Finally, the equations of motion of charged particles are given, in covariant form, by 



m- 



dr 



(4.50) 



(4.51) 



where e is the charge of the particle and w m (t) its four-velocity as a function of the proper time. 

The physical role of the vector potential becomes manifest only in Quantum Mechanics. Us- 
ing the prescription of minimal substitution p — >■ p — eA, the Schrodinger equation describing a 
particle with charge e moving in an electromagnetic field is 



2m V 



A 2 

ieA ) + e^> 



(4.52) 



Because of the explicit dependence on the electromagnetic potentials ip and A, this equation seems 
to change under the gauge transformations (4.48). This is physically acceptable only if the ambi- 
guity does not affect the probability density given by \^(t, x) | 2 . Therefore, a gauge transformation 
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of the electromagnetic potential should amount to a change in the (unobservable) phase of the wave 
function. This is indeed what happens: the Schrodinger equation (4.52) is invariant under the gauge 
transformations (4.48) provided the phase of the wave function is transformed at the same time 
according to 



Aharonov-Bohm effect. This interplay between gauge transformations and the phase of the 
wave function give rise to surprising phenomena. The first evidence of the role played by the 
electromagnetic potentials at the quantum level was pointed out by Yakir Aharonov and David 
Bohm [19]. Let us consider a double slit experiment as shown in Fig. 7, where we have placed a 
shielded solenoid just behind the first screen. Although the magnetic field is confined to the interior 
of the solenoid, the vector potential is nonvanishing also outside. Of course the value of A outside 
the solenoid is a pure gauge, i.e. V x A = , however because the region outside the solenoid is not 
simply connected the vector potential cannot be gauged to zero everywhere. If we denote by \T/^ 
and ^ 0) the wave functions for each of the two electron beams in the absence of the solenoid, the 
total wave function once the magnetic field is switched on can be written as 



where Ti and T 2 are two curves surrounding the solenoid from different sides, and T is any closed 
loop surrounding it. Therefore the relative phase between the two beams gets an extra term depend- 
ing on the value of the vector potential outside the solenoid as 



Because of the change in the relative phase of the electron wave functions, the presence of the 
vector potential becomes observable even if the electrons do not feel the magnetic field. If we 
perform the double-slit experiment when the magnetic field inside the solenoid is switched off we 
will observe the usual interference pattern on the second screen. However if now the magnetic field 
is switched on, because of the phase (4.54), a change in the interference pattern will appear. This is 
the Aharonov-Bohm effect. 

The first question that comes up is what happens with gauge invariance. Since we said that 
A can be changed by a gauge transformation it seems that the resulting interference patters might 
depend on the gauge used. Actually, the phase U in (4.55) is independent of the gauge although, 
unlike other gauge-invariant quantities like E and B, is nonlocal. Notice that, since Vxi = 
outside the solenoid, the value of U does not change under continuous deformations of the closed 
curve T, so long as it does not cross the solenoid. 

The Dirac monopole. It is very easy to check that the vacuum Maxwell equations remain 
invariant under the transformation 




(4.53) 



e ie Iv 1 A-dx^G) e iefr 2 A-dx^(0) 
e ie f Ti A-dx ^(0) + e ie £, A-dx^{0) 



(4.54) 




(4.55) 





(4.56) 
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Electron 
source 




Screen 

Fig. 7: Illustration of an interference experiment to show the Aharonov-Bohm effect. S represent the solenoid 
in whose interior the magnetic field is confined. 

which, in particular, for 9 = | interchanges the electric and the magnetic fields: E — >■ B, B —s- —E. 
This duality symmetry is however broken in the presence of electric sources. Nevertheless the 
Maxwell equations can be "completed" by introducing sources for the magnetic field (p m ,j m ) m 
such a way that the duality (4.56) is restored when supplemented by the transformation 

p-ip m — > e t0 (p - ip m ), j — i] m — > e l6 (j- ij m ). (4.57) 

Again for 9 = n/2 the electric and magnetic sources get interchanged. 

In 1931 Dirac [20] studied the possibility of finding solutions of the completed Maxwell 
equation with a magnetic monopoles of charge g, i.e. solutions to 

V-B = g5(x). (4.58) 

Away from the position of the monopole V ■ B = and the magnetic field can be still derived 
locally from a vector potential A according to B = V x A. However, the vector potential cannot 
be regular everywhere since otherwise Gauss law would imply that the magnetic flux threading a 
closed surface around the monopole should vanish, contradicting (4.58). 

We look now for solutions to Eq. (4.58). Working in spherical coordinates we find 

B r = jir,, B^ = B e = 0. (4.59) 

fI 

Away from the position of the monopole (x ^ 0) the magnetic field can be derived from the vector 
potential 

a 9 

= ^ tan -, A r = A e = 0. (4.60) 

\x\ 2 
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Dirac string 



Fig. 8: The Dirac monopole. 

As expected we find that this vector potential is actually singular around the half-line 9 = tt (see 
Fig. 8). This singular line starting at the position of the monopole is called the Dirac string and 
its position changes with a change of gauge but cannot be eliminated by any gauge transformation. 
Physically we can see it as an infinitely thin solenoid confining a magnetic flux entering into the 
magnetic monopole from infinity that equals the outgoing magnetic flux from the monopole. 

Since the position of the Dirac string depends on the gauge chosen it seems that the presence 
of monopoles introduces an ambiguity. This would be rather strange, since Maxwell equations 
are gauge invariant also in the presence of magnetic sources. The solution to this apparent riddle 
lies in the fact that the Dirac string does not pose any consistency problem as far as it does not 
produce any physical effect, i.e. if its presence turns out to be undetectable. From our discussion 
of the Aharonov-Bohm effect we know that the wave function of charged particles pick up a phase 
(4.55) when surrounding a region where magnetic flux is confined (for example the solenoid in the 
Aharonov-Bohm experiment). As explained above, the Dirac string associated with the monopole 
can be seen as a infinitely thin solenoid. Therefore the Dirac string will be unobservable if the phase 
picked up by the wave function of a charged particle is equal to one. A simple calculation shows 
that this happens if 

e i£ 9 = 1 => eg = 2nn with n e Z. (4.61) 

Interestingly, this discussion leads to the conclusion that the presence of a single magnetic mono- 
poles somewhere in the Universe implies for consistency the quantization of the electric charge in 
units of — , where g the magnetic charge of the monopole. 
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Quantization of the electromagnetic field. We now proceed to the quantization of the elec- 
tromagnetic field in the absence of sources p = 0, j — 0. In this case the Maxwell equations (4.49) 
can be derived from the Lagrangian density 

^Maxwell = ~\f^ = \{e 2 - B 2 ) . (4.62) 

Although in general the procedure to quantize the Maxwell Lagrangian is not very different from 
the one used for the Klein-Gordon or the Dirac field, here we need to deal with a new ingredient: 
gauge invariance. Unlike the cases studied so far, here the photon field is not unambiguously 
defined because the action and the equations of motion are insensitive to the gauge transformations 
— > + <9 M £. A first consequence of this symmetry is that the theory has less physical degrees 
of freedom than one would expect from the fact that we are dealing with a vector field. 

The way to tackle the problem of gauge invariance is to fix the freedom in choosing the 
electromagnetic potential before quantization. This can be done in several ways, for example by 
imposing the Lorentz gauge fixing condition 

d^A" = 0. (4.63) 

Notice that this condition does not fix completely the gauge freedom since Eq. (4.63) is left invariant 
by gauge transformations satisfying d^d^e = 0. One of the advantages, however, of the Lorentz 
gauge is that it is covariant and therefore does not pose any danger to the Lorentz invariance of the 
quantum theory. Besides, applying it to the Maxwell equation d^F^ = one finds 

= d^A" - d v (dpA") = 8^ A", (4.64) 

which means that since A^ satisfies the massless Klein-Gordon equation the photon, the quantum 
of the electromagnetic field, has zero mass. 

Once gauge invariance is fixed A^ is expanded in a complete basis of solutions to (4.64) and 
the canonical commutation relations are imposed 



^ f H 3 k 1 r - 

A fl (t,x) = /T^^n ^(kAMkA)e-^ t+ ^ + e,(k,XTa\k,X)e 



A=±l ' 



i\k\t—ik-x 



(4.65) 



where A = ±1 represent the helicity of the photon, and e M (fc, A) are solutions to the equations of 
motion with well defined momentum an helicity. Because of (4.63) the polarization vectors have to 
be orthogonal to 

k^e^k, A) = fc%(jfe, A)* = 0. (4.66) 
The canonical commutation relations imply that 

\a(k, \),at(k', A')] = i5{k-k')5 xy 

[a(k,X),a(k',X')} = [a\k,X),a\k',X')] = 0. (4.67) 
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Therefore a(k, \),a^(k, A) form a set of creation- annihilation operators for photons with momentum 

— * 

k and helicity A. 

Behind the simple construction presented above there are a number of subleties related with 
gauge invariance. In particular the gauge freedom seem to introduce states in the Hilbert space with 
negative probability. A careful analysis shows that when gauge invariance if properly handled these 
spurious states decouple from physical states and can be eliminated. The details can be found in 
standard textbooks [1-10]. 

Coupling gauge fields to matter. Once we know how to quantize the electromagnetic field 
we consider theories containing electrically charged particles, for example electrons. To couple 
the Dirac Lagrangian to electromagnetism we use as guiding principle what we learned about the 
Schrodinger equation for a charged particle. There we saw that the gauge ambiguity of the electro- 
magnetic potential is compensated with a U(l) phase shift in the wave function. In the case of the 
Dirac equation we know that the Lagrangian is invariant under ip — y e'^ip, with e a constant. How- 
ever this invariance is broken as soon as one identifies e with the gauge transformation parameter of 
the electromagnetic field which depends on the position. 

Looking at the Dirac Lagrangian (4.27) it is easy to see that in order to promote the global 
U(l) symmetry into a local one, ip — > e te£ ^ip, it suffices to replace the ordinary derivative <9 M by a 
covariant one satisfying 

£>„ [e^ty] = e ie£{x) D^. (4.68) 
This covariant derivative can be constructed in terms of the gauge potential as 

D, = d»- teA,. (4.69) 
The Lagrangian of a spin-| field coupled to electromagnetism is written as 

£qed = -\f^ + - m)if>, (4.70) 
invariant under the gauge transformations 

^_^ e -(^ A^^A^ + d^x). (4.71) 

Unlike the theories we have seen so far, the Lagrangian (4.70) describe an interacting theory. 
By plugging (4.69) into the Lagrangian we find that the interaction between fermions and photons 
to be 

4^ = -eA^ 7 ^. (4-72) 

As advertised above, in the Dirac theory the electric current four- vector is given by = eip'j^ip. 

The quantization of interacting field theories poses new problems that we did not meet in the 
case of the free theories. In particular in most cases it is not possible to solve the theory exactly. 
When this happens the physical observables have to be computed in perturbation theory in powers 



34 



of the coupling constant. An added problem appears when computing quantum corrections to the 
classical result, since in that case the computation of observables are plagued with infinities that 
should be taken care of. We will go back to this problem in section 8. 

Nonabelian gauge theories. Quantum electrodynamics (QED) is the simplest example of a 
gauge theory coupled to matter based in the abelian gauge symmetry of local U(l) phase rotations. 
However, it is possible also to construct gauge theories based on nonabelian groups. Actually, our 
knowledge of the strong and weak interactions is based on the use of such nonabelian generalizations 
of QED. 

Let us consider a gauge group G with generators T a , a = 1, . . . , dimG satisfying the Lie 
algebra 6 

[T a ,T b }=if abc T c . (4.73) 

A gauge field taking values on the Lie algebra of Q can be introduced = A^T a which transforms 
under a gauge transformations as 

A u — ► —UdM- 1 + UAM~\ U = e^W, (4.74) 
ig 

where g is the coupling constant. The associated field strength is defined as 

F; v = d,A a v - %Al - gf ahc A\Al. (4.75) 

Notice that this definition of the F£ v reduces to the one used in QED in the abelian case when 
jabc — o. In general, however, unlike the case of QED the field strength is not gauge invariant. In 
terms of = F^JT a it transforms as 

— ► UF^U-K (4.76) 

The coupling of matter to a nonabelian gauge field is done by introducing again a covariant 
derivative. For a field in a representation of Q 

$ — ► [/$ (4.77) 

the covariant derivative is given by 

= - tgA«T a <S>. (4.78) 

With the help of this we can write a generic Lagrangian for a nonabelian gauge field coupled to 
scalars and spinors ip as 

C = -ii^F^ + i^pi; + DrfD»<f> - i> [Mi(0) + *7 5 M 2 (0)] ^ - V {<!>). (4.79) 

In order to keep the theory renormalizable we have to restrict Mi (0) and M 2 (4>) to be at most linear 
in 4> whereas V(<j>) have to be at most of quartic order. The Lagrangian of the Standard Model is of 
the form (4.79). 

6 Some basics facts about Lie groups have been summarized in Appendix A. 



35 



4.4 Understanding gauge symmetry 

In classical mechanics the use of the Hamiltonian formalism starts with the replacement of general- 
ized velocities by momenta 

dL 

Pi = => Qi = qi(q,p)- (4.80) 

% 

Most of the times there is no problem in inverting the relations = Pi(q, q). However in some 
systems these relations might not be invertible and result in a number of constraints of the type 

f a (q,p) = 0, a = l,...,N 1 . (4.81) 

These systems are called degenerate or constrained [22,23]. 

The presence of constraints of the type (4.81) makes the formulation of the Hamiltonian for- 
malism more involved. The first problem is related to the ambiguity in defining the Hamiltonian, 
since the addition of any linear combination of the constraints do not modify its value. Secondly, 
one has to make sure that the constraints are consistent with the time evolution in the system. In the 
language of Poisson brackets this means that further constraints have to be imposed in the form 

{f a ,H}^0. (4.82) 

Following [22] we use the symbol ps to indicate a "weak" equality that holds when the constraints 
fa(q,p) = are satisfied. Notice however that since the computation of the Poisson brackets 
involves derivatives, the constraints can be used only after the bracket is computed. In principle 
the conditions (4.82) can give rise to a new set of constraints gb(q,p) = 0, b = 1, . . . , N 2 . Again 
these constraints have to be consistent with time evolution and we have to repeat the procedure. 
Eventually this finishes when a set of constraints is found that do not require any further constraint 
to be preserved by the time evolution 7 . 

Once we find all the constraints of a degenerate system we consider the so-called first class 
constraints 4> a (q,p) = 0, a = 1, . . . , M, which are those whose Poisson bracket vanishes weakly 

{(p a , <j) h } = c abc (p c pa 0. (4.83) 

The constraints that do not satisfy this condition, called second class constraints, can be eliminated 
by modifying the Poisson bracket [22]. Then the total Hamiltonian of the theory is defined by 

M 

H T = Piqi -L + Y,Ht)<l>a- (4-84) 



What has all this to do with gauge invariance? The interesting answer is that for a singular 
system the first class constraints a generate gauge transformations. Indeed, because {<fi a , <fi b } 

7 In principle it is also possible that the procedure finishes because some kind of inconsistent identity is found. In 
this case the system itself is inconsistent as it is the case with the Lagrangian L(q, q) = q. 
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{(p a , H} the transformations 



M 
a 

M 

Pi > Pi + ^2£a(t){pi,(f)a} 



(4.85) 



leave invariant the state of the system. This ambiguity in the description of the system in terms of the 
generalized coordinates and momenta can be traced back to the equations of motion in Lagrangian 
language. Writing them in the form 



d 2 L 



dqidq 



d 2 L . dL 



(4.86) 



we find that order to determine the accelerations in terms of the positions and velocities the matrix 
Q^.Qq. has to be invertible. However, the existence of constraints (4.81) precisely implies that the 
determinant of this matrix vanishes and therefore the time evolution is not uniquely determined in 
terms of the initial conditions. 

Let us apply this to Maxwell electrodynamics described by the Lagrangian 



(4.87) 



The generalized momentum conjugate to A„ is given by 



^ = 



5L 



In particular for the time component we find the constraint n° = 0. The Hamiltonian is given by 



(4.88) 



H = d 6 x [tt^o^ -£]= d 



x 



-(E 2 + B^+n°d A + A V-E 



Requiring the consistency of the constraint it = we find a second constraint 

{tt , H} w <9 7t° + V • E = 0. 



(4.89) 



(4.90) 



Together with the first constraint it = this one implies Gauss' law V-E — 0. These two constrains 
have vanishing Poisson bracket and therefore they are first class. Therefore the total Hamiltonian is 
given by 



H T = H+ d 3 x 



/' 



A^tt + A 2 (x)V • E 



(4.91) 



where we have absorbed A in the definition of the arbitrary functions X±(x) and A 2 (a;). Actually, 
we can fix part of the ambiguity taking Ai = 0. Notice that, because A has been included in the 
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multipliers, fixing Ai amounts to fixing the value of A and therefore it is equivalent to taking a 
temporal gauge. In this case the Hamiltonian is 



H T = J d 3 x i (^E 2 + B 2 ^j + e(x)V ■ E 



(4.92) 



and we are left just with Gauss' law as the only constraint. Using the canonical commutation 
relations 

{Ai{t,S),Ej{t,x')} = S^x-x') (4.93) 
we find that the remaining gauge transformations are generated by Gauss' law 

5Ai = {A u [ d 3 x' eV-E} = d^e, (4.94) 



while leaving A invariant, so for consistency with the general gauge transformations the function 
e(x) should be independent of time. Notice that the constraint V • E — can be implemented by 

— * — * — * 

demanding V • A = which reduces the three degrees of freedom of A to the two physical degrees 
of freedom of the photon. 

— * — * 

So much for the classical analysis. In the quantum theory the constraint V • E — has to be 
imposed on the physical states |phys). This is done by defining the following unitary operator on 
the Hilbert space 

m^ijtosvi.a). (4.95, 

By definition, physical states should not change when a gauge transformations is performed. This 
is implemented by requiring that the operator U(e) acts trivially on a physical state 

U{e) |phys) = |phys) (V • E) |phys> = 0. (4.96) 

In the presence of charge density p, the condition that physical states are annihilated by Gauss' law 
changes to (V • E — p)|phys) = 0. 

The role of gauge transformations in the quantum theory is very illuminating in understanding 
the real role of gauge invariance [24]. As we have learned, the existence of a gauge symmetry in a 
theory reflects a degree of redundancy in the description of physical states in terms of the degrees 
of freedom appearing in the Lagrangian. In Classical Mechanics, for example, the state of a system 
is usually determined by the value of the canonical coordinates (qi,Pi). We know, however, that 
this is not the case for constrained Hamiltonian systems where the transformations generated by the 
first class constraints change the value of q { and pi withoug changing the physical state. In the case 

— * 

of Maxwell theory for every physical configuration determined by the gauge invariant quantities E, 

— » 

B there is an infinite number of possible values of the vector potential that are related by gauge 
transformations 5A^ = d^e. 

In the quantum theory this means that the Hilbert space of physical states is defined as the 
result of identifying all states related by the operator U{e) with any gauge function e(x) into a 
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(a) (b) 

Fig. 9: Compactification of the real line (a) into the circumference S 1 (b) by adding the point at infinity. 

single physical state |phys). In other words, each physical state corresponds to a whole orbit of 
states that are transformed among themselves by gauge transformations. 

This explains the necessity of gauge fixing. In order to avoid the redundancy in the states a 
further condition can be given that selects one single state on each orbit. In the case of Maxwell 
electrodynamics the conditions A = 0, V • A = selects a value of the gauge potential among all 
possible ones giving the same value for the electric and magnetic fields. 

Since states have to be identified by gauge transformations the topology of the gauge group 
plays an important physical role. To illustrate the point let us first deal with a toy model of a U(l) 
gauge theory in 1+1 dimensions. Later we will be more general. In the Hamiltonian formalism 
gauge transformations g(x) are functions defined on R with values on the gauge group U(l) 

g:R — > 17(1). (4.97) 

We assume that g(x) is regular at infinity. In this case we can add to the real line R the point at 
infinity to compactify it into the circumference S 1 (see Fig. 9). Once this is done g(x) are functions 
defined on S 1 with values on U(l) = S 1 that can be parametrized as 

<7 : S 1 17(1), g(x)=e ia{ - x \ (4.98) 

withx G [0,2tt]. 

Because S* 1 does have a nontrivial topology, g(x) can be divided into topological sectors. 
These sectors are labelled by an integer number n G Z and are defined by 

a(2vr) = a(0) + 2nn . (4.99) 

Geometrically n gives the number of times that the spatial S 1 winds around the S* 1 defining the 
gauge group U(l). This winding number can be written in a more sophisticated way as 

& g{xY l dg(x) = 2vrn , (4.100) 
Js 1 

where the integral is along the spatial S 1 . 
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In IR 3 a similar situation happens with the gauge group 8 SU(2). If we demand g(x) G SU(2) 
to be regular at infinity \x\ — > oo we can compactify M 3 into a three-dimensional sphere S* 3 , exactly 
as we did in 1+1 dimensions. On the other hand, the function g(x) can be written as 

g(x) = a°(x)l + a{x) -a (4.101) 

and the conditions g(x) j( g(x) = 1, det g = 1 implies that (a ) 2 + a 2 = 1. Therefore SU(2) is a 
three-dimensional sphere and g(x) defines a function 

g : S 3 — > S 3 . (4.102) 

As it was the case in 1+1 dimensions here the gauge transformations g(x) are also divided into 
topological sectors labelled this time by the winding number 

n = 24^ J s3 d3xe ^ kTT [(^9) (g^^g) {g-%g)] e z. (4.103) 

In the two cases analyzed we find that due to the nontrivial topology of the gauge group 
manifold the gauge transformations are divided into different sectors labelled by an integer n. Gauge 
transformations with different values of n cannot be smoothly deformed into each other. The sector 
with n = corresponds to those gauge transformations that can be connected with the identity. 

Now we can be a bit more formal. Let us consider a gauge theory in 3+1 dimensions with 
gauge group G and let us denote by Q the set of all gauge transformations Q = {g : S 3 — > G}. At 
the same time we define Q as the set of transformations in Q that can be smoothly deformed into 
the identity. Our theory will have topological sectors if 

G/Go ^ 1. (4.104) 

In the case of the electromagnetism we have seen that Gauss' law annihilates physical states. For a 
nonabelian theory the analysis is similar and leads to the condition 



W(#o)|phys) = exp 



J d 3 x X a (x)V ■ E a 



|phys) = |phys), (4.105) 



where g (x) = e l * a{ - x ) Ta is in the connected component of the identity Q . The important point 
to realize here is that only the elements of Q Q can be written as exponentials of the infinitesimal 
generators. Since this generators annihilate the physical states this implies that W(g )|phys) = 
|phys) only when go £ Go . 

What happens then with the other topological sectors? If g e Q/Go there is still a unitary 
operator U(g) that realizes gauge transformations on the Hilbert space of the theory. However since 
g is not in the connected component of the identity, it cannot be written as the exponential of Gauss' 
law. Still gauge invariance is preserved if U(g) only changes the overall global phase of the physical 
states. For example, if g 1 is a gauge transformation with winding number n — 1 

W(<7i)|phys) =e ie |phys>. (4-106) 



8 Although we present for simplicity only the case of SU(2), similar arguments apply to any simple group. 
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It is easy to convince oneself that all transformations with winding number n — 1 have the same 
value of 9 modulo 2n. This can be shown by noticing that if g(x) has winding number n — 1 then 
c/(x) _1 has opposite winding number n = — 1. Since the winding number is additive, given two 
transformations g±, g 2 with winding number 1, g± 1 g 2 has winding number n — 0. This implies that 

|phys) = W(^-^2)|phys) = U( gi yU(g 2 )\phjs) = e^^phys) (4.107) 

and we conclude that 9 1 = 9 2 mod 2ir. Once we know this it is straightforward to conclude that a 
gauge transformation g n (x) with winding number n has the following action on physical states 

Z% n )|phys) = e in *|phys), n e Z. (4.108) 

To find a physical interpretation of this result we are going to look for similar things in other 
physical situations. One of then is borrowed from condensed matter physics and refers to the quan- 
tum states of electrons in the periodic potential produced by the ion lattice in a solid. For simplicity 
we discuss the one-dimensional case where the minima of the potential are separated by a distance 
a. When the barrier between consecutive degenerate vacua is high enough we can neglect tunnel- 
ing between different vacua and consider the ground state \na) of the potential near the minimum 
located at x = na (n e Z) as possible vacua of the theory. This vacuum state is, however, not 
invariant under lattice translations 

e iap \na) = \(n+ l)a). (4.109) 
However, it is possible to define a new vacuum state 

\k) = J2 e ~ ikna \ na )> (4.110) 

neZ 

which under e mP transforms by a global phase 

jaPiy = J2e- tkna \{n+l)a) =e lka \k). (4.111) 



e 



This ground state is labelled by the momentum k and corresponds to the Bloch wave function. 

This looks very much the same as what we found for nonabelian gauge theories. The vacuum 
state labelled by 9 plays a role similar to the Bloch wave function for the periodic potential with 
the identification of 9 with the momentum k. To make this analogy more precise let us write the 
Hamiltonian for nonabelian gauge theories 

H = \J d 3 x (7? a .7? + S -4) = \J (E a -E a + B a -B a y (4.112) 

where we have used the expression of the canonical momenta 7r* and we assume that the Gauss' law 
constraint is satisfied. Looking at this Hamiltonian we can interpret the first term within the brackets 
as the kinetic energy T = ^n a ■ n a and the second term as the potential energy V = \B a ■ B a . 

— » 

Since V > we can identify the vacua of the theory as those A for which V = 0, modulo gauge 
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transformations. This happens wherever A is a pure gauge. However, since we know that the gauge 
transformations are labelled by the winding number we can have an infinite number of vacua which 
cannot be continuously connected with one another using trivial gauge transformations. Taking a 
representative gauge transformation g n (x) in the sector with winding number n, these vacua will be 
associated with the gauge potentials 

A = -g n (x)- 1 Vg n (x), (4.113) 
W 

modulo topologically trivial gauge transformations. Therefore the theory is characterized by an 
infinite number of vacua \n) labelled by the winding number. These vacua are not gauge invariant. 
Indeed, a gauge transformation with n — 1 will change the winding number of the vacua in one unit 

U{ gi )\n) = \n+l). (4.114) 



Nevertheless a gauge invariant vacuum can be defined as 

\0) = ^e- in6 \n), withtfeM (4.115) 

satisfying 

U( 9l )\9)=e»\9). (4.116) 



We have concluded that the nontrivial topology of the gauge group have very important phys- 
ical consequences for the quantum theory. In particular it implies an ambiguity in the definition of 
the vacuum. Actually, this can also be seen in a Lagrangian analysis. In constructing the Lagrangian 
for the nonabelian version of Maxwell theory we only consider the term F" F M "°. However this is 
not the only Lorentz and gauge invariant term that contains just two derivatives. We can write the 
more general Lagrangian 



1 9 

4 c" 327T 2 M 



where F®„ is the dual of the field strength defined by 

F% = \^ a xF aX . (4.118) 

— * — * 

The extra term in (4. 1 17), proportional to E a ■ B a , is actually a total derivative and does not change 
the equations of motion or the quantum perturbation theory. Nevertheless it has several important 
physical consequences. One of them is that it violates both parity P and the combination of charge 
conjugation and parity CP. This means that since strong interactions are described by a nonabelian 
gauge theory with group SU(3) there is an extra source of CP violation which puts a strong bound 
on the value of 9. One of the consequences of a term like (4.117) in the QCD Lagrangian is a 
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nonvanishing electric dipole moment for the neutron [25]. The fact that this is not observed impose 
a very strong bound on the value of the ^-parameter 

|0| < 1(T 9 (4.119) 

From a theoretical point of view it is still to be fully understood why 9 either vanishes or has a very 
small value. 

Finally, the ^-vacuum structure of gauge theories that we found in the Hamiltonian formalism 
can be also obtained using path integral techniques form the Lagrangian (4.1 17). The second term 
in Eq. (4.1 17) gives then a contribution that depends on the winding number of the corresponding 
gauge configuration. 



5 Towards computational rules: Feynman diagrams 

As the basic tool to describe the physics of elementary particles, the final aim of Quantum Field 
Theory is the calculation of observables. Most of the information we have about the physics of 
subatomic particles comes from scattering experiments. Typically, these experiments consist of 
arranging two or more particles to collide with a certain energy and to setup an array of detectors, 
sufficiently far away from the region where the collision takes place, that register the outgoing 
products of the collision and their momenta (together with other relevant quantum numbers). 

Next we discuss how these cross sections can be computed from quantum mechanical ampli- 
tudes and how these amplitudes themselves can be evaluated in perturbative Quantum Field Theory. 
We keep our discussion rather heuristic and avoid technical details that can be found in standard 
texts [1-10]. The techniques described will be illustrated with the calculation of the cross section 
for Compton scattering at low energies. 



5.1 Cross sections and S-matrix amplitudes 

In order to fix ideas let us consider the simplest case of a collision experiment where two particles 
collide to produce again two particles in the final state. The aim of such an experiments is a direct 
measurement of the number of particles per unit time ^ (#,</?) registered by the detector flying 
within a solid angle d£l in the direction specified by the polar angles 9, (p (see Fig. 10). On general 
grounds we know that this quantity has to be proportional to the flux of incoming particles 9 , f m . 
The proportionality constant defines the differential cross section 

In natural units f m has dimensions of (length) -3 , and then the differential cross section has dimen- 
sions of (length) 2 . It depends, apart from the direction (9, ip), on the parameters of the collision 
(energy, impact parameter, etc.) as well as on the masses and spins of the incoming particles. 

Differential cross sections measure the angular distribution of the products of the collision. 
It is also physically interesting to quantify how effective the interaction between the particles is to 



9 This is defined as the number of particles that enter the interaction region per unit time and per unit area perpen- 
dicular to the direction of the beam. 
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detector 




detector 



Fig. 10: Schematic setup of a two-to-two-particles single scattering event in the center of mass reference 
frame. 

produce a nontrivial dispersion. This is measured by the total cross section, which is obtained by 
integrating the differential cross section over all directions 



To get some physical intuition of the meaning of the total cross section we can think of the classical 
scattering of a point particle off a sphere of radius R. The particle undergoes a collision only when 
the impact parameter is smaller than the radius of the sphere and a calculation of the total cross 
section yields a = tiR 2 . This is precisely the cross area that the sphere presents to incoming 
particles. 

In Quantum Mechanics in general and in Quantum Field Theory in particular the starting point 
for the calculation of cross sections is the probability amplitude for the corresponding process. In a 
scattering experiment one prepares a system with a given number of particles with definite momenta 
Pi, . . . ,p n . In the Heisenberg picture this is described by a time independent state labelled by the 
incoming momenta of the particles (to keep things simple we consider spinless particles) that we 
denote by 



On the other hand, as a result of the scattering experiment a number k of particles with momenta 
Pi, . . . , pk are detected. Thus, the system is now in the "out" Heisenberg picture state 




(5.2) 



■ ■ ■ ,P«;in). 



(5.3) 



...,p fe ';out) 



(5.4) 
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labelled by the momenta of the particles detected at late times. The probability amplitude of detect- 
ing k particles in the final state with momenta pY, . . . , p k ' in the collision of n particles with initial 
momenta pi , . . . , p„ defines the ^-matrix amplitude 

S^in ->• out) = (pi, . . . ,p fe ';out|pi, ...,p n ; in). (5.5) 

It is very important to keep in mind that both the (5.3) and (5.4) are time-independent states 
in the Hilbert space of a very complicated interacting theory. However, since both at early and late 
times the incoming and outgoing particles are well apart from each other, the "in" and "out" states 
can be thought as two states |pi, . . . ,p n ) and |pY, . . . ,p k ) of the Fock space of the corresponding 
free theory in which the coupling constants are zero. Then, the overlaps (5.5) can be written in terms 
of the matrix elements of an ^-matrix operator S acting on the free Fock space 

(pi',...,p fc ';out|pi,...,p n ;iii) = (pi,...,Pk'\S\pi,...,p n ). (5.6) 

The operator S is unitary, = S^ 1 , and its matrix elements are analytic in the external momenta. 

In any scattering experiment there is the possibility that the particles do not interact at all and 
the system is left in the same initial state. Then it is useful to write the ^-matrix operator as 

S = l + iT, (5.7) 

where 1 represents the identity operator. In this way, all nontrivial interactions are encoded in the 
matrix elements of the T-operator (pY, . . . ,p k '\iT\pi, . . . ,p n ). Since momentum has to be con- 
served, a global delta function can be factored out from these matrix elements to define the invariant 
scattering amplitude iAi 

(pi', . . . ,p k '\if |pi, ...,p n ) = (2n) 4 5^ ( J2 Pi - Y.P'f] iM ^ ■ ■ ■iPni&i ■ ■ • ( 5 - 8 ) 

V initial final / 

Total and differential cross sections can be now computed from the invariant amplitudes. Here 
we consider the most common situation in which two particles with momenta pi and p 2 collide to 
produce a number of particles in the final state with momenta p/. In this case the total cross section 
is given by 



(2u Pl )(2u P2 )\v 12 \ 



jj d 3 p' t 1 



final 
states 



(2tt)3 2cu p , 



M 



! (27r)V 4 )L+p 2 - £ A (5.9) 

\ final / 



where v 12 is the relative velocity of the two scattering particles. The corresponding differential cross 
section can be computed by dropping the integration over the directions of the final momenta. We 
will use this expression later in Section 5.3 to evaluate the cross section of Compton scattering. 

We seen how particle cross sections are determined by the invariant amplitude for the corre- 
sponding proccess, i.e. ^-matrix amplitudes. In general, in Quantum Field Theory it is not possible 
to compute exactly these amplitudes. However, in many physical situations it can be argued that in- 
teractions are weak enough to allow for a perturbative evaluation. In what follows we will describe 
how S'-matrix elements can be computed in perturbation theory using Feynman diagrams and rules. 
These are very convenient bookkeeping techniques allowing both to keep track of all contributions 
to a process at a given order in perturbation theory, and computing the different contributions. 
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5.2 Feynman rules 

The basic quantities to be computed in Quantum Field Theory are vacuum expectation values of 
products of the operators of the theory. Particularly useful are time-ordered Green functions, 



(n\T\o 1 (x 1 )...o n (x n )\\n), 

where is the the ground state of the theory and the time ordered product is defined 



T 



OiWGM = 9(x° - y )O t {x)O 3 {y) + % u - x^O^O^x). 



™0\ 



(5.10) 



(5.11) 



The generalization to products with more than two operators is straightforward: operators are always 
multiplied in time order, those evaluated at earlier times always to the right. The interest of these 
kind of correlation functions lies in the fact that they can be related to S-matrix amplitudes through 
the so-called reduction formula. To keep our discussion as simple as possible we will not derived 
it or even write it down in full detail. Its form for different theories can be found in any textbook. 
Here it suffices to say that the reduction formula simply states that any ^-matrix amplitude can be 
written in terms of the Fourier transform of a time-ordered correlation function. Morally speaking 



(pi', . . . , Prrl] OUt ...,p n ] hi) 



4 



d 4 Xl ... d 4 y n (n\TU(x^ . . . <Kz m )V(yi) • • • <Kvnj\ |n> ^' X1 ■ ■ ■ e~ M ", 



(5.12) 



where 4>(x) is the field whose elementary excitations are the particles involved in the scattering. 

The reduction formula reduces the problem of computing ^-matrix amplitudes to that of eval- 
uating time-ordered correlation functions of field operators. These quantities are easy to compute 
exactly in the free theory. For an interacting theory the situation is more complicated, however. 
Using path integrals, the vacuum expectation value of the time-ordered product of a number of 
operators can be expressed as 



(n\T o 1 (x 1 )...o n (x n ) \n) 



@<j>&<fi 1 {x 1 ) . . . O n (x n ) e lS[<t> ^ ] 



t JS[<j>,^} 



(5.13) 



For an theory with interactions, neither the path integral in the numerator or in the denominator is 
Gaussian and they cannot be calculated exactly. However, Eq. (5.13) is still very useful. The action 
S[4>, 0t] can be split into the free (quadratic) piece and the interaction part 

5[0,0 t ] = 5 o [0,0 t ] + 5 int [0,0 t ]. (5.14) 

All dependence in the coupling constants of the theory comes from the second piece. Expanding 
now exp [iSint] in power series of the coupling constant we find that each term in the series expansion 
of both the numerator and the denominator has the structure 



/ 



(5.15) 
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where ". . ." denotes certain monomial of fields. The important point is that now the integration 
measure only involves the free action, and the path integral in (5.15) is Gaussian and therefore can 
be computed exactly. The same conclusion can be reached using the operator formalism. In this case 
the correlation function (5.10) can be expressed in terms of correlation functions of operators in the 
interaction picture. The advantage of using this picture is that the fields satisfy the free equations of 
motion and therefore can be expanded in creation- annihilation operators. The correlations functions 
are then easily computed using Wick's theorem. 

Putting together all the previous ingredients we can calculate S-matrix amplitudes in a per- 
turbative series in the coupling constants of the field theory. This can be done using Feynman 
diagrams and rules, a very economical way to compute each term in the perturbative expansion of 
the S'-matrix amplitude for a given process. We will not detail the the construction of Feynman rules 
but just present them heuristically. 

For the sake of concreteness we focus on the case of QED first. Going back to Eq. (4.70) we 
expand the covariant derivative to write the action 



'QED 



d A x 



(5.16) 



The action contains two types of particles, photons and fermions, that we represent by straight and 
wavy lines respectively 



The arrow in the fermion line does not represent the direction of the momentum but the flux of 
(negative) charge. This distinguishes particles form antiparticles: if the fermion propagates from 
left to right (i.e. in the direction of the charge flux) it represents a particle, whereas when it does 
from right to left it corresponds to an antiparticle. Photons are not charged and therefore wavy lines 
do not have orientation. 

Next we turn to the interaction part of the action containing a photon field, a spinor and its 
conjugate. In a Feynman diagram this corresponds to the vertex 



Now, in order to compute an S'-matrix amplitude to a given order in the coupling constant e for a 
process with certain number of incoming and outgoing asymptotic states one only has to draw all 
possible diagrams with as many vertices as the order in perturbation theory, and the corresponding 
number and type of external legs. It is very important to keep in mind that in joining the fermion 
lines among the different building blocks of the diagram one has to respect their orientation. This 
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reflects the conservation of the electric charge. In addition one should only consider diagrams 
that are topologically non-equivalent, i.e. that they cannot be smoothly deformed into one another 
keeping the external legs fixed 10 . 

To show in a practical way how Feynman diagrams are drawn, we consider Bhabha scattering, 
i.e. the elastic dispersion of an electron and a positron: 

e + + e~ — > e + + e~. 

Our problem is to compute the ^-matrix amplitude to the leading order in the electric charge. Be- 
cause the QED vertex contains a photon line and our process does not have photons either in the 
initial or the final states we find that drawing a Feynman diagram requires at least two vertices. 
In fact, the leading contribution is of order e 2 and comes from the following two diagrams, each 
containing two vertices: 




e e e 



Incoming and outgoing particles appear respectively on the left and the right of this diagram. Notice 
how the identification of electrons and positrons is done comparing the direction of the charge 
flux with the direction of propagation. For electrons the flux of charges goes in the direction of 
propagation, whereas for positrons the two directions are opposite. These are the only two diagrams 
that can be drawn at this order in perturbation theory. It is important to include a relative minus sign 
between the two contributions. To understand the origin of this sign we have to remember that in 
the operator formalism Feynman diagrams are just a way to encode a particular Wick contraction 
of field operators in the interaction picture. The factor of —1 reflects the relative sign in Wick 
contractions represented by the two diagrams, due to the fermionic character of the Dirac field. 

We have learned how to draw Feynman diagrams in QED. Now one needs to compute the 
contribution of each one to the corresponding amplitude using the so-called Feynman rules. The 
idea is simple: given a diagram, each of its building blocks (vertices as well as external and internal 
lines) has an associated contribution that allows the calculation of the corresponding diagram. In the 
case of QED in the Feynman gauge, we have the following correspondence for vertices and internal 
propagators: 



From the point of view of the operator formalism, the requirement of considering only diagrams that are topo- 
logically nonequivalent comes from the fact that each diagram represents a certain Wick contraction in the correlation 
function of interaction-picture operators. 
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a 



j> — m + ie 




A change in the gauge would reflect in an extra piece in the photon propagator. The delta function 
implementing conservation of momenta is written using the convention that all momenta are enter- 
ing the vertex. In addition, one has to perform an integration over all momenta running in internal 
lines with the measure 




(5.17) 



and introduce a factor of —1 for each fermion loop in the diagram 11 . 

In fact, some of the integrations over internal momenta can actually be done using the delta 
function at the vertices, leaving just a global delta function implementing the total momentum con- 
servation in the diagram [cf. Eq. (5.8)]. It is even possible that all integrations can be eliminated 
in this way. This is the case when we have tree level diagrams, i.e. those without closed loops. 
In the case of diagrams with loops there will be as many remaining integrations as the number of 
independent loops in the diagram. 

The need to perform integrations over internal momenta in loop diagrams has important con- 
sequences in Quantum Field Theory. The reason is that in many cases the resulting integrals are 
ill-defined, i.e. are divergent either at small or large values of the loop momenta. In the first case 
one speaks of infrared divergences and usually they cancel once all contributions to a given pro- 
cess are added together. More profound, however, are the divergences appearing at large internal 
momenta. These ultraviolet divergences cannot be cancelled and have to be dealt through the renor- 
malization procedure. We will discuss this problem in some detail in Section 8. 

Were we computing time-ordered (amputated) correlation function of operators, this would 
be all. However, in the case of S'-matrix amplitudes this is not the whole story. In addition to the 

1 1 The contribution of each diagram comes also multiplied by a degeneracy factor that takes into account in how many 
ways a given Wick contraction can be done. In QED, however, these factors are equal to 1 for many diagrams. 
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previous rules here one needs to attach contributions also to the external legs in the diagram. These 
are the wave functions of the corresponding asymptotic states containing information about the spin 
and momenta of the incoming and outgoing particles. In the case of QED these contributions are: 



Incoming fermion: a 




Incoming antifermion: a 




Outgoing fermion: 




u a (p,s) 



Outgoing antifermion: 




Incoming photon: ^ 




Outgoing photon: 




Here we have assumed that the momenta for incoming (resp. outgoing) particles are entering (resp. 
leaving) the diagram. It is important also to keep in mind that in the computation of S'-matrix 
amplitudes all external states are on-shell. In Section 5.3 we illustrate the use of the Feynman rules 
for QED with the case of the Compton scattering. 

The application of Feynman diagrams to carry out computations in perturbation theory is ex- 
tremely convenient. It provides a very useful bookkeeping technique to account for all contributions 
to a process at a given order in the coupling constant. This does not mean that the calculation of 
Feynman diagrams is an easy task. The number of diagrams contributing to the process grows very 
fast with the order in perturbation theory and the integrals that appear in calculating loop diagrams 
also get very complicated. This means that, generically, the calculation of Feynman diagrams be- 
yond the first few orders very often requires the use of computers. 



50 



Above we have illustrated the Feynman rules with the case of QED. Similar rules can be 
computed for other interacting quantum field theories with scalar, vector or spinor fields. In the 
case of the nonabelian gauge theories introduced in Section 4.3 we have: 




permutations 



A*; a u,b 



It is not our aim here to give a full and detailed description of the Feynman rules for nonabelian 
gauge theories. It suffices to point out that, unlike the case of QED, here the gauge fields can interact 
among themselves. Indeed, the three and four gauge field vertices are a consequence of the cubic 
and quartic terms in the action 

S = ~\J d A xF^ u F flua , (5.18) 
where the nonabelian gauge field strength F° is given in Eq. (4.75). The self-interaction of the 
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nonabelian gauge fields has crucial dynamical consequences and its at the very heart of its success 
in describing the physics of elementary particles. 

5.3 An example: Compton scattering 

To illustrate the use of Feynman diagrams and Feynman rules we compute the cross section for the 
dispersion of photons by free electrons, the so-called Compton scattering: 

7 (fc, A) + e-(p, s) — ► i(k', X') + e-(p', s'). 

In brackets we have indicated the momenta for the different particles, as well as the polarizations and 
spins of the incoming and outgoing photon and electrons respectively. The first step is to identify 
all the diagrams contributing to the process at leading order. Taking into account that the vertex 
of QED contains two fermion and one photon leg, it is straightforward to realize that any diagram 
contributing to the process at hand must contain at least two vertices. Hence the leading contribution 
is of order e 2 . A first diagram we can draw is: 

P, s p', s ' 




k, X k', X 

This is, however, not the only possibility. Indeed, there is a second possible diagram: 

p, s k', X' 




k, X p', s' 

It is important to stress that these two diagrams are topologically nonequivalent, since deforming 
one into the other would require changing the label of the external legs. Therefore the leading 0(e 2 ) 
amplitude has to be computed adding the contributions from both of them. 

Using the Feynman rules of QED we find 

+ {ie) 2 u(pi, s>)0, A) /~5v 2 +m V (fr, X'yu(p, s). (5.19) 

[p — k ) — m e 

Because the leading order contributions only involve tree-level diagrams, there is no integration over 
internal momenta and therefore we are left with a purely algebraic expression for the amplitude. To 
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get an explicit expression we begin by simplifying the numerators. The following simple identity 
turns out to be very useful for this task 

= + 2(a ■ b)l. (5.20) 

Indeed, looking at the first term in Eq. (5.19) we have 

{j) + ft + m e y(k, X)u(p, s) = -0,X)(p , -m e )u(p,s) + ^(k,X)u(p,s) 

+ 2p-e(k,X)u(p,s), (5.21) 

where we have applied the identity (5.20) on the first term inside the parenthesis. The first term 
on the right-hand side of this equation vanishes identically because of Eq. (4.35). The expression 
can be further simplified if we restrict our attention to the Compton scattering at low energy when 
electrons are nonrelativistic. This means that all spatial momenta are much smaller than the electron 
mass 

\p\, \k\, \p'\, \k'\ < m e . (5.22) 

In this approximation we have that p'^ ps (m e , 0) and therefore 

p-e(jfe,A) = 0. (5.23) 

This follows from the absence of temporal photon polarization. Then we conclude that at low 
energies 

{$ + ft + m e )0, X)u(p, s) = $0, X)u(p, s) (5.24) 
and similarly for the second term in Eq. (5.19) 

+ m e )f(k', X')*u(p, s) = -$'i'{k', X')*u(p, s). (5.25) 

Next, we turn to the denominators in Eq. (5.19). As it was explained in Section 5.2, in 
computing scattering amplitudes incoming and outgoing particles should have on- shell momenta, 



p z = m z e = p' 2 and k = = k . (5.26) 
Then, the two denominator in Eq. (5.19) simplify respectively to 

(p + kf - ml = p 2 + k 2 + 2p ■ k - m\ = 2p ■ k = 2u p \k\ - 2p ■ k (5.27) 

and 

(p - k'f -m 2 e =p 2 + k' 2 + 2p ■ k' - m\ = -2p ■ k! = -2u p \k'\ + 2p ■ k'. (5.28) 
Working again in the low energy approximation (5.22) these two expressions simplify to 

(p + k) 2 - m 2 w 2m e \k\, (p - k') 2 -m\^ -2m e \k'\. (5.29) 
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Putting together all these expressions we find that at low energies 



[iey 
2m f 



ikZi\i\* $ ft 



u(p', s') W'(k'X')*^re(k, A) + e(k, \)-Z-f(k'\' s ' 



\k 



k' 



u(p,s). 



(5.30) 



Using now again the identity (5.20) a number of times as well as the transversality condition of the 
polarization vectors (4.66) we end up with a handier equation 



+ 



e 

m e 
2m, 



e(k, A) • e'(k', X')* u(p', s')^u(p, s) 



\k\ \k' 



u(p,s). (5.31) 



With a little bit of effort we can show that the second term on the right-hand side vanishes. First 
we notice that in the low energy limit \k\ ~ \k'\. If in addition we make use the conservation of 
momentum k — k' = p' — p and the identity (4.35) 



u(p', s ')0,\y(k',\y 



\k\ \k' 



u{p, s) 



-=ru{p', s')i{k, X)f(k', A')*(/ - m e )u(P, s). 
\k\ 



(5.32) 



Next we use the identity (5.20) to take the term (ft — m e ) to the right. Taking into account that in 
the low energy limit the electron four-momenta are orthogonal to the photon polarization vectors 
[see Eq. (5.23)] we conclude that 



u{p\ s'){{k, \)i\k\ \y(ft - m e )u(P, s) 

= u(p', s'){ft - m e )0, \)4'{k', X')*u(p, s) = 



(5.33) 



where the last identity follows from the equation satisfied by the conjugate positive-energy spinor, 

u(p' i s')(ft — m e ) = 0. 

After all these lengthy manipulations we have finally arrived at the expression of the invariant 
amplitude for the Compton scattering at low energies 



iM = 



ik,X)-e'(k\Xr}u(p',s')^-u(p,s). 

-I \k\ 



(5.34) 



The calculation of the cross section involves computing the modulus squared of this quantity. For 
many physical applications, however, one is interested in the dispersion of photons with a given po- 
larization by electrons that are not polarized, i.e. whose spins are randomly distributed. In addition 
in many situations either we are not interested, or there is no way to measure the final polarization 
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of the outgoing electron. This is for example the situation in cosmology, where we do not have 
any information about the polarization of the free electrons in the primordial plasma before or after 
the scattering with photons (although we have ways to measure the polarization of the scattered 
photons). 

To describe this physical situations we have to average over initial electron polarization (since 
we do not know them) and sum over all possible final electron polarization (because our detector is 
blind to this quantum number), 



\iM\ 



1 



2 \ mJk\ 



I 2 

e(k, A) • e'(k', A')* E E \u{p' , s')jtu{p, s) 



(5.35) 



The factor of \ comes from averaging over the two possible polarizations of the incoming electrons. 
The sums in this expression can be calculated without much difficulty. Expanding the absolute value 
explicitly 

E E |u(MMP>)f= E E [<P,s)^W',s'y][u(p',s'^u(P,s)], (5.36) 



s =±h s '=±h 



using that 7^ = 7°7^7° and after some manipulation one finds that 



E E \u{p' , s')jtu{p, s) 



s=±i s'=±i 



E U a (p,s)Uf3(p,s) 



Pa 



pa 



Tr [(p l + m e )lk{p" + m e )ltj, 



(5.37) 



where the final expression has been computed using the completeness relations in Eq. (4.38). The 
final evaluation of the trace can be done using the standard Dirac matrices identities. Here we 
compute it applying again the relation (5.20) to commute ft and p. Using that k 2 = and that we 
are working in the low energy limit we have 12 



Tr 



{ft) + m e )W + m e )f = 2(p • k)(p' ■ k)Tr 1 w 8m 2 e \k\ 2 . 



(5.38) 



This gives the following value for the invariant amplitude 



\iM\ 2 = 4e 4 



e(k,X) ■e'(fc',A')* 



(5.39) 



Plugging \iM | 2 into the formula for the differential cross section we get 

.2 \ 2 



da 



dtt QA7r 2 m 2 e 



\iM\ 



Aura, 



e(k,X) ■e'ik',^)* 



(5.40) 



12 We use also the fact that the trace of the product of an odd number of Dirac matrices is always zero. 
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The prefactor of the last equation is precisely the square of the so-called classical electron radius 
r cl . In fact, the previous differential cross section can be rewritten as 



da 3 

dn = 8^ aT 

where a T is the total Thomson cross section 



e(k, X)-e'{k',X'y 



8tt 



(Jt 



cl- 



(5.41) 



(5.42) 



The result (5.41) is relevant in many areas of Physics, but its importance is paramount in the 
study of the cosmological microwave background (CMB). Just before recombination the universe is 
filled by a plasma of electrons interacting with photons via Compton scattering, with temperatures 
of the order of 1 keV. Electrons are then nonrelativistic (m e ~ 0.5 MeV) and the approximations 
leading to Eq. (5.41) are fully valid. Because we do not know the polarization state of the photons 
before being scattered by electrons we have to consider the cross section averaged over incoming 
photon polarizations. From Eq. (5.41) we see that this is proportional to 



\ \«k, X) ■ e'(k>, A')* 2 = \ Y1 e A A) 



A=l,2 



A=l,2 



(5.43) 



The sum inside the brackets can be computed using the normalization of the polarization vectors, 
\e(k, A)| 2 = 1, and the transversality condition k ■ e(k, A) = 



lYl kk,X).e'(k',X') 



A=l,2 



1 r 

2 



k\ 2 

l-\£-e'(k\X')\ 2 



(5.44) 



where £ = 4- is the unit vector in the direction of the incoming photon. 

\k\ 

From the last equation we conclude that Thomson scattering suppresses all polarizations par- 
allel to the direction of the incoming photon I, whereas the differential cross section reaches the 
maximum in the plane normal to I. If photons would collide with the electrons in the plasma with 
the same intensity from all directions, the result would be an unpolarized CMB radiation. The 
fact that polarization is actually measured in the CMB carries crucial information about the physics 
of the plasma before recombination and, as a consequence, about the very early universe (see for 
example [21] for a throughout discussion). 



6 Symmetries 

6.1 Noether's theorem 

In Classical Mechanics and Classical Field Theory there is a basic result that relates symmetries and 
conserved charges. This is called Noether's theorem and states that for each continuous symmetry 
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of the system there is conserved current. In its simplest version in Classical Mechanics it can be 
easily proved. Let us consider a Lagrangian L{q i ,q i ) which is invariant under a transformation 
Qi(f) ~~ ^ Qi(ti e ) labelled by a parameter e. This means that L(q',q') = L(q,q) without using the 
equations of motion 13 . If e 1 we can consider an infinitesimal variation of the coordinates 8 e qi(t) 
and the invariance of the Lagrangian implies 



= 5 t L(q h qi 



<9L dL _ . 

a~ S eQi + w^O e qi 



dL 

dqi 



d dL 

dt dqi 



d (dL s 
* e * + dt {d* 5 * 



(6.1) 



When 5 e qi is applied on a solution to the equations of motion the term inside the square brackets 
vanishes and we conclude that there is a conserved quantity 



dL 

Q = with Q = —S e qi 
dqi 



(6.2) 



Notice that in this derivation it is crucial that the symmetry depends on a continuous parameter since 
otherwise the infinitesimal variation of the Lagrangian in Eq. (6.1) does not make sense. 

In Classical Field Theory a similar result holds. Let us consider for simplicity a theory of 
a single field <j){x). We say that the variations 5 t (p depending on a continuous parameter e are a 
symmetry of the theory if, without using the equations of motion, the Lagrangian density changes 

by 



5X = d u K». 



(6.3) 



If this happens then the action remains invariant and so do the equations of motion. Working out 
now the variation of C under 5 e <fi we find 



dC 



d„ 



dC 



d(d,<P) 



5 e cf>. (6.4) 



If <j)(x) is a solution to the equations of motion the last terms disappears, and we find that there is a 
conserved current 



d^ = with J» 



dC 



-5 e (j) - K". 



did^) 

Actually a conserved current implies the existence of a charge 



(6.5) 



which is conserved 



^ = J d 3 xd J°(t,x) = - J cPxdiJtfag) = 0, 



(6.6) 



(6.7) 



13 The following result can be also derived a more general situations where the Lagrangian changes by a total time 
derivative. 
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provided the fields vanish at infinity fast enough. Moreover, the conserved charge Q is a Lorentz 
scalar. After canonical quantization the charge Q defined by Eq. (6.6) is promoted to an operator 
that generates the symmetry on the fields 

6(/> = i[(l>,Q]. (6.8) 

As an example we can consider a scalar field <j>(x) which under a coordinate transformation 
x — > x' changes as (j)'(x') = 4>(x). In particular performing a space-time translation x^' = x^ + a M 
we have 

<j>'{x) - <j>(x) = -a"d^<f) + 0(a 2 ) 5<P = -a^d^. (6.9) 

Since the Lagrangian density is also a scalar quantity, it transforms under translations as 

5C = -a^d^C. (6.10) 
Therefore the corresponding conserved charge is 

dC 

where we introduced the energy-momentum tensor 

dC 

We find that associated with the invariance of the theory with respect to space-time translations 
there are four conserved currents defined by T^ v with v — 0, . . . , 3, each one associated with the 
translation along a space-time direction. These four currents form a rank-two tensor under Lorentz 
transformations satisfying 

d li T' u ' = 0. (6.13) 
The associated conserved charges are given by 

P v = [ d 3 xT 0u (6.14) 



and correspond to the total energy-momentum content of the field configuration. Therefore the 
energy density of the field is given by T 00 while T° l is the momentum density. In the quantum 
theory the P M are the generators of space-time translations. 

Another example of a symmetry related with a physically relevant conserved charge is the 
global phase invariance of the Dirac Lagrangian (4.27), ip — > e t9 ^jj. For small 9 this corresponds to 
variations Sgip = i9t/j, = —iO^p which by Noether's theorem result in the conserved charge 

f = ^f^, d,f = 0. (6.15) 
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Thus implying the existence of a conserved charge 



Q = J d 3 xtp~f°ip = J (fx^i). (6.16) 

In physics there are several instances of global U(l) symmetries that act as phase shifts on spinors. 
This is the case, for example, of the baryon and lepton number conservation in the Standard Model. 
A more familiar case is the U(l) local symmetry associated with electromagnetism. Notice that 
although in this case we are dealing with a local symmetry, 9 — > ea(x), the invariance of the 
Lagrangian holds in particular for global transformations and therefore there is a conserved current 
ji* = eip^ip. In Eq. (4.72) we saw that the spinor is coupled to the photon field precisely through 
this current. Its time component is the electric charge density p, while the spatial components are 
the current density vector j. 

This analysis can be carried over also to nonabelian unitary global symmetries acting as 

fr^HJitfj, U*U = 1 (6.17) 

and leaving invariant the Dirac Lagrangian when we have several fermions. If we write the matrix 
U in terms of the hermitian group generators T a as 

U = exp(ia a T a ), (T a )t = T a , (6.18) 

we find the conserved current 

j^ = ^7%, 9^ = 0. (6.19) 

This is the case, for example of the approximate flavor symmetries in hadron physics. The simplest 
example is the isospin symmetry that mixes the quarks u and d 

U d ^j — ► M ( ^ , Me SU(2). (6.20) 

Since the proton is a bound state of two quarks u and one quark d while the neutron is made out 
of one quark u and two quarks d, this isospin symmetry reduces at low energies to the well known 
isospin transformations of nuclear physics that mixes protons and neutrons. 

6.2 Symmetries in the quantum theory 

We have seen that in canonical quantization the conserved charges Q a associated to symmetries 
by Noether's theorem are operators implementing the symmetry at the quantum level. Since the 
charges are conserved they must commute with the Hamiltonian 

[Q a ,H] = 0. (6.21) 

There are several possibilities in the quantum mechanical realization of a symmetry: 
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Wigner-Weyl realization. In this case the ground state of the theory |0) is invariant under the 
symmetry. Since the symmetry is generated by Q a this means that 

U(a)\0) = e iaaQa \0) = |0) =^ Q a \0) = 0. (6.22) 

At the same time the fields of the theory have to transform according to some irreducible represen- 
tation of the group generated by the Q a . From Eq. (6.8) it is easy to prove that 

U(a)faU(a)~ 1 = Uij(a)<l>j, (6.23) 

where £4, (a) is an element of the representation in which the field fa transforms. If we consider 
now the quantum state associated with the operator fa 

\i) = fa\0) (6-24) 

we find that because of the invariance of the vacuum (6.22) the states \i) transform in the same 
representation as fa 

U(a)\i) = U(a)faU(a)- l U(a)\Q) = U i:j (a)fa\0) = U t ,(a)\j). (6.25) 

Therefore the spectrum of the theory is classified in multiplets of the symmetry group. In addition, 
since [H, 14(a)] = all states in the same multiplet have the same energy. If we consider one- 
particle states, then going to the rest frame we conclude that all states in the same multiplet have 
exactly the same mass. 

Nambu-Goldstone realization. In our previous discussion the result that the spectrum of 
the theory is classified according to multiplets of the symmetry group depended crucially on the 
invariance of the ground state. However this condition is not mandatory and one can relax it to 
consider theories where the vacuum state is not left invariant by the symmetry 

e iaaQa \0) ^\0) Q a \0)^0. (6.26) 

In this case it is also said that the symmetry is spontaneously broken by the vacuum. 

To illustrate the consequences of (6.26) we consider the example of a number scalar fields Lp l 
(i — 1, . . . , N) whose dynamics is governed by the Lagrangian 

C = ±d li <p i d»<p i -V(<p), (6.27) 

where we assume that V(fa) is bounded from below. This theory is globally invariant under the 
transformations 

V = e a (T a )y, (6.28) 
with T a , a = 1, . . . , ^N(N - 1) the generators of the group SO(N). 



To analyze the structure of vacua of the theory we construct the Hamiltonian 



H = I d 3 x 



ivrV + ivy?* • + V(ip) 



(6.29) 
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and look for the minimum of 



d 3 x 



• Vy? 1 + V{<p) 



(6.30) 



Since we are interested in finding constant field configurations, Vy? = to preserve translational 
invariance, the vacua of the potential V((f) coincides with the vacua of V(<p). Therefore the minima 
of the potential correspond to the vacuum expectation values 14 



W» = o, 



dV 



dip* 



= 0. 



(6.31) 



<pi = {<pi) 



We divide the generators T a of SO(iV) into two groups: Those denoted by H a (a — 1, ... ,h) 
that satisfy 



(/n;v> = o. 



(6.32) 



This means that the vacuum configuration is left invariant by the transformation generated 
by H a . For this reason we call them unbroken generators. Notice that the commutator of two 
unbroken generators also annihilates the vacuum expectation value, [H a , H^]ij{(p j ) = 0. Therefore 
the generators {H a } form a subalgebra of the algebra of the generators of SO(iV). The subgroup of 
the symmetry group generated by them is realized a la Wigner-Weyl. 

The remaining generators K A , with A = 1, . . . , |iV(iV — l)—h, by definition do not preserve 
the vacuum expectation value of the field 



(6.33) 



These will be called the broken generators. Next we prove a very important result concerning 
the broken generators known as the Goldstone theorem: for each generator broken by the vacuum 
expectation value there is a massless excitation. 

The mass matrix of the excitations around the vacuum (</?*) is determined by the quadratic 
part of the potential. Since we assumed that V({ip)) = and we are expanding around a minimum, 
the first term in the expansion of the potential V((p) around the vacuum expectation values is given 
by 



V{tp) 



and the mass matrix is: 



d 2 V 



¥>=<¥>} 



tf-tf)W-(<p?)) + 0[{<p-(<p)Y\ 



(6.34) 



M ■ 



d 2 V 



(6.35) 



<p={tp) 



In order to avoid a cumbersome notation we do not show explicitly the dependence of the mass 
matrix on the vacuum expectation values (</? 1 ). 



14 



For simplicity we consider that the minima of V(<p) occur at zero potential. 
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To extract some information about the possible zero modes of the mass matrix, we write down 
the conditions that follow from the invariance of the potential under 5(p % = e a (T a )'j(p j . At first order 
in e a 

5V(<p)=e a ^(T a )y = 0. (6.36) 
Differentiating this expression with respect to (p k we arrive at 

+ W (T ° )l = °- <6 ' 37) 

Now we evaluate this expression in the vacuum tp 1 = (</?*). Then the derivative in the second term 
cancels while the second derivative in the first one gives the mass matrix. Hence we find 

A4(n;v> = o. (6.38) 

Now we can write this expression for both broken and unbroken generators. For the unbroken ones, 
since (iP)* (y^) = 0, we find a trivial identity = 0. On the other hand for the broken generators 
we have 

Ml (K A )){^) = 0. (6.39) 

Since {K A ) l -{^) ^ this equation implies that the mass matrix has as many zero modes as broken 
generators. Therefore we have proven Goldstone's theorem: associated with each broken symmetry 
there is a massless mode in the theory. Here we have presented a classical proof of the theorem. In 
the quantum theory the proof follows the same lines as the one presented here but one has to consider 
the effective action containing the effects of the quantum corrections to the classical Lagrangian. 

As an example to illustrate this theorem, we consider a SO(3) invariant scalar field theory with 
a "mexican hat" potential 

V{0) = \{$ 2 -<?)\ (6.40) 

The vacua of the theory correspond to the configurations satisfying (<p) 2 = a 2 . In field space 
this equation describes a two-dimensional sphere and each solution is just a point in that sphere. 
Geometrically it is easy to visualize that a given vacuum field configuration, i.e. a point in the 
sphere, is preserved by SO(2) rotations around the axis of the sphere that passes through that point. 
Hence the vacuum expectation value of the scalar field breaks the symmetry according to 

(0) : SO(3) — > SO(2). (6.41) 

Since SO(3) has three generators and SO(2) only one we see that two generators are broken and 
therefore there are two massless Goldstone bosons. Physically this massless modes can be thought 
of as corresponding to excitations along the surface of the sphere (<p) 2 = a 2 . 

Once a minimum of the potential has been chosen we can proceed to quantize the excitations 
around it. Since the vacuum only leaves invariant a SO(2) subgroup of the original SO(3) symmetry 
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group it seems that the fact that we are expanding around a particular vacuum expectation value of 
the scalar field has resulted in a lost of symmetry. This is however not the case. The full quantum 
theory is symmetric under the whole symmetry group SO(3). This is reflected in the fact that the 
physical properties of the theory do not depend on the particular point of the sphere (<p) 2 = a 2 that 
we have chosen. Different vacua are related by the full SO(3) symmetry and therefore should give 
the same physics. 

It is very important to realize that given a theory with a vacuum determined by (if) all other 
possible vacua of the theory are unaccessible in the infinite volume limit. This means that two 
vacuum states |0i), |0 2 ) corresponding to different vacuum expectation values of the scalar field are 
orthogonal (0i|0 2 ) = and cannot be connected by any local observable <&(x), (Oi|$(:r) |0 2 ) = 0. 
Heuristically this can be understood by noticing that in the infinite volume limit switching from one 
vacuum into another one requires changing the vacuum expectation value of the field everywhere 
in space at the same time, something that cannot be done by any local operator. Notice that this is 
radically different to our expectations based on the Quantum Mechanics of a system with a finite 
number of degrees of freedom. 

In High Energy Physics the typical example of a Goldstone boson is the pion, associated with 
the spontaneous breaking of the global chiral isospin SU(2) L x SU(2) R symmetry. This symmetry 
acts independently in the left- and right-handed spinors as 



Presumably since the quarks are confined at low energies this symmetry is spontaneously broken 
down to the diagonal SU(2) acting in the same way on the left- and right-handed components of 
the spinors. Associated with this symmetry breaking there is a Goldstone mode which is identified 
as the pion. Notice, nevertheless, that the SU(2) L xSU(2) R would be an exact global symmetry of 
the QCD Lagrangian only in the limit when the masses of the quarks are zero m u , m d — > 0. Since 
these quarks have nonzero masses the chiral symmetry is only approximate and as a consequence 
the corresponding Goldstone boson is not massless. That is why pions have masses, although they 
are the lightest particle among the hadrons. 

Symmetry breaking appears also in many places in condensed matter. For example, when a 
solid crystallizes from a liquid the translational invariance that is present in the liquid phase is broken 
to a discrete group of translations that represent the crystal lattice. This symmetry breaking has 
Goldstone bosons associated which are identified with phonons which are the quantum excitation 
modes of the vibrational degrees of freedom of the lattice. 

The Higgs mechanism. Gauge symmetry seems to prevent a vector field from having a mass. 
This is obvious once we realize that a term in the Lagrangian like m 2 A fl A fi is incompatible with 
gauge invariance. 

However certain physical situations seem to require massive vector fields. This happened for 
example during the 1960s in the study of weak interactions. The Glashow model gave a common 
description of both electromagnetic and weak interactions based on a gauge theory with group 
SU(2)xU(l) but, in order to reproduce Fermi's four-fermion theory of the /3-decay it was necessary 
that two of the vector fields involved would be massive. Also in condensed matter physics massive 





(6.42) 
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vector fields are required to describe certain systems, most notably in superconductivity. 

The way out to this situation is found in the concept of spontaneous symmetry breaking dis- 
cussed previously. The consistency of the quantum theory requires gauge invariance, but this in- 
variance can be realized a la Nambu-Goldstone. When this is the case the full gauge symmetry is 
not explicitly present in the effective action constructed around the particular vacuum chosen by the 
theory. This makes possible the existence of mass terms for gauge fields without jeopardizing the 
consistency of the full theory, which is still invariant under the whole gauge group. 

To illustrate the Higgs mechanism we study the simplest example, the Abelian Higgs model: 
a U(l) gauge field coupled to a self-interacting charged complex scalar field $ with Lagrangian 

C = —F^F^ + D^D»$ - j ($$ - /J 2 ) 2 , (6.43) 

where the covariant derivative is given by Eq. (4.69). This theory is invariant under the gauge 
transformations 

$_> e **(*)$, A^A^ + d^x). (6.44) 

The minimum of the potential is defined by the equation |$| = /x. We have a continuum of different 
vacua labelled by the phase of the scalar field. None of these vacua, however, is invariant under the 
gauge symmetry 

($) = ite?* ->• fie Mo+ia{x) (6.45) 

and therefore the symmetry is spontaneously broken Let us study now the theory around one of 
these vacua, for example ($} = fx, by writing the field $ in terms of the excitations around this 
particular vacuum 

(6.46) 

Independently of whether we are expanding around a particular vacuum for the scalar field we 
should keep in mind that the whole Lagrangian is still gauge invariant under (6.44). This means that 
performing a gauge transformation with parameter a(x) = —i3(x) we can get rid of the phase in 
Eq. (6.46). Substituting then — pL+ -^<r(x) in the Lagrangian we find 

C = -- A F, u F^ + e 2 ii 2 A^+ l -d,ad»a-^ 2 a 2 

- A/i<T 3 — ^<T 4 + e 2 fiA fM A fl a + e 2 A fl A fJl a 2 . (6.47) 

What are the excitation of the theory around the vacuum ($} = /i? First we find a massive real 
scalar field a{x). The important point however is that the vector field A^ now has a mass given by 

ml = 2e 2 /i 2 . (6.48) 
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The remarkable thing about this way of giving a mass to the photon is that at no point we have given 
up gauge invariance. The symmetry is only hidden. Therefore in quantizing the theory we can still 
enjoy all the advantages of having a gauge theory but at the same time we have managed to generate 
a mass for the gauge field. 

It is surprising, however, that in the Lagrangian (6.47) we did not found any massless mode. 
Since the vacuum chosen by the scalar field breaks the U(l) generator of U(l) we would have 
expected one masless particle from Goldstone's theorem. To understand the fate of the missing 
Goldstone boson we have to revisit the calculation leading to Eq. (6.47). Were we dealing with a 
global U(l) theory, the Goldstone boson would correspond to excitation of the scalar field along the 
valley of the potential and the phase would be the massless Goldstone boson. However we have 
to keep in mind that in computing the Lagrangian we managed to get rid of $(x) by shifting it into 
using a gauge transformation. Actually by identifying the gauge parameter with the Goldstone 
excitation we have completely fixed the gauge and the Lagrangian (6.47) does not have any gauge 
symmetry left. 

— # — # 

A massive vector field has three polarizations: two transverse ones k • e(k, ±1) = plus a 
longitudinal one ex,(fc) ~ k. In gauging away the massless Goldstone boson we have trans- 
formed it into the longitudinal polarization of the massive vector field. In the literature this is usually 
expressed saying that the Goldstone mode is "eaten up" by the longitudinal component of the gauge 
field. It is important to realize that in spite of the fact that the Lagrangian (6.47) looks pretty dif- 
ferent from the one we started with we have not lost any degrees of freedom. We started with the 
two polarizations of the photon plus the two degrees of freedom associated with the real and imag- 
inary components of the complex scalar field. After symmetry breaking we end up with the three 
polarizations of the massive vector field and the degree of freedom of the real scalar field a(x). 

We can also understand the Higgs mechanism in the light of our discussion of gauge symme- 
try in section 4.4. In the Higgs mechanism the invariance of the theory under infinitesimal gauge 
transformations is not explicitly broken, and this implies that Gauss' law is satisfied quantum me- 
chanically, V • E a \phys) = 0. The theory remains invariant under gauge transformations in the 
connected component of the identity Go, the ones generated by Gauss' law. This does not pose any 
restriction on the possible breaking of the invariance of the theory with respect to transformations 
that cannot be continuously deformed to the identity. Hence in the Higgs mechanism the invariance 
under gauge transformation that are not in the connected component of the identity, G /Go, can be 
broken. Let us try to put it in more precise terms. As we learned in section 4.4, in the Hamiltonian 
formulation of the theory finite energy gauge field configurations tend to a pure gauge at spatial 
infinity 

A^x^-gix^Vgix), \x\ -> oo (6.49) 
%g 

The set transformations go(x) G Go that tend to the identity at infinity are the ones generated by 
Gauss' law. However, one can also consider in general gauge transformations g(x) which, as \x\ — >■ 
oo, approach any other element g e G. The quotient Goo = G /Go gives a copy of the gauge group 
at infinity. There is no reason, however, why this group should not be broken, and in general it is if 
the gauge symmetry is spontaneously broken. Notice that this is not a threat to the consistency of 
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the theory. Properties like the decoupling of unphysical states are guaranteed by the fact that Gauss' 
law is satisfied quantum mechanically and are not affected by the breaking of Qoo. 

The Abelian Higgs model discussed here can be regarded as a toy model of the Higgs mech- 
anism responsible for giving mass to the W ± and Z G gauge bosons in the Standard Model. In 
condensed matter physics the symmetry breaking described by the nonrelativistic version of the 
Abelian Higgs model can be used to characterize the onset of a superconducting phase in the BCS 
theory, where the complex scalar field <3> is associated with the Cooper pairs. In this case the param- 
eter fj 2 depends on the temperature. Above the critical temperature T c , yU 2 (T) > and there is only 
a symmetric vacuum ($) = 0. When, on the other hand, T < T c then fi 2 (T) < and symmetry 
breaking takes place. The onset of a nonzero mass of the photon (6.48) below the critical tem- 
perature explains the Meissner effect: the magnetic fields cannot penetrate inside superconductors 
beyond a distance of the order — . 

7 Anomalies 

So far we did not worry too much about how classical symmetries of a theory are carried over to the 
quantum theory. We have implicitly assumed that classical symmetries are preserved in the process 
of quantization, so they are also realized in the quantum theory. 

This, however, does not have to be necessarily the case. Quantizing an interacting field theory 
is a very involved process that requires regularization and renormalization and sometimes, it does 
not matter how hard we try, there is no way for a classical symmetry to survive quantization. When 
this happens one says that the theory has an anomaly (for a review see [27]). It is important to avoid 
here the misconception that anomalies appear due to a bad choice of the way a theory is regularized 
in the process of quantization. When we talk about anomalies we mean a classical symmetry that 
cannot be realized in the quantum theory, no matter how smart we are in choosing the regularization 
procedure. 

In the following we analyze some examples of anomalies associated with global and local 
symmetries of the classical theory. In Section 8 we will encounter yet another example of an 
anomaly, this time associated with the breaking of classical scale invariance in the quantum the- 
ory. 

7.1 Axial anomaly 

Probably the best known examples of anomalies appear when we consider axial symmetries. If we 
consider a theory of two Weyl spinors u± 



the Lagrangian is invariant under two types of global U(l) transformations. In the first one both 
helicities transform with the same phase, this is a vector transformation: 



C = iip^tp = iu' + a1d^,u + 





(7.1) 




(7.2) 



66 



whereas in the second one, the axial U(l), the signs of the phases are different for the two chiralities 

U(1) A : u± — ► e ±ia u±. (7.3) 
Using Noether's theorem, there are two conserved currents, a vector current 

J£ = = u ] + a^u + + ulatu- => 0^ = (7.4) 

and an axial vector current 

J» = ^f lf $ = u \.(T%u+ - ulatu- 8^ = 0. (7.5) 

The theory described by the Lagrangian (7.1) can be coupled to the electromagnetic field. 
The resulting classical theory is still invariant under the vector and axial U(l) symmetries (7.2) and 
(7.3). Surprisingly, upon quantization it turns out that the conservation of the axial current (7.5) is 
spoiled by quantum effects 

d^-HE-B. (7.6) 

To understand more clearly how this result comes about we study first a simple model in 
two dimensions that captures the relevant physics involved in the four-dimensional case [28]. We 
work in Minkowski space in two dimensions with coordinates (x ,^ 1 ) = (t,x) and where the 
spatial direction is compactified to a circle S 1 . In this setup we consider a fermion coupled to the 
electromagnetic field. Notice that since we are living in two dimensions the field strength F^ v only 
has one independent component that corresponds to the electric field along the spatial direction, 
F 01 = E (in two dimensions there are no magnetic fields!). 

To write the Lagrangian for the spinor field we need to find a representation of the algebra of 
7-matrices 

Wf} = *r with 77= J (7.7) 
In two dimensions the dimension of the representation of the 7-matrices is 2^1 = 2. Here take 

/«>-(:;), tw-(_;;)- 

This is a chiral representation since the matrix 75 is diagonal 15 

75 - -7V =(\ _l) (7-9) 



Writing the two-component spinor ip as 



(7.10) 



15 In any even number of dimensions 75 is denned to satisfy the conditions 7! = 1 and {75, 7^} = 0. 
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Fig. 1 1 : Spectrum of the massless two-dimensional Dirac field. 



and defining as usual the projectors P± = |(1 ± 75) we find that the components u± of ip are 
respectively a right- and left-handed Weyl spinor in two dimensions. 

Once we have a representation of the 7-matrices we can write the Dirac equation. Expressing 
it in terms of the components u± of the Dirac spinor we find 

(d o -d 1 )u+ = 0, (d o + d 1 )u- = 0. (7.11) 

The general solution to these equations can be immediately written as 

u + = u + (x° + u_ = u_(x° -x 1 ). (7.12) 

Hence u± are two wave packets moving along the spatial dimension respectively to the left (u+) 
and to the right Notice that according to our convention the left-moving u + is a right-handed 
spinor (positive helicity) whereas the right-moving w_ is a left-handed spinor (negative helicity). 

If we want to interpret (7.1 1) as the wave equation for two-dimensional Weyl spinors we have 
the following wave functions for free particles with well defined momentum p M = (E,p). 

u { f\x° ± x l ) = -Le-tf^i* 1 ) W ith p = T E. (7.13) 
y/L 

As it is always the case with the Dirac equation we have both positive and negative energy solutions. 
For u + , since E = —p, we see that the solutions with positive energy are those with negative 
momentum p < 0, whereas the negative energy solutions are plane waves with p > 0. For the 
left-handed spinor w_ the situation is reversed. Besides, since the spatial direction is compact with 
length L the momentum p is quantized according to 

p=^, neZ. (7.14) 
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Fig. 12: Vacuum of the theory. 



The spectrum of the theory is represented in Fig. 1 1 . 

Once we have the spectrum of the theory the next step is to obtain the vacuum. As with the 
Dirac equation in four dimensions we fill all the states with E < (Fig. 12). Exciting of a particle 
in the Dirac see produces a positive energy fermion plus a hole that is interpreted as an antiparticle. 
This gives us the clue on how to quantize the theory. In the expansion of the operator u± in terms of 
the modes (7.13) we associate positive energy states with annihilation operators whereas the states 
with negative energy are associated with creation operators for the corresponding antiparticle 



The operator a±(E) acting on the vacuum |0, ±) annihilates a particle with positive energy E and 
momentum ^fE. In the same way b±(E) creates out of the vacuum an antiparticle with positive 
energy E and spatial momentum =pE. In the Dirac sea picture the operator b±{E)^ is originally an 
annihilation operator for a state of the sea with negative energy — E. As in the four-dimensional 
case the problem of the negative energy states is solved by interpreting annihilation operators for 
negative energy states as creation operators for the corresponding antiparticle with positive energy 
(and vice versa). The operators appearing in the expansion of u± in Eq. (7.15) satisfy the usual 
algebra 



where we have introduced the label A, A' = ±. Also, a\{E), a\{E) anticommute with by(E'), 




(7.15) 



E>0 



{a x (E),a{,(E?)} = {b x (E),b[,(E')} = 5 E , E ,5 XX ,, 



(7.16) 



b[,(E<). 



The Lagrangian of the theory 



C = iu\_(d + d\)u + + iu_(d — d\)u 



(7.17) 
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is invariant under both U(l)y, Eq. (7.2), and U(l)^, Eq. (7.3). The associated Noether currents are 
in this case 

\ —u' + u + + u_u^ j \ —u' + u + — U_U- J 

The associated conserved charges are given, for the vector current by 

Q v = J dx 1 (u\u + + «Lu_) (7.19) 

and for the axial current 

Q A = dx 1 (u\u+ - wl«_) . (7.20) 

Using the orthonormality relations for the modes v±\x) 

c 

dx 1 vf\x) vf'\x) = 8 e ,e> (7.21) 
we find for the conserved charges: 

Qv = Y, [al(E)a + (E)-bi(E)b + (E) + al(E)a-(E)-bl(E)b-(E) , 

E>0 

Qa = Y [al(E)a+(E)-bl(E)b + (E)-al(E)a-(E) + bl(E)b4E) . 



E>0 



(7.22) 



We see that Q v counts the net number (particles minus antiparticles) of positive helicity states plus 
the net number of states with negative helicity. The axial charge, on the other hand, counts the 
net number of positive helicity states minus the number of negative helicity ones. In the case of 
the vector current we have subtracted a formally divergent vacuum contribution to the charge (the 
"charge of the Dirac sea"). 

In the free theory there is of course no problem with the conservation of either Qv or Qa, 
since the occupation numbers do not change. What we want to study is the effect of coupling the 
theory to electric field £. We work in the gauge A = 0. Instead of solving the problem exactly we 
are going to simulate the electric field by adiabatically varying in a long time r the vector potential 
A 1 from zero value to — £t . From our discussion in section 4.3 we know that the effect of the 
electromagnetic coupling in the theory is a shift in the momentum according to 

p^p- eA u (7.23) 

where e is the charge of the fermions. Since we assumed that the vector potential varies adiabati- 
cally, we can assume it to be approximately constant at each time. 

Then, we have to understand what is the effect of (7.23) on the vacuum depicted in Fig. 
(12). What we find is that the two branches move as shown in Fig. (13) resulting in some of 
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Fig. 13: Effect of the electric field. 



the negative energy states of the v + branch acquiring positive energy while the same number of 
the empty positive energy states of the other branch V- will become empty negative energy states. 
Physically this means that the external electric field £ creates a number of particle-antiparticle pairs 
out of the vacuum. Denoting by N ~ e£ the number of such pairs created by the electric field per 
unit time, the final values of the charges Q v and Q A are 

Qa(t ) = (N-0) + (0-N) = 0, 

Q v (t ) = (N - 0) - (0 - N) = 2N. (7.24) 

Therefore we conclude that the coupling to the electric field produces a violation in the conservation 
of the axial charge per unit time given by AQa ~ e£. This implies that 

eh£, (7.25) 

where we have restored h to make clear that the violation in the conservation of the axial current is 
a quantum effect. At the same time AQ V = guarantees that the vector current remains conserved 
also quantum mechanically, d^Jy = 0. 

We have just studied a two-dimensional example of the Adler-Bell-Jackiw axial anomaly [29]. 
The heuristic analysis presented here can be made more precise by computing the quantity 



C7'- = (o|r[^(x)^(o)]|o> = 




(7.26) 
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The anomaly is given then by d^C^. A careful calculation yields the numerical prefactor missing 
in Eq. (7.25) leading to the result 

P R 

withe 01 = -e 10 = 1. 

The existence of an anomaly in the axial symmetry that we have illustrated in two dimensions 
is present in all even dimensional of space-times. In particular in four dimensions the axial anomaly 
it is given by 

d,r A = -^e^ x F^F aX . (7.28) 

This result has very important consequences in the physics of strong interactions as we will see in 
what follows 



7.2 Chiral symmetry in QCD 

Our knowledge of the physics of strong interactions is based on the theory of Quantum Chromo- 
dynamics (QCD) [31]. This is a nonabelian gauge theory with gauge group SU(iV c ) coupled to a 
number Nf of quarks. These are spin-| particles Q lf labelled by two quantum numbers: color 
i — 1, . . . , N c and flavor / = 1, . . . , Nf. The interaction between them is mediated by the N% — 1 
gauge bosons, the gluons A a , a = 1, . . . , N% — 1. In the real world N c = 3 and the number of 
flavors is six, corresponding to the number of different quarks: up (u), down (d), charm (c), strange 
(s), top (t) and bottom (b). 

For the time being we are going to study a general theory of QCD with N c colors and Nf 
flavors. Also, for reasons that will be clear later we are going to work in the limit of vanishing quark 
masses, rrif — > 0. In this cases the Lagrangian is given by 



(7.29) 



Nf 

/=i 

where the subscripts L and R indicate respectively left and right-handed spinors, R = P±Q f , 
and the field strength F* v and the covariant derivative are respectively defined in Eqs. (4.75) 
and (4.78). Apart from the gauge symmetry, this Lagrangian is also invariant under a global 
\J(Nf) L x\J(Nf) R acting on the flavor indices and defined by 



Qi EAU L )ffQ f L { Qi Q 



L 



V(N f ) L : < \J(N f ) R : <^ (7.30) 

Qi -»> Qi ( Qr Ef,(U R )f f ,Q f R 

with U L ,U R E U(Nf). Actually, since U(A r )=U(l)xSU(A^) this global symmetry group can be 
written as SU(A^/) L x SU(A^/) R x U(1) L x U(l) fl . The abelian subgroup U(1) L x U(l) fl can be now 
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decomposed into their vector U(1) B and axial U(1)a subgroups defined by the transformations 

(7.31) 



Q[ -+ e ia Q f L { Q[ e ia Q f L 

U(l) fl : < U(1)a ■ 



Qr z ia Q f R 



Qr e- fa Q£ 



According to Noether's theorem, associated with these two abelian symmetries we have two con- 
served currents: 



^ = E Q ^ Qf > j a = J2 Q V75 Q f . 

/=1 /=1 



(7.32) 



The conserved charge associated with vector charge J v is actually the baryon number defined as 
the number of quarks minus number of antiquarks. 

The nonabelian part of the global symmetry group S\J(Nf) L xS\J(Nf) R can also be decom- 
posed into its vector and axial subgroups, S\J(Nf) v x S\J(Nf) A , defined by the following transfor- 
mations of the quarks fields 

/' 



Qi ZfWfrQt 



Ql Ef>(u L )ff>Q f L 

SV(N f ) v : <( SU(N f ) A : <j (7.33) 

Q R { Q f R E^VO* 

Again, the application of Noether's theorem shows the existence of the following nonabelian con- 
served charges 



N f Nf 

4" = E Qf-f^frQ* 1 , JiT = E Q'-f'KiT'hrQr. 

/,/'=! /,/'=! 



(7.34) 



To summarize, we have shown that the initial chiral symmetry of the QCD Lagrangian (7.29) can 
be decomposed into its chiral and vector subgroups according to 



U(N f ) L x U(N f ) R = SU(jV,V x S\J(N f ) A x U(1) B x U(1) A . 



(7.35) 



The question to address now is which part of the classical global symmetry is preserved by the 
quantum theory. 

As argued in section 7.1, the conservation of the axial currents J A and J^ 1 can in principle be 
spoiled due to the presence of an anomaly. In the case of the abelian axial current J A the relevant 
quantity is the correlation function 



= (0|T [J^(x)f g : use (x')j 



gauge 



N f 

(0)] |o> = E 
/=1 




(7.36) 



symmetric 



73 



Here j| a ^ gc is the nonabelian conserved current coupling to the gluon field 



.7 gauge — ^ ] Q Q 



(7.37) 



/=1 



where, to avoid confusion with the generators of the global symmetry we have denoted by r a the 
generators of the gauge group SU(iV c ). The anomaly can be read now from d^C^" . If we impose 
Bose symmetry with respect to the interchange of the two outgoing gluons and gauge invariance of 
the whole expression, d^C^ 1 " 7 = = d a C^ ua , we find that the axial abelian global current has an 
anomaly given by 16 



In the case of the nonabelian axial global symmetry S\J(Nf)A the calculation of the anomaly 
is made as above. The result, however, is quite different since in this case we conclude that the non- 
abelian axial current is not anomalous. This can be easily seen by noticing that associated with 
the axial current vertex we have a generator T 1 of SU (Nf), whereas for the two gluon vertices we 
have the generators r a of the gauge group SU(iV c ). Therefore, the triangle diagram is proportional 
to the group-theoretic factor 



which vanishes because the generators of SU (Nf) are traceless. 

From here we would conclude that the nonabelian axial symmetry SU(Nf)A is nonanomalous. 
However this is not the whole story since quarks are charged particles that also couple to photons. 
Hence there is a second potential source of an anomaly coming from the the one-loop triangle 




(7.38) 




J symmetric 



trT 1 tr {r a ,r 6 } = 



(7.39) 



diagram coupling J 1 / 1 to two photons 




Q f 



7 



7 



J symmetric 



(7.40) 



where jg m is the electromagnetic current 




(7.41) 



/=i 



l6 The normalization of the generators T 1 of the global SU(iV/) is given by tr (T I T J ) 
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with qj the electric charge of the /-th quark flavor. A calculation of the diagram in (7.40) shows the 
existence of an Adler-Bell-Jackiw anomaly given by 

e^ x F^F aX , (7.42) 

where F^ is the field strength of the electromagnetic field coupling to the quarks. The only chance 
for the anomaly to cancel is that the factor between brackets in this equation be identically zero. 

Before proceeding let us summarize the results found so far. Because of the presence of 
anomalies the axial part of the global chiral symmetry, SU(A^/)a and U(1)a are not realized quantum 
mechanically in general. We found that U(l)/i is always affected by an anomaly. However, because 
the right-hand side of the anomaly equation (7.38) is a total derivative, the anomalous character of 
J\ does not explain the absence of U(l)/i multiplets in the hadron spectrum, since a new current 
can be constructed which is conserved. In addition, the nonexistence of candidates for a Goldstone 
boson associated with the right quantum numbers indicates that U(1)a is not spontaneously broken 
either, so it has be explicitly broken somehow. This is the so-called U(l)-problem which was solved 
by 't Hooft [32], who showed how the contribution of quantum transitions between vacua with 
topologically nontrivial gauge field configurations (instantons) results in an explicit breaking of this 
symmetry. 

Due to the dynamics of the SU(iV c ) gauge theory the axial nonabelian symmetry is sponta- 
neously broken due to the presence at low energies of a vacuum expectation value for the fermion 
bilinear Q^Qf 

(0\Q f Q f \0) ^ (No summation in /!). (7.43) 

This nonvanishing vacuum expectation value for the quark bilinear actually breaks chiral invari- 
ance spontaneously to the vector subgroup SU(Nf) v , so the only subgroup of the original global 
symmetry that is realized by the full theory at low energy is 

U(N f ) L x \J(N f ) R — ► S\J(N f ) v x U(1) B . (7.44) 

Associated with this breaking a Goldstone boson should appear with the quantum numbers of the 
broken nonabelian current. For example, in the case of QCD the Goldstone bosons associated 
with the spontaneously symmetry breaking induced by the vacuum expectation values (uu), (dd) 
and ((ud — du)) have been identified as the pions n°, tx ± . These bosons are not exactly massless 
because of the nonvanishing mass of the u and d quarks. Since the global chiral symmetry is already 
slightly broken by mass terms in the Lagrangian, the associated Goldstone bosons also have masses 
although they are very light compared to the masses of other hadrons. 

In order to have a better physical understanding of the role of anomalies in the physics of 
strong interactions we particularize now our analysis of the case of real QCD. Since the u and d 
quarks are much lighter than the other four flavors, QCD at low energies can be well described 
by including only these two flavors and ignoring heavier quarks. In this approximation, from our 
previous discussion we know that the low energy global symmetry of the theory is SU(2)y xU(l)s, 



16tt 2 



£(*")//«? 
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where now the vector group SU(2)y is the well-known isospin symmetry. The axial U(l)^ current 
is anomalous due to Eq. (7.38) with Nf = 2. In the case of the nonabelian axial symmetry SU(2) A , 
taking into account that q u = |e and q d = — |e and that the three generators of SU(2) can be written 



in terms of the Pauli matrices as T K = \a K we find 



e 2 



E ( Tl )// 9? = E ( Tl )// 9? = °» E ( T3 )// 9? = ^ (7-45) 

/=«,d /="M /=i*,d 

Therefore is anomalous. 

Physically, the anomaly in the axial current has an important consequence. In the quark 
model, the wave function of the neutral pion ir° is given in terms of those for the u and d quark by 



7T 



0\ 



1 
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(\u)\u) - \d)\d)) . (7.46) 



The isospin quantum numbers of are those of the generator T 3 . Actually the analogy goes 
further since d^J^ is the operator creating a pion 7r° out of the vacuum 

~ d,J 3 /\0). (7.47) 

This leads to the physical interpretation of the triangle diagram (7.40) with as the one loop 
contribution to the decay of a neutral pion into two photons 

7T° — > 2 7 . (7.48) 



This is an interesting piece of physics. In 1967 Sutherland and Veltman [33] presented a 
calculation, using current algebra techniques, according to which the decay of the pion into two 
photons should be suppressed. This however contradicted the experimental evidence that showed 
the existence of such a decay. The way out to this paradox, as pointed out in [29], is the axial 
anomaly. What happens is that the current algebra analysis overlooks the ambiguities associated 
with the regularization of divergences in Quantum Field Theory. A QED evaluation of the triangle 
diagram leads to a divergent integral that has to be regularized somehow. It is in this process that 
the Adler-Bell-Jackiw axial anomaly appears resulting in a nonvanishing value for the 7r° — > 27 
amplitude 17 . 

The existence of anomalies associated with global currents does not necessarily mean diffi- 
culties for the theory. On the contrary, as we saw in the case of the axial anomaly it is its existence 
what allows for a solution of the Sutherland- Veltman paradox and an explanation of the electromag- 
netic decay of the pion. The situation, however, is very different if we deal with local symmetries. 
A quantum mechanical violation of gauge symmetry leads to all kinds of problems, from lack of 
renormalizability to nondecoupling of negative norm states. This is because the presence of an 
anomaly in the theory implies that the Gauss' law constraint V • E a = p a cannot be consistently 

17 An early computation of the triangle diagram for the electromagnetic decay of the pion was made by Steinberger 
in [30]. 
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implemented in the quantum theory. As a consequence states that classically are eliminated by the 
gauge symmetry become propagating fields in the quantum theory, thus spoiling the consistency of 
the theory. 

Anomalies in a gauge symmetry can be expected only in chiral theories where left and right- 
handed fermions transform in different representations of the gauge group. Physically, the most 
interesting example of such theories is the electroweak sector of the Standard Model where, for 
example, left handed fermions transform as doublets under SU(2) whereas right-handed fermions 
are singlets. On the other hand, QCD is free of gauge anomalies since both left- and right-handed 
quarks transform in the fundamental representation of SU(3). 

We consider the Lagrangian 

1 N+ _ N - _ 

C = --F a ^F; u + i X» + * £>1 # , (V-49) 

i=i j=i 

where the chiral fermions ip± transform according to the representations r" ± of the gauge group G 
(a = 1, . . . , dim G). The covariant derivatives are then defined by 

= d^ ± + igA^r^i. (7.50) 

As for global symmetries, anomalies in the gauge symmetry appear in the triangle diagram with one 
axial and two vector gauge current vertices 




(7.51) 

J symmetric 



where gauge vector and axial currents jy M , j"^ 1 are given by 

N+ N— 

i=i j=i 

N + AT_ 

j7 = E^-E^- T -^-' ( 7 - 52 ) 

i=i i=i 

Luckily, we do not have to compute the whole diagram in order to find an anomaly cancellation 
condition, it is enough if we calculate the overall group theoretical factor. In the case of the diagram 
in Eq. (7.51) for every fermion species running in the loop this factor is equal to 

tr [r« ± {r^,r^}], (7.53) 
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where the sign ± corresponds respectively to the generators of the representation of the gauge group 
for the left and right-handed fermions. Hence the anomaly cancellation condition reads 

N+ AT_ 

E tr K+K+^+i] "E tr VlM-^t-}] =0. (7.54) 

Knowing this we can proceed to check the anomaly cancellation in the Standard Model 
SU(3)xSU(2)xU(l). Left handed fermions (both leptons and quarks) transform as doublets with 
respect to the SU(2) factor whereas the right-handed components are singlets. The charge with 
respect to the U(l) part, the hypercharge Y, is determined by the Gell-Mann-Nishijima formula 

Q = T 3 + Y, (7.55) 

where Q is the electric charge of the corresponding particle and T 3 is the eigenvalue with respect 
to the third generator of the SU(2) group in the corresponding representation: T 3 = |o- 3 for the 
doublets and T 3 = for the singlets. For the first family of quarks (u, d) and leptons (e, v e ) we have 
the following field content 



quarks: [ ^ j u a R ,_ d° R 



e 



leptons: I Ue ] e R _i (7.56) 



L,-i 



where a = 1,2,3 labels the color quantum number and the subscript indicates the value of the weak 
hypercharge Y. Denoting the representations of SU(3)xSU(2)xU(l) by (n c , n w )y, with n c and n w 
the representations of SU(3) and SU(2) respectively and Y the hypercharge, the matter content of 
the Standard Model consists of a three family replication of the representations: 

left-handed fermions: (3, 2){ (1,2)^ 

6 2 

(7.57) 

right-handed fermions: (3,1)? (3,l)fi (1, 1)^. 

3 3 

In computing the triangle diagram we have 10 possibilities depending on which factor of the gauge 
group SU(3)xSU(2)xU(l) couples to each vertex: 

SU(3) 3 SU(2) 3 U(l) 3 

SU(3) 2 SU(2) SU(2)U(1) 

SU(3) 2 U(1) SU(2)U(1) 2 
SU(3) SU(2) 2 



SU(3)SU(2)U(1) 
SU(3)U(1) 2 
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It is easy to check that some of them do not give rise to anomalies. For example the anomaly for 
the SU(3) 3 case cancels because left and right-handed quarks transform in the same representation. 
In the case of SU(2) 3 the cancellation happens term by term because of the Pauli matrices identity 

a a a b = § ab + i£ abc a c ^ tQ 

tr [a a {a\ a c }] = 2 (tr a a ) 5 bc = 0. (7.58) 

However the hardest anomaly cancellation condition to satisfy is the one with three U(l)'s. In this 
case the absence of anomalies within a single family is guaranteed by the nontrivial identity 

Yl-^Yl = 3x2x(i) 3 + 2x(-I) 3 -3x(f) 3 -3x(-I) 3 -,- 1)3 



left right 



I) + G 



0. (7.59) 



It is remarkable that the anomaly exactly cancels between leptons and quarks. Notice that this result 
holds even if a right-handed sterile neutrino is added since such a particle is a singlet under the whole 
Standard Model gauge group and therefore does not contribute to the triangle diagram. Therefore 
we see how the matter content of the Standard Model conspires to yield a consistent quantum field 
theory. 

In all our discussion of anomalies we only considered the computation of one-loop diagrams. 
It may happen that higher loop orders impose additional conditions. Fortunately this is not so: the 
Adler-Bardeen theorem [34] guarantees that the axial anomaly only receives contributions from one 
loop diagrams. Therefore, once anomalies are canceled (if possible) at one loop we know that there 
will be no new conditions coming from higher-loop diagrams in perturbation theory. 

The Adler-Bardeen theorem, however, only applies in perturbation theory. It is nonetheless 
possible that nonperturbative effects can result in the quantum violation of a gauge symmetry. This 
is precisely the case pointed out by Witten [35] with respect to the SU(2) gauge symmetry of the 
Standard Model. In this case the problem lies in the nontrivial topology of the gauge group SU(2). 
The invariance of the theory with respect to gauge transformations which are not in the connected 
component of the identity makes all correlation functions equal to zero. Only when the number of 
left-handed SU(2) fermion doublets is even gauge invariance allows for a nontrivial theory. It is 
again remarkable that the family structure of the Standard Model makes this anomaly to cancel 

3x^j +lx^ e j =4 SU(2)-doublets, (7.60) 

where the factor of 3 comes from the number of colors. 



8 Renormalization 

8.1 Removing infinities 

From its very early stages, Quantum Field Theory was faced with infinities. They emerged in the 
calculation of most physical quantities, such as the correction to the charge of the electron due to 
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the interactions with the radiation field. The way these divergences where handled in the 1940s, 
starting with Kramers, was physically very much in the spirit of the Quantum Theory emphasis in 
observable quantities: since the observed magnitude of physical quantities (such as the charge of the 
electron) is finite, this number should arise from the addition of a "bare" (unobservable) value and 
the quantum corrections. The fact that both of these quantities were divergent was not a problem 
physically, since only its finite sum was an observable quantity. To make thing mathematically 
sound, the handling of infinities requires the introduction of some regularization procedure which 
cuts the divergent integrals off at some momentum scale A. Morally speaking, the physical value of 
an observable p h y sicai is given by 

^physical = lim [0(A)bare + AO(A) h ] , (8.1) 

A^oo 

where AO(A) h represents the regularized quantum corrections. 

To make this qualitative discussion more precise we compute the corrections to the elec- 
tric charge in Quantum Electrodynamics. We consider the process of annihilation of an electron- 
positron pair to create a muon-antimuon pair e~e + — > To lowest order in the electric charge 
e the only diagram contributing is 




However, the corrections at order e 4 to this result requires the calculation of seven more diagrams 




In order to compute the renormalization of the charge we consider the first diagram which 
takes into account the first correction to the propagator of the virtual photon interchanged between 
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the pairs due to vacuum polarization. We begin by evaluating 




—irf 
q 2 + ie 



a 



O 



-irf v 
q 2 + ie' 



(8.2) 



where the diagram between brackets is given by 



a 



Qp = ^{q)=i\-ie)\-l) J ^ 



d A k Tr(^ + m e )7 a (^ + ^ + m e )7 /3 



[k 2 — ml + ie] [(k + q) 2 — m 2 e + ie] 



(8.3) 



Physically this diagram includes the correction to the propagator due to the polarization of the vac- 
uum, i.e. the creation of virtual electron-positron pairs by the propagating photon. The momentum 
q is the total momentum of the electron-positron pair in the intermediate channel. 

It is instructive to look at this diagram from the point of view of perturbation theory in non- 
relativistic Quantum Mechanics. In each vertex the interaction consists of the annihilation (resp. 
creation) of a photon and the creation (resp. annihilation) of an electron-positron pair. This can be 
implemented by the interaction Hamiltonian 



H iQt = e J d 6 x^fij}A^ (8.4) 

All fields inside the integral can be expressed in terms of the corresponding creation-annihilation 
operators for photons, electrons and positrons. In Quantum Mechanics, the change in the wave 
function at first order in the perturbation H mt is given by 

| 7 ,m) = |7,m) + 2^ — p p \n) (8.5) 

£/i n — £j n 

n 

and similarly for | 7 , out), where we have denoted symbolically by \n) all the possible states of the 
electron-positron pair. Since these states are orthogonal to | 7 , in) , | 7 , out) , we find torder e 2 

/ • i / ,\ / • i / a i o(7,in|^mt|n) (n| J ff int | 7 / ,out)o , 4 . 

( 7 , m| 7 , out) = 0(7, m|7 , out) + > 7^ r—— — + 0(e). (8.6) 

\tj iTl — £/ n J(£/ out - tL n ) 

Hence, we see that the diagram of Eq. (8.2) really corresponds to the order-e 2 correction to the 
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photon propagator (7, in I7', out) 



7' 




o(7,in|7',out)( 



(7, in|ff in t|n) (n\H- mt \j, out) 
7 " " y " 4" (E in -E n )(E out -E n ) 



Once we understood the physical meaning of the Feynman diagram to be computed we pro- 
ceed to its evaluation. In principle there is no problem in computing the integral in Eq. (8.2) for 
nonzero values of the electron mass. However since here we are going to be mostly interested in 
seeing how the divergence of the integral results in a scale-dependent renormalization of the electric 
charge, we will set m e = 0. This is something safe to do, since in the case of this diagram we are 
not inducing new infrared divergences in taking the electron as massless. Doing some 7-matrices 
gymnastics it is not complicated to show that the polarization tensor II ^(q) defined in Eq. (8.3) can 
be written as 

IW?) = (q\u ~ q^qu) n(g 2 ) (8.8) 

with 

2 _ 4e 2 f d 4 k k 2 + k-q 

[<1 J 3q 2 J (2tt) 4 [k 2 + ie] [{k + q) 2 + ie] ' 1 ' ' 

Although by naive power counting we could conclude that the previous integral is quadratically 
divergent, it can be seen that the quadratic divergence actually cancels leaving behind only a log- 
arithmic one. In order to handle this divergent integral we have to figure out some procedure to 
render it finite. This can be done in several ways, but here we choose to cut the integrals off at a 
high energy scale A, where new physics might be at work, \p\ < A. This gives the result 

IT(g 2 ) ~ ^ log (j^j + finite terms. (8.10) 

If we would send the cutoff to infinity A — > 00 the divergence blows up and something has to be 
done about it. 

If we want to make sense out of this, we have to go back to the physical question that led us 
to compute Eq. (8.2). Our primordial motivation was to compute the corrections to the annihilation 
of two electrons into two muons. Including the correction to the propagator of the virtual photon 
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we have 




4 1 + li l0g (^)]} ( ^ 7V) - (8 ' 11} 



(vAm) + rj a p (v e -f a u e ) — —U(q 2 ) (tvW) 



Now let us imagine that we are performing a e e + — > ji fi + with a center of mass energy /x. From 
the previous result we can identify the effective charge of the particles at this energy scale e(/x) as 



This charge, e(/i), is the quantity that is physically measurable in our experiment. Now we can make 
sense of the formally divergent result (8.11) by assuming that the charge appearing in the classical 
Lagrangian of QED is just a "bare" value that depends on the scale A at which we cut off the theory, 
e = e(A)b a rc- In order to reconcile (8.11) with the physical results (8.12) we must assume that the 
dependence of the bare (unobservable) charge e(A) barc on the cutoff A is determined by the identity 



If we still insist in removing the cutoff, A — > oo we have to send the bare charge to zero e (A) bare — > 
in such a way that the effective coupling has the finite value given by the experiment at the energy 
scale pi. It is not a problem, however, that the bare charge is small for large values of the cutoff, since 
the only measurable quantity is the effective charge that remains finite. Therefore all observable 
quantities should be expressed in perturbation theory as a power series in the physical coupling 
e(fx) 2 and not in the unphysical bare coupling e(A) bar e- 

8.2 The beta-function and asymptotic freedom 

We can look at the previous discussion, an in particular Eq. (8.13), from a different point of view. 
In order to remove the ambiguities associated with infinities we have been forced to introduce a 
dependence of the coupling constant on the energy scale at which a process takes place. From the 
expression of the physical coupling in terms of the bare charge (8.13) we can actually eliminate 
the cutoff A, whose value after all should not affect the value of physical quantities. Taking into 
account that we are working in perturbation theory in e(/i) 2 , we can express the bare charge e(A)^ arc 
in terms of e(/x) 2 as 




(8.12) 




(8.13) 



e(A) 2 = e(/i) 2 1 + 



e(/x) 2 



log(^)]+0[e(^) 6 ]. 



(8.14) 



12tt 2 
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This expression allow us to eliminate all dependence in the cutoff in the expression of the effective 
charge at a scale /i by replacing e(A) bare in Eq. (8.13) by the one computed using (8.14) at a given 
reference energy scale /jl 



e(» 2 = e(nof 



e(/io) 2 , ( ^ 



1 + T2^ l0g V^;j 



(8.15) 



From this equation we can compute, at this order in perturbation theory, the effective value of 
the coupling constant at an energy //, once we know its value at some reference energy scale /j, . In 
the case of the electron charge we can use as a reference Thompson's scattering at energies of the 
order of the electron mass m e ~ 0.5 MeV, at where the value of the electron charge is given by the 
well known value 

eK) 2 ~ (8.16) 

With this we can compute e(/i) 2 at any other energy scale applying Eq. (8.15), for example at the 
electron mass /j, = m e ~ 0.5 MeV. However, in computing the electromagnetic coupling constant 
at any other scale we must take into account the fact that other charged particles can run in the loop 
in Eq. (8.11). Suppose, for example, that we want to calculate the fine structure constant at the 
mass of the Z°-boson \x = M z = 92 GeV. Then we should include in Eq. (8.15) the effect of other 
fermionic Standard Model fields with masses below M z . Doing this, we find 18 



e{M z f = e(m e ) 2 




(8.17) 



where is the charge in units of the electron charge of the i-th fermionic species running in the 
loop and we sum over all fermions with masses below the mass of the Z° boson. This expression 
shows how the electromagnetic coupling grows with energy. However, in order to compare with the 
experimental value of e(M z ) 2 it is not enough with including the effect of fermionic fields, since 
also the W ± bosons can run in the loop (M w < M z ). Taking this into account, as well as threshold 
effects, the value of the electron charge at the scale M z is found to be [36] 

e(M z f ~ — . (8.18) 
v ; 128.9 

This growing of the effective fine structure constant with energy can be understood heuris- 
tically by remembering that the effect of the polarization of the vacuum shown in the diagram of 
Eq. (8.2) amounts to the creation of a plethora of electron-positron pairs around the location of the 
charge. These virtual pairs behave as dipoles that, as in a dielectric medium, tend to screen this 
charge and decreasing its value at long distances (i.e. lower energies). 



18 In the first version of these notes the argument used to show the growing of the electromagnetic coupling constant 
could have led to confusion to some readers. To avoid this potential problem we include in the equation for the running 
coupling e(p) 2 the contribution of all fermions with masses below Mz- We thank Lubos Motl for bringing this issue to 
our attention. 
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The variation of the coupling constant with energy is usually encoded in Quantum Field The- 
ory in the beta function defined by 



0(9) = V 



dg_ 



(8.19) 



In the case of QED the beta function can be computed from Eq. (8.15) with the result 



/5(e) QED 




(8.20) 



The fact that the coefficient of the leading term in the beta-function is positive f3 = ^ > gives 
us the overall behavior of the coupling as we change the scale. Eq. (8.20) means that, if we start 
at an energy where the electric coupling is small enough for our perturbative treatment to be valid, 
the effective charge grows with the energy scale. This growing of the effective coupling constant 
with energy means that QED is infrared safe, since the perturbative approximation gives better and 
better results as we go to lower energies. Actually, because the electron is the lighter electrically 
charged particle and has a finite nonvanishing mass the running of the fine structure constant stops 
at the scale m e in the well-known value Would other charged fermions with masses below m e 
be present in Nature, the effective value of the fine structure constant in the interaction between 
these particles would run further to lower values at energies below the electron mass. 

On the other hand if we increase the energy scale e(/x) 2 grows until at some scale the coupling 
is of order one and the perturbative approximation breaks down. In QED this is known as the 
problem of the Landau pole but in fact it does not pose any serious threat to the reliability of QED 
perturbation theory: a simple calculation shows that the energy scale at which the theory would 
become strongly coupled is A La iidau — 10 277 GeV. However, we know that QED does not live that 
long! At much lower scales we expect electromagnetism to be unified with other interactions, and 
even if this is not the case we will enter the uncharted territory of quantum gravity at energies of the 
order of 10 19 GeV. 

So much for QED. The next question that one may ask at this stage is whether it is possible 
to find quantum field theories with a behavior opposite to that of QED, i.e. such that they become 
weakly coupled at high energies. This is not a purely academic question. In the late 1960s a series of 
deep-inelastic scattering experiments carried out at SLAC showed that the quarks behave essentially 
as free particles inside hadrons. The apparent problem was that no theory was known at that time 
that would become free at very short distances: the example set by QED seem to be followed by 
all the theories that were studied. This posed a very serious problem for Quantum Field Theory 
as a way to describe subnuclear physics, since it seemed that its predictive power was restricted to 
electrodynamics but failed miserably when applied to describe strong interactions. 

Nevertheless, this critical time for Quantum Field Theory turned out to be its finest hour. In 
1973 David Gross and Frank Wilczek [37] and David Politzer [38] showed that nonabelian gauge 
theories can actually display the required behavior. For the QCD Lagrangian in Eq. (7.29) the beta 
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Fig. 14: Beta function for a hypothetical theory with three fixed points g*, g\ and g\. A perturbative analysis 
would capture only the regions shown in the boxes. 

function is given by 19 



In particular, for real QCD (N c = 3, Nf = 6) we have that (3(g) = < 0. This means that 

for a theory that is weakly coupled at an energy scale fi the coupling constant decreases as the 
energy increases /i — > oo. This explain the apparent freedom of quarks inside the hadrons: when 
the quarks are very close together their effective color charge tend to zero. This phenomenon is 
called asymptotic freedom. 

Asymptotic free theories display a behavior that is opposite to that found above in QED. At 
high energies their coupling constant approaches zero whereas at low energies they become strongly 
coupled (infrared slavery). This features are at the heart of the success of QCD as a theory of strong 
interactions, since this is exactly the type of behavior found in quarks: they are quasi-free particles 
inside the hadrons but the interaction potential potential between them increases at large distances. 

Although asymptotic free theories can be handled in the ultraviolet, they become extremely 
complicated in the infrared. In the case of QCD it is still to be understood (at least analytically) how 
the theory confines color charges and generates the spectrum of hadrons, as well as the breaking of 
the chiral symmetry (7.43). 

In general, the ultraviolet and infrared properties of a theory are controlled by the fixed points 
of the beta function, i.e. those values of the coupling constant g for which it vanishes 



The expression of the beta function of QCD was also known to 't Hooft [39]. There are even earlier computations 
in the russian literature [40]. 




(8.21) 



(3(g*) = 0. 



(8.22) 
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Using perturbation theory we have seen that for both QED and QCD one of such fixed points 
occurs at zero coupling, g* = 0. However, our analysis also showed that the two theories present 
radically different behavior at high and low energies. From the point of view of the beta function, 
the difference lies in the energy regime at which the coupling constant approaches its critical value. 
This is in fact governed by the sign of the beta function around the critical coupling. 

We have seen above that when the beta function is negative close to the fixed point (the case of 
QCD) the coupling tends to its critical value, g* = 0, as the energy is increased. This means that the 
critical point is ultraviolet stable, i.e. it is an attractor as we evolve towards higher energies. If, on 
the contrary, the beta function is positive (as it happens in QED) the coupling constant approaches 
the critical value as the energy decreases. This is the case of an infrared stable fixed point. 

This analysis that we have motivated with the examples of QED and QCD is completely 
general and can be carried out for any quantum field theory. In Fig. 14 we have represented the 
beta function for a hypothetical theory with three fixed points located at couplings g\, g* 2 and g\. 
The arrows in the line below the plot represent the evolution of the coupling constant as the energy 
increases. From the analysis presented above we see that g* — and g% are ultraviolet stable fixed 
points, while the fixed point g\ is infrared stable. 

In order to understand the high and low energy behavior of a quantum field theory it is then 
crucial to know the structure of the beta functions associated with its couplings. This can be a very 
difficult task, since perturbation theory only allows the study of the theory around "trivial" fixed 
points, i.e. those that occur at zero coupling like the case of g\ in Fig. 14. On the other hand, any 
"nontrivial" fixed point occurring in a theory (like g\ and g%) cannot be captured in perturbation 
theory and requires a full nonperturbative analysis. 

The moral to be learned from our discussion above is that dealing with the ultraviolet di- 
vergences in a quantum field theory has the consequence, among others, of introducing an energy 
dependence in the measured value of the coupling constants of the theory (for example the elec- 
tric charge in QED). This happens even in the case of renormalizable theories without mass terms. 
These theories are scale invariant at the classical level because the action does not contain any di- 
mensionful parameter. In this case the running of the coupling constants can be seen as resulting 
from a quantum breaking of classical scale invariance: different energy scales in the theory are dis- 
tinguished by different values of the coupling constants. Remembering what we learned in Section 
7, we conclude that classical scale invariance is an anomalous symmetry. One heuristic way to see 
how the conformal anomaly comes about is to notice that the regularization of an otherwise scale 
invariant field theory requires the introduction of an energy scale (e.g. a cutoff). This breaking of 
scale invariance cannot be restored after renormalization. 

Nevertheless, scale invariance is not lost forever in the quantum theory. It is recovered at the 
fixed points of the beta function where, by definition, the coupling does not run. To understand how 
this happens we go back to a scale invariant classical field theory whose field <f>(x) transform under 
coordinate rescalings as 

a? — -» Ax", 4>{x) — > \- A (f)(X- 1 x), (8.23) 
where A is called the canonical scaling dimension of the field. An example of such a theory is a 
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massless 4 theory in four dimensions 



£ = i<9 M 0^0-|0 4 , (8.24) 

where the scalar field has canonical scaling dimension A = 1. The Lagrangian density transforms 
as 

C — >• A~ 4 £[0] (8.25) 

and the classical action remains invariant 20 . 

We look at the free theory g = for a moment. Now there are no divergences and all corre- 
lation functions can be exactly computed. In particular we consider the momentum space n-point 
correlation function 

G (p 1 ,...,p n )(2ir)^\p 1 + ...+p n ) 

= J rf 4 x 1 ...rf 4 x n e^- :E1+ - + ^(O|T[0 o ( a ; 1 )...0 o (x n )]|O), (8.26) 

where by (j>o(x) we denote the free field operator. Applying the rescaling (8.23) we find the follow- 
ing transformation for the correlation function 

G (pi, ...,?„)—► X 4 ^- 1) - nA G (X Pl , \p n ). (8.27) 

For the free theory the only relevant correlation function is the two-point function, where we have 
(remember that we are dealing with a massless theory) 

G (P 2 ) = - 2 ^ A 2 G (AV) = - r (8.28) 

The transformation of any other correlation function follows from this result and Wick's theorem, 
that allows to write any the 2n-correlation function as sum of products of n 2-point correlation 
functions (correlations functions with an odd number of fields are identically zero). 

We turn to the interacting theory. Things now get much more complicated, since correlation 
functions cannot be exactly computed in general. However, when the theory sits at the critical 
coupling we can use a few useful facts. For example, since the critical theory is scale invariant, it 
should either contain only massless one-particle states or have continuous spectrum. To keep the 
argument simple, we consider the first possibility. Hence, the exact two-point function should have 
a pole at p 2 = 0, and close to this pole the correlation function has the form 

G(p 2 ;/,)-^, (8-29) 
p 1 



20 In a D-dimensional theory the canonical scaling dimensions of the fields coincide with its engineering dimension: 
A = ^=1 for bosonic fields and A = for fermionic ones. For a Lagrangian with no dimensionful parameters 
classical scale invariance follows then from dimensional analysis. 
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where Z ([/,), called the field renormalization, depends on the scale. The anomalous dimension 7(g) 
is then defined by the equation 

l(g) = l^\ogZ. (8.30) 

This new function is the analog of the beta function (8.19) for the field renormalization Z(fj). 
Moreover, at the critical point <?(//) = g* and the anomalous dimension is independent of the energy, 
7* — lid*)- m this case Eq. (8.30) can be integrated to give 

/ \ 27* 

Z(n) = Z Q UL\ , (8.31) 

where Zq and fj, are some reference values. Then, we find that the two-point function at the critical 
point is invariant under the rescaling 

G(p 2 ; 11) — >■ ^^^(AV; A/i). (8.32) 

Here we have presented a rather sketchy and heuristic argument. A more thorough analysis 
(using for example the Callan-Symanzik equation [1-10]) shows that at the critical point all n-point 
correlation functions are invariant under the rescaling 

G(p u . . . ,p n ; n) — > \« n -V- n W)G{\ Pl , . . . , \p n , A/i). (8.33) 

Comparing (8.32) and (8.33) with (8.27) we see that this invariance is analogous to the one of the 
free (scale invariant) theory. Now, however, the fields transform under rescalings with an anomalous 
scaling dimension given by 

A anom = A + 7*, (8.34) 

with A the canonical scaling dimension of the corresponding field. This justifies the name given to 
the function 7(g) defined in Eq. (8.30). Notice, however, that strictly speaking 7(g) only represents 
an anomalous dimension for the theory at the critical coupling g*. 

The previous discussion clarifies a little bit the high-energy properties of an asymptotically 
free theories like QCD. The fact that the fixed point occurs at zero coupling might give the wrong 
impression that the theory at the critical point is just the one obtains by setting g — in the action. 
Life, however, is more complicated than that. What we have seen above shows that although the 
critical theory is a free scale invariant field theory, the fields have anomalous scaling dimensions 
which are different from the ones of the "naive" free theory. These anomalous dimensions carry the 
dynamical information about the high-energy behavior of the asymptotically free theory. 

8.3 The renormalization group 

In spite of its successes, the renormalization procedure presented above can be seen as some kind 
of prescription or recipe to get rid of the divergences in an ordered way. This discomfort about 
renormalization was expressed in occasions by comparing it with "sweeping the infinities under 



89 



Fig. 15: Systems of spins in a two-dimensional square lattice. 

the rug". However thanks to Ken Wilson to a large extent [41] the process of renormalization is 
now understood in a very profound way as a procedure to incorporate the effects of physics at high 
energies by modifying the value of the parameters that appear in the Lagrangian. 

Statistical mechanics. Wilson's ideas are both simple and profound and consist in thinking 
about Quantum Field Theory as the analog of a thermodynamical description of a statistical system. 
To be more precise, let us consider an Ising spin system in a two-dimensional square lattice as the 
one depicted in Fig 15. In terms of the spin variables s; = ±|, where % labels the lattice site, the 
Hamiltonian of the system is given by 

H = -jJ^SiSj, (8.35) 

(id) 

where indicates that the sum extends over nearest neighbors and J is the coupling constant 
between neighboring spins (here we consider that there is no external magnetic field). The starting 
point to study the statistical mechanics of this system is the partition function defined as 

3 = 5%-** (8.36) 

{*} 

where the sum is over all possible configurations of the spins and (3 = ^ is the inverse temperature. 
For J > the Ising model presents spontaneous magnetization below a critical temperature T c , 
in any dimension higher than one. Away from this temperature correlations between spins decay 
exponentially at large distances 

(siSj) ~ e~¥, (8.37) 

with \xij\ the distance between the spins located in the 2-th and j-th sites of the lattice. This ex- 
pression serves as a definition of the correlation length £ which sets the characteristic length scale 
at which spins can influence each other by their interaction through their nearest neighbors. 

Suppose now that we are interested in a macroscopic description of this spin system. We can 
capture the relevant physics by integrating out somehow the physics at short scales. A way in which 



90 




Fig. 16: Decimation of the spin lattice. Each block in the upper lattice is replaced by an effective spin 
computed according to the rule (8.39). Notice also that the size of the lattice spacing is doubled in the 
process. 



this can be done was proposed by Leo Kadanoff [42] and consists in dividing our spin system in 
spin-blocks like the ones showed in Fig 16. Now we can construct another spin system where each 
spin-block of the original lattice is replaced by an effective spin calculated according to some rule 
from the spins contained in each block B a 



{si-.ie B a } 



(i) 



(8.38) 



For example we can define the effective spin associated with the block B a by taking the majority 
rule with an additional prescription in case of a draw 



(i) 




(8.39) 



where we have used the sign function, sign(x) = rjjr, with the additional definition sgn(O) = 1. This 
procedure is called decimation and leads to a new spin system with a doubled lattice space. 

The idea now is to rewrite the partition function (8.36) only in terms of the new effective spins 
s a ■ Then we start by splitting the sum over spin configurations into two nested sums, one over the 
spin blocks and a second one over the spins within each block 



Z 



{s} 



-PH[s 



E E * 

{s(l)} {s&B a } 



s a {1) ~ si g n 




-PH[si] 



(8.40) 
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The interesting point now is that the sum over spins inside each block can be written as the expo- 
nential of a new effective Hamiltonian depending only on the effective spins, [s^] 



E ' 

{seB a } 



s (1) 



sign 



\ieB a 



e -pH[si] = e -PHW[ S V]_ ( g_ 41 ) 



The new Hamiltonian is of course more complicated 

jy<i> = -jW£ a < 1 > a < 1 > + ... (8.42) 

(iJ) 

where the dots stand for other interaction terms between the effective block spins. This new terms 
appear because in the process of integrating out short distance physics we induce interactions be- 
tween the new effective degrees of freedom. For example the interaction between the spin block 
variables sf^ will in general not be restricted to nearest neighbors in the new lattice. The impor- 
tant point is that we have managed to rewrite the partition function solely in terms of this new 
(renormalized) spin variables s ^ interacting through a new Hamiltonian 

Z=J2 e~^''\ (8.43) 

Let us now think about the space of all possible Hamiltonians for our statistical system includ- 
ing all kinds of possible couplings between the individual spins compatible with the symmetries of 
the system. If denote by 1Z the decimation operation, our previous analysis shows that 1Z defines a 
map in this space of Hamiltonians 

K : H ->■ ff (8.44) 

At the same time the operation 1Z replaces a lattice with spacing a by another one with double 
spacing 2a. As a consequence the correlation length in the new lattice measured in units of the 
lattice spacing is divided by two, 1Z : £ — > |. 

Now we can iterate the operation 1Z an indefinite number of times. Eventually we might reach 
a Hamiltonian ff* that is not further modified by the operation TZ 

H if CD A H (2) A ... A if*. (8.45) 

The fixed point Hamiltonian ff* is scale invariant because it does not change as 1Z is performed. 
Notice that because of this invariance the correlation length of the system at the fixed point do not 
change under 1Z. This fact is compatible with the transformation £ — > | only if £ = or £ = oo. 
Here we will focus in the case of nontrivial fixed points with infinite correlation length. 

The space of Hamiltonians can be parametrized by specifying the values of the coupling 
constants associated with all possible interaction terms between individual spins of the lattice. If 
we denote by O a [si] these (possibly infinite) interaction terms, the most general Hamiltonian for the 
spin system under study can be written as 

oo 

H[si] =5> a a [si], ( 8 - 46 ) 

a=l 
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where are the coupling constants for the corresponding operators. These constants can be 

thought of as coordinates in the space of all Hamiltonians. Therefore the operation 1Z defines a 
transformation in the set of coupling constants 



For example, in our case we started with a Hamiltonian in which only one of the coupling constants 
is different from zero (say X L = — J). As a result of the decimation Ai = — J — > —J^ while some 
of the originally vanishing coupling constants will take a nonzero value. Of course, for the fixed 
point Hamiltonian the coupling constants do not change under the scale transformation 1Z. 

Physically the transformation 1Z integrates out short distance physics. The consequence for 
physics at long distances is that we have to replace our Hamiltonian by a new one with different 
values for the coupling constants. That is, our ignorance of the details of the physics going on 
at short distances result in a renormalization of the coupling constants of the Hamiltonian that 
describes the long range physical processes. It is important to stress that although TZ is sometimes 
called a renormalization group transformation in fact this is a misnomer. Transformations between 
Hamiltonians defined by 1Z do not form a group: since these transformations proceed by integrating 
out degrees of freedom at short scales they cannot be inverted. 

In statistical mechanics fixed points under renormalization group transformations with £ = oo 
are associated with phase transitions. From our previous discussion we can conclude that the space 
of Hamiltonians is divided in regions corresponding to the basins of attraction of the different fixed 
points. We can ask ourselves now about the stability of those fixed points. Suppose we have a 
statistical system described by a fixed-point Hamiltonian H* and we perturb it by changing the 
coupling constant associated with an interaction term O. This is equivalent to replace H* by the 
perturbed Hamiltonian 



where 5X is the perturbation of the coupling constant corresponding to O (we can also consider per- 
turbations in more than one coupling constant). At the same time thinking of the A a 's as coordinates 
in the space of all Hamiltonians this corresponds to moving slightly away from the position of the 
fixed point. 

The question to decide now is in which direction the renormalization group flow will take the 
perturbed system. Working at first order in 5\ there are three possibilities: 

- The renormalization group flow takes the system back to the fixed point. In this case the 
corresponding interaction O is called irrelevant. 

- 1Z takes the system away from the fixed point. If this is what happens the interaction is called 
relevant. 

- It is possible that the perturbation actually does not take the system away from the fixed point 
at first order in 5\. In this case the interaction is said to be marginal and it is necessary to go 
to higher orders in SX in order to decide whether the system moves to or away the fixed point, 
or whether we have a family of fixed points. 




(8.47) 



H = H* + 5\0, 



(8.48) 
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Fig. 17: Example of a renormalization group flow. 



Therefore we can picture the action of the renormalization group transformation as a flow in 
the space of coupling constants. In Fig. 17 we have depicted an example of such a flow in the case 
of a system with two coupling constants Ai and A 2 . In this example we find two fixed points, one at 
the origin O and another at F for a finite value of the couplings. The arrows indicate the direction 
in which the renormalization group flow acts. The free theory at A x = A 2 = is a stable fix point 
since any perturbation 5\i, 5\ 2 > makes the theory flow back to the free theory at long distances. 
On the other hand, the fixed point F is stable with respect to certain type of perturbations (along the 
line with incoming arrows) whereas for any other perturbations the system flows either to the free 
theory at the origin or to a theory with infinite values for the couplings. 

Quantum field theory. Let us see now how these ideas of the renormalization group apply to 
Field Theory. Let us begin with a quantum field theory defined by the Lagrangian 



where £o[0j is the kinetic part of the Lagrangian and g,i are the coupling constants associated 
with the operators Oi [<f> a ] . In order to make sense of the quantum theory we introduce a cutoff in 
momenta A. In principle we include all operators Oi compatible with the symmetries of the theory. 

In section 8.2 we saw how in the cases of QED and QCD, the value of the coupling constant 
changed with the scale from its value at the scale A. We can understand now this behavior along the 
lines of the analysis presented above for the Ising model. If we would like to compute the effective 
dynamics of the theory at an energy scale /i < A we only have to integrate out all physical models 
with energies between the cutoff A and the scale of interest fi. This is analogous to what we did in 
the Ising model by replacing the original spins by the block spins. In the case of field theory the 




(8.49) 
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effective action S[(f> a , ji\ at scale ji can be written in the language of functional integration as 



J^<P<A a 



(8.50) 



Here S[<f> a , A] is the action at the cutoff scale 



S[4> a , A] = d 4 x I £ O [0J + $>( A )<WJ 



(8.51) 



and the functional integral in Eq. (8.50) is carried out only over the field modes with momenta in 
the range p, < p < A. The action resulting from integrating out the physics at the intermediate 
scales between A and p depends not on the original field variable 4> a but on some renormalized 
field (f)' a . At the same time the couplings gi(p) differ from their values at the cutoff scale <?j(A). 
This is analogous to what we learned in the Ising model: by integrating out short distance physics 
we ended up with a new Hamiltonian depending on renormalized effective spin variables and with 
renormalized values for the coupling constants. Therefore the resulting effective action at scale p 
can be written as 



This Wilsonian interpretation of renormalization sheds light to what in section 8.1 might have 
looked just a smart way to get rid of the infinities. The running of the coupling constant with 
the energy scale can be understood now as a way of incorporating into an effective action at scale p 
the effects of field excitations at higher energies E > p. 

As in statistical mechanics there are also quantum field theories that are fixed points of the 
renormalization group flow, i.e. whose coupling constants do not change with the scale. We have 
encountered them already in Section 8.2 when studying the properties of the beta function. The 
most trivial example of such theories are massless free quantum field theories, but there are also 
examples of four-dimensional interacting quantum field theories which are scale invariant. Again 
we can ask the question of what happens when a scale invariant theory is perturbed with some 
operator. In general the perturbed theory is not scale invariant anymore but we may wonder whether 
the perturbed theory flows at low energies towards or away the theory at the fixed point. 

In quantum field theory this can be decided by looking at the canonical dimension d[0] of the 
operator 0[<f> a ] used to perturb the theory at the fixed point. In four dimensions the three possibilities 
are defined by: 

- d[0] > 4: irrelevant perturbation. The running of the coupling constants takes the theory 
back to the fixed point. 

- d[0] < 4: relevant perturbation. At low energies the theory flows away from the scale- 
invariant theory. 

- d[0\ = 4: marginal deformation. The direction of the flow cannot be decided only on dimen- 
sional grounds. 




(8.52) 
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As an example, let us consider first a massless fermion theory perturbed by a four-fermion 
interaction term 



£ = ##-^(^) 2 - (8-53) 

This is indeed a perturbation by an irrelevant operator, since in four-dimensions [ip] = §. Inter- 
actions generated by the extra term are suppressed at low energies since typically their effects are 
weighted by the dimensionless factor jfa, where E is the energy scale of the process. This means 
that as we try to capture the relevant physics at lower and lower energies the effect of the pertur- 
bation is weaker and weaker rendering in the infrared limit E — >■ again a free theory. Hence, the 
irrelevant perturbation in (8.53) makes the theory flow back to the fixed point. 

On the other hand relevant operators dominate the physics at low energies. This is the case, 
for example, of a mass term. As we lower the energy the mass becomes more important and once 
the energy goes below the mass of the field its dynamics is completely dominated by the mass term. 
This is, for example, how Fermi's theory of weak interactions emerges from the Standard Model at 
energies below the mass of the W ± boson 




At energies below Mw = 80.4 GeV the dynamics of the W + boson is dominated by its mass term 
and therefore becomes nonpropagating, giving rise to the effective four-fermion Fermi theory. 

To summarize our discussion so far, we found that while relevant operators dominate the dy- 
namics in the infrared, taking the theory away from the fixed point, irrelevant perturbations become 
suppressed in the same limit. Finally we consider the effect of marginal operators. As an example 
we take the interaction term in massless QED, O = ip^ip A^. Taking into account that in d = 4 the 
dimension of the electromagnetic potential is [A^] = 1 the operator O is a marginal perturbation. 
In order to decide whether the fixed point theory 

C Q = ~F lw F^ + i^p^ (8.54) 

is restored at low energies or not we need to study the perturbed theory in more detail. This we have 
done in section 8.1 where we learned that the effective coupling in QED decreases at low energies. 
Then we conclude that the perturbed theory flows towards the fixed point in the infrared. 

As an example of a marginal operator with the opposite behavior we can write the Lagrangian 
for a SU(iV c ) gauge theory, C = -\F* v F a » v , as 

C = -^(dpAZ-d v A*)(d»A av -drA a »)-4gf abc A°A!ldi i A cv 

+ g 2 f abc f ade A b ll AlA dll A ev = Co + O g , (8.55) 
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i.e. a marginal perturbation of the free theory described by C , which is obviously a fixed point 
under renormalization group transformations. Unlike the case of QED we know that the full theory 
is asymptotically free, so the coupling constant grows at low energies. This implies that the operator 
O g becomes more and more important in the infrared and therefore the theory flows away the fixed 
point in this limit. 

It is very important to notice here that in the Wilsonian view the cutoff is not necessarily 
regarded as just some artifact to remove infinities but actually has a physical origin. For example in 
the case of Fermi's theory of /3-decay there is a natural cutoff A = M w at which the theory has to 
be replaced by the Standard Model. In the case of the Standard Model itself the cutoff can be taken 
at Planck scale A ~ 10 19 GeV or the Grand Unification scale A ~ 10 16 GeV, where new degrees 
of freedom are expected to become relevant. The cutoff serves the purpose of cloaking the range of 
energies at which new physics has to be taken into account. 

Provided that in the Wilsonian approach the quantum theory is always defined with a physical 
cutoff, there is no fundamental difference between renormalizable and nonrenormalizable theories. 
Actually, a renormalizable field theory, like the Standard Model, can generate nonrenormalizable 
operators at low energies such as the effective four-fermion interaction of Fermi's theory. They 
are not sources of any trouble if we are interested in the physics at scales much below the cutoff, 
EcA, since their contribution to the amplitudes will be suppressed by powers of 



9 Special topics 

9.1 Creation of particles by classical fields 

Particle creation by a classical source. In a free quantum field theory the total number of particles 
contained in a given state of the field is a conserved quantity. For example, in the case of the quantum 
scalar field studied in section 3 we have that the number operator commutes with the Hamiltonian 



n = 



d 3 k 
J2nf 



o?{k)a{k), [H,n] = 0. 



(9.1) 



This means that any states with a well-defined number of particle excitations will preserve this 
number at all times. The situation, however, changes as soon as interactions are introduced, since in 
this case particles can be created and/or destroyed as a result of the dynamics. 

Another case in which the number of particles might change is if the quantum theory is cou- 
pled to a classical source. The archetypical example of such a situation is the Schwinger effect, 
in which a classical strong electric field produces the creation of electron-positron pairs out of the 
vacuum. However, before plunging into this more involved situation we can illustrate the relevant 
physics involved in the creation of particles by classical sources with the help of the simplest ex- 
ample: a free scalar field theory coupled to a classical external source J(x). The action for such a 
theory can be written as 



S 



d x 



I 771 

-d^{x)d"ct>{x) - —<P{xf + J(x)cf>(x) 



(9.2) 



where J(x) is a real function of the coordinates. Its identification with a classical source is obvious 
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once we calculate the equations of motion 



(V 2 + m 2 ) <f)(x) = J{x) 



(9.3) 



Our plan is to quantize this theory but, unlike the case analyzed in section 3, now the presence of 
the source J(x) makes the situation a bit more involved. The general solution to the equations of 
motion can be written in terms of the retarded Green function for the Klein-Gordon equation as 



(9.4) 



(f)(x) = <f) (x) + i J (Fx' G R (x - x')J(x'), 

where 4>q(x) is a general solution to the homogeneous equation and 

d A k i _ ik . x . „,. N f d 3 k 1 



G R (t,x) = 



(2tt) 4 k 2 - m 2 



= i6(t) 



(2tt)3 2u k 



e -iu) k t+k-x _ ^jjj k t—ip-x 



(9.5) 



with 9(x) the Heaviside step function. The integration contour to evaluate the integral over p° 
surrounds the poles at p° = ±u k from above. Since Gn(t,x) = for t < 0, the function <f>o(x) 
corresponds to the solution of the field equation at t — > — oo, before the interaction with the external 
source 21 

To make the argument simpler we assume that J(x) is switched on at t = 0, and only last for 
a time r, that is 



J(t,x) = ift<0ort>r. 



(9.6) 



We are interested in a solution of (9.3) for times after the external source has been switched off, 
t > r. In this case the expression (9.5) can be written in terms of the Fourier modes J(u, k) of the 
source as 



(f>(t,x) = (f>o(x) + i 



d 3 k 1 



J(u k , k)e 



-iuj^t+ik-x 



J(u k , k)*e 



* iui^t—ik-x 



(9.7) 



On the other hand, the general solution <j>o(x) has been already computed in Eq. (3.53). Combining 
this result with Eq. (9.7) we find the following expression for the late time general solution to the 
Klein-Gordon equation in the presence of the source 



(f)(t,x) 



d 3 k 



1 



+ 



(2vr) 3 j%T k 



a*(k) - 



a(k) + 



J(u k , k)* 



J(uj k , k) 



-iiOfrt+ik-x 



D iLo^t—ik-x 



(9.8) 



We should not forget that this is a solution valid for times t > r, i.e. once the external source has 
been disconnected. On the other hand, for t < we find from Eqs. (9.4) and (9.5) that the general 
solution is given by Eq. (3.53). 



21 We could have taken instead the advanced propagator Ga (x) in which case <j>o (x) would correspond to the solution 
to the equation at large times, after the interaction with J (a;). 
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Now we can proceed to quantize the theory. The conjugate momentum ir(x) = d <j)(x) can 
be computed from Eqs. (3.53) and (9.8). Imposing the canonical equal time commutation relations 
(3.50) we find that a(k), cfi(k) satisfy the creation- annihilation algebra (3.27). From our previous 
calculation we find that for t > r the expansion of the operator <f>(x) in terms of the creation- 
annihilation operators a(k), a\k) can be obtained from the one for t < by the replacement 

a(k) — > (3{k) = a(k) + —^=J(u k ,k), 

a\k) — ► ^(k) = a^k)--^=J(uj k ,kr. (9.9) 

Actually, since J(co k , k) is a c-number, the operators f3(k), P\k) satisfy the same algebra as a(k), 
a^(k) and therefore can be interpreted as well as a set of creation-annihilation operators. This means 
that we can define two vacuum states, |0_), |0 + ) associated with both sets of operators 



a(fc)|0_) = 
/3$|0 + ) = 



V k. (9.10) 



— * — * 

For an observer at t < 0, a(k) and a(k) are the natural set of creation-annihilation operators 
in terms of which to expand the field operator <f)(x). After the usual zero-point energy subtraction 
the Hamiltonian is given by 

H { - ] = J d 3 ku k a\k)a{k) (9.11) 

and the ground state of the spectrum for this observer is the vacuum |0_). At the same time, a 
second observer at t > t will also see a free scalar quantum field (the source has been switched 
off at t — r) and consequently will expand in terms of the second set of creation- annihilation 
operators (5(k), (5\k). In terms of this operators the Hamiltonian is written as 

H {+) = J d 3 ku k ft(k)[3(k). (9.12) 

Then for this late-time observer the ground state of the Hamiltonian is the second vacuum state |0 + ). 

In our analysis we have been working in the Heisenberg picture, where states are time- 
independent and the time dependence comes in the operators. Therefore the states of the theory 
are globally defined. Suppose now that the system is in the "in" ground state |0_). An observer at 
t < will find that there are no particles 

C-) |0_> = 0. (9.13) 



n 



However the late-time observer will find that the state |0_) contains an average number of particles 
given by 

d 3 k 1 - 2 

(9.14) 



(tL|^)|0_) = /^|JK5) 
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Fig. 18: Pair creation by a electric field in the Dirac sea picture. 



Moreover, |0_) is no longer the ground state for the "out" observer. On the contrary, this state have 



The key to understand what is going on here lies in the fact that the external source breaks 
the invariance of the theory under space-time translations. In the particular case we have studied 
here where J(x) has support over a finite time interval < t < r, this implies that the vacuum is 
not invariant under time translations, so observers at different times will make different choices of 
vacuum that will not necessarily agree with each other. This is clear in our example. An observer in 
t < t will choose the vacuum to be the lowest energy state of her Hamiltonian, |0_). On the other 
hand, the second observer at late times t > r will naturally choose |0+) as the vacuum. However, 
for this second observer, the state |0_) is not the vacuum of his Hamiltonian, but actually an excited 
state that is a superposition of states with well-defined number of particles. In this sense it can be 
said that the external source has the effect of creating particles out of the "in" vacuum. Besides, 
this breaking of time translation invariance produces a violation in the energy conservation as we 
see from Eq. (9.15). Particles are actually created from the energy pumped into the system by the 
external source. 

The Schwinger effect. A classical example of creation of particles by a external field was 
pointed out by Schwinger [43] and consists of the creation of electron-positron pairs by a strong 
electric field. In order to illustrate this effect we are going to follow a heuristic argument based on 
the Dirac sea picture and the WKB approximation. 

In the absence of an electric field the vacuum state of a spin-| field is constructed by filling 
all the negative energy states as depicted in Fig. 2. Let us now connect a constant electric field 



a vacuum expectation value for 




(9.15) 
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£ = — £u x in the range < x < L created by a electrostatic potential 

x < 

V(f)={ £(x-x ) 0<x<L (9.16) 

x > L 

After the field has been switched on, the Dirac sea looks like in Fig. 18. In particular we find that 
if £L > 2m there are negative energy states at x > L with the same energy as the positive energy 
states in the region x < 0. Therefore it is possible for an electron filling a negative energy state 
with energy close to —2m to tunnel through the forbidden region into a positive energy state. The 
interpretation of such a process is the production of an electron-positron pair out of the electric field. 

We can compute the rate at which such pairs are produced by using the WKB approximation. 
Focusing for simplicity on an electron on top of the Fermi surface near x = L with energy E , the 
transmission coefficient in this approximation is given by 22 



?wkb = exp 



cxp 



/•^(Bo+V^W) I — 

2 / dx y m 2 — [E — eS(x — x )\ + p^ 



{PT+m) 



(9.17) 



where pj, = p 2 + p z . This gives the transition probability per unit time and per unit cross section 
dydz for an electron in the Dirac sea with transverse momentum p T and energy E . To get the total 
probability per unit time and per unit volume we have to integrate over all possible values of pr and 
E . Actually, in the case of the energy, because of the relation between E and the coordinate x at 
which the particle penetrates into the barrier we can write 4|f = j^dx and the total probability per 
unit time and per unit volume for the creation of a pair is given by 

^ 2 (i)/^-* wt - , -^-< 

where the factor of 2 accounts for the two polarizations of the electron. 

Then production of electron-positron pairs is exponentially suppressed and it is only sizeable 
for strong electric fields. To estimate its order of magnitude it is useful to restore the powers of c 
and h in (9.18) 

,=2 02 2 3 

W=-^i—e-V- (9.19) 
Allien 2 

The exponential suppression of the pair production disappears when the electric field reaches the 
critical value £ ciit at which the exponent is of order one 

™2 3 

£ C rit = - 1-3 x 10 16 VcnT 1 . (9.20) 

ne 



22 Notice that the electron satisfy the relativistic dispersion relation E = \/f> 2 + m 2 + V and therefore — p 2 
m 2 — (E — V) 2 + Pt- The integration limits are set by those values of x at which p x = 0. 
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This is indeed a very strong field which is extremely difficult to produce. A similar effect, however, 
takes place also in a time- varying electric field [44] and there is the hope that pair production could 
be observed in the presence of the alternating electric field produced by a laser. 

The heuristic derivation that we followed here can be made more precise in QED. There the 
decay of the vacuum into electron-positron pairs can be computed from the imaginary part of the 
effective action T[A^] in the presence of a classical gauge potential 

1 - iefi- 



= log det 



— m 



(9.21) 



This determinant can be computed using the standard heat kernel techniques. The probability of 
pair production is proportional to the imaginary part of «r[A M ] and gives 



W = 



JC2 ^ 1 2 



4tt 3 ^ n 

n=l 



Our simple argument based on tunneling in the Dirac sea gave only the leading term of Schwinger's 
result (9.22). The remaining terms can be also captured in the WKB approximation by taking into 
account the probability of production of several pairs, i.e. the tunneling of more than one electron 
through the barrier. 

Here we have illustrated the creation of particles by semiclassical sources in Quantum Field 
Theory using simple examples. Nevertheless, what we learned has important applications to the 
study of quantum fields in curved backgrounds. In Quantum Field Theory in Minkowski space-time 
the vacuum state is invariant under the Poincare group and this, together with the covariance of 
the theory under Lorentz transformations, implies that all inertial observers agree on the number of 
particles contained in a quantum state. The breaking of such invariance, as happened in the case of 
coupling to a time-varying source analyzed above, implies that it is not possible anymore to define 
a state which would be recognized as the vacuum by all observers. 

This is precisely the situation when fields are quantized on curved backgrounds. In particular, 
if the background is time-dependent (as it happens in a cosmological setup or for a collapsing star) 
different observers will identify different vacuum states. As a consequence what one observer call 
the vacuum will be full of particles for a different observer. This is precisely what is behind the 
phenomenon of Hawking radiation [45]. The emission of particles by a physical black hole formed 
from gravitational collapse of a star is the consequence of the fact that the vacuum state in the 
asymptotic past contain particles for an observer in the asymptotic future. As a consequence, a 
detector located far away from the black hole detects a stream of thermal radiation with temperature 

he 3 

THawking = MG N kM (9 ' 23) 
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where M is the mass of the black hole, G N is Newton's constant and k is Boltzmann's constant. 
There are several ways in which this results can be obtained. A more heuristic way is perhaps to 
think of this particle creation as resulting from quantum tunneling of particles across the potential 
barrier posed by gravity [46]. 

9.2 Supersymmetry 

One of the things that we have learned in our journey around the landscape of Quantum Field 
Theory is that our knowledge of the fundamental interactions in Nature is based on the idea of 
symmetry, and in particular gauge symmetry. The Lagrangian of the Standard Model can be written 
just including all possible renormalizable terms (i.e. with canonical dimension smaller o equal to 4) 
compatible with the gauge symmetry SU(3)xSU(2)xU(l) and Poincare invariance. All attempts to 
go beyond start with the question of how to extend the symmetries of the Standard Model. 

As explained in Section 5.1, in a quantum field theoretical description of the interaction of 
elementary particles the basic observable quantity to compute is the scattering or S-matrix giving the 
probability amplitude for the scattering of a number of incoming particles with a certain momentum 
into some final products 

Aim — > out) = (p/, . . . ; out|pi, . . . ; in). (9.24) 

An explicit symmetry of the theory has to be necessarily a symmetry of the S-matrix. Hence it is 
fair to ask what is the largest symmetry of the S'-matrix. 

Let us ask this question in the simple case of the scattering of two particles with four-momenta 
Pi and p 2 in the t-channel 



We will make the usual assumptions regarding positivity of the energy and analyticity. Invariance 
of the theory under the Poincare group implies that the amplitude can only depend on the scattering 
angle d through 

t = (pi - Pi) 2 = 2 (ml - pi • p[) =2 (ml- E X E' X + \p x \\p{\ costf) . (9.25) 

If there would be any extra bosonic symmetry of the theory it would restrict the scattering angle to a 
set of discrete values. In this case the S'-matrix cannot be analytic since it would vanish everywhere 
except for the discrete values selected by the extra symmetry. 

Actually, the only way to extend the symmetry of the theory without renouncing to the ana- 
lyticity of the scattering amplitudes is to introduce "fermionic" symmetries, i.e. symmetries whose 
generators are anticommuting objects [47]. This means that in addition to the generators of the 
Poincare group 23 P M , M^ v and the ones for the internal gauge symmetries G, we can introduce a 

23 The generators M pi/ are related with the ones for boost and rotations introduced in section 4.1 by J 1 = M°\ 
M l = \e % i k M3 k . In this section we also use the "dotted spinor" notation, in which spinors in the (|, 0) and (0, |) 
representations of the Lorentz group are indicated respectively by undotted (a, 6, . . .) and dotted (d, b, . . .) indices. 
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number of fermionic generators Q J a , Q hI (I = 1, . . . ,Af), where Q hI = (QiY ■ The most general 
algebra that these generators satisfy is the AA-extended supersymmetry algebra [48] 

{QlQ hj } = 2a^P ll 6 I J , 

{QlQl} = 2e ah Z 1J , (9.26) 
{QIqI) = -2e ab Z U , (9.27) 

where Z IJ E C commute with any other generator and satisfies Z IJ = —Z JI . Besides we have the 
commutators that determine the Poincare transformations of the fermionic generators Qi,Qaj 

iQin = iQai,n = o, 

[Ql,Mn = \(<rnaQl (9-28) 

where a 0i = —ia\ a 1 ^ = e^ k a k and a^ u = (a^Y . These identities simply mean that Q^, Q dJ 
transform respectively in the (|, 0) and (0, |) representations of the Lorentz group. 

We know that the presence of a global symmetry in a theory implies that the spectrum can be 
classified in multiplets with respect to that symmetry. In the case of supersymmetry start with the 
case case J\f = 1 in which there is a single pair of supercharges Q a , Q h satisfying the algebra 

{Q a ,Q b } = 2^ ab P„ {Q a ,Q b } = {Qa,Q b } = 0. (9.29) 

Notice that in the J\f = 1 case there is no possibility of having central charges. 

We study now the representations of the supersymmetry algebra (9.29), starting with the mass- 
less case. Given a state | k) satisfying k 2 = 0, we can always find a reference frame where the 
four- vector k 11 takes the form k^ = (E, 0, 0, E). Since the theory is Lorentz covariant we can obtain 
the representation of the supersymmetry algebra in this frame where the expressions are simpler. In 
particular, the right-hand side of the first anticommutator in Eq. (9.29) is given by 

K, P > = 2 ( P ° " ff3p3 ) = ( o ^E ) ■ (930) 
Therefore the algebra of supercharges in the massless case reduces to 

{Q U Q\} = {Q 1 ,Ql} = 0, 

{Q 2 ,Q\} = AE. (9.31) 

The commutator {Qi,Q\} = implies that the action of Qi on any state gives a zero-norm state of 
the Hilbert space \\Qi || = 0. If we want the theory to preserve unitarity we must eliminate these 
null states from the spectrum. This is equivalent to setting Qi = 0. On the other hand, in terms of 
the second generator Q 2 we can define the operators 

Q 2 , at = -±=Ql (9.32) 



2VE ' 2y/E 
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which satisfy the algebra of a pair of fermionic creation- annihilation operators, {a, al} = 1, a 2 = 
(a)) 2 = 0. Starting with a vacuum state a\ A) = with helicity A we can build the massless multiplet 

|A>, |A+ J> = a+IA). (9.33) 

Here we consider two important cases: 

- Scalar multiplet: we take the vacuum state to have zero helicity |0 + ) so the multiplet consists 
of a scalar and a helicity- | state 

|0 + ), |i) = a t |0 + >. (9.34) 

However, this multiplet is not invariant under the CPT transformation which reverses the sign 
of the helicity of the states. In order to have a CPT-invariant theory we have to add to this 
multiplet its CPT-conjugate which can be obtain from a vacuum state with helicity A = — \ 

|0->, |-|). (9.35) 

Putting them together we can combine the two zero helicity states with the two fermionic 
ones into the degrees of freedom of a complex scalar field and a Weyl (or Majorana) spinor. 

- Vector multiplet: now we take the vacuum state to have helicity A = \, so the multiplet 
contains also a massless state with helicity A = 1 

||), |l) = a f ||)- (9-36) 

As with the scalar multiplet we add the CPT conjugated obtained from a vacuum state with 
helicity A = — 1 

|-|), 1-1), (9-37) 
which together with (9.36) give the propagating states of a gauge field and a spin-| gaugino. 

In both cases we see the trademark of supersymmetric theories: the number of bosonic and fermionic 
states within a multiplet are the same. 

In the case of extended supersymmetry we have to repeat the previous analysis for each su- 
persymmetry charge. At the end, we have hf sets of fermionic creation- annihilation operators 
{a 7 ,a|} = S 1 j, (a/) 2 = (a|) 2 = 0. Let us work out the case of J\f = 8 supersymmetry. Since 
for several reasons we do not want to have states with helicity larger than 2, we start with a vacuum 
state | — 2) of helicity A = —2. The rest of the states of the supermultiplet are obtained by applying 



105 



the eight different creation operators a\ to the vacuum: 



A = 2 : a\ ... 4| - 2) ^ g J = 1 state, 

A = ^: a\ 1 ...a\ r \-2) Q = 8 states, 



A = 1 : a\ . . . a} 6 \ - 2} =28 states, 



1 



A = ^ : a ] h . . . a\ 5 \ - 2) I " J = 56 states, 



2 11 J5 ' ' V5 



A = : . . . a{j - 2) ( " ) = 70 states, (9.38) 



A = ~ : alx^O/gl - 2 ) ( 3 ) = 56 states > 



A = -1 : a} x a} 2 1 - 2) ( ° ) = 28 states, 



A = -- : a{j - 2) ( n J = 8 states, 



3 

"2" V 1 
A = -2 : I - 2) 1 state. 

Putting together the states with opposite helicity we find that the theory contains: 

- 1 spin- 2 field g^ u (a graviton), 

- 8 spin- 1 gravitino fields ip 1 , 

- 28 gauge fields A [ "\ 

- 56 spin- \ fermions ^ IJK \ 

- 70 scalars ^ IJKL \ 

where by [IJ...] we have denoted that the indices are antisymmetrized. We see that, unlike the 
massless multiplets of M = 1 supersymmetry studied above, this multiplet is CPT invariant by 
itself. As in the case of the massless Af = 1 multiplet, here we also find as many bosonic as 
fermionic states: 

bosons: 1 + 28 + 70 + 28 + 1 = 128 states, 
fermions: 8 + 56 + 56 + 8 = 128 states. 

Now we study briefly the case of massive representations \k), k 2 = M 2 . Things become 
simpler if we work in the rest frame where P° = M and the spatial components of the momentum 
vanish. Then, the supersymmetry algebra becomes: 

{QlQ, J } = 2M5.5 I J . (9.39) 
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We proceed now in a similar way to the massless case by defining the operators 



< s im Q '- ai - s im Q "- (9 - 40) 

The multiplets are found by choosing a vacuum state with a definite spin. For example, for M = 1 
and taking a spin-0 vacuum |0) we find three states in the multiplet transforming irreducibly with 
respect to the Lorentz group: 

|0>, 4|0), e^4ot|0>, (9.41) 

which, once transformed back from the rest frame, correspond to the physical states of two spin-0 
bosons and one spin-| fermion. For A/"-extended supersymmetry the corresponding multiplets can 
be worked out in a similar way. 

The equality between bosonic and fermionic degrees of freedom is at the root of many of 
the interesting properties of supersymmetric theories. For example, in section 4 we computed the 
divergent vacuum energy contributions for each real bosonic or fermionic propagating degree of 
freedom is 24 

E yac = ±^5(6) J d 3 pu p , (9.42) 

where the sign ± corresponds respectively to bosons and fermions. Hence, for a supersymmet- 
ric theory the vacuum energy contribution exactly cancels between bosons and fermions. This 
boson-fermion degeneracy is also responsible for supersymmetric quantum field theories being less 
divergent than nonsupersymmetric ones. 



Appendix: A crash course in Group Theory 

In this Appendix we summarize some basic facts about Group Theory. Given a group G a represen- 
tation of G is a correspondence between the elements of G and the set of linear operators acting on 
a vector space V, such that for each element of the group g e G there is a linear operator D(g) 

D(g) :V^V (A.43) 

satisfying the group operations 

D( 9l )D(g 2 ) = D(g l92 ), D(g^) = D(g 1 )~ 1 , 9l ,g 2 G Q. (A.44) 

The representation D{g) is irreducible if and only if the only operators A : V — > V commuting with 
all the elements of the representation D{g) are the ones proportional to the identity 

[D(g),A] — 0, V<? <^ A = X1, AeC (A.45) 

24 For a boson, this can be read off Eq. (3.56). In the case of fermions, the result of Eq. (4.44) gives the vacuum 
energy contribution of the four real propagating degrees of freedom of a Dirac spinor. 
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More intuitively, we can say that a representation is irreducible if there is no proper subspace U C V 
(i.e. U ^ V and U ^ 0) such that D(g)U C U for every element g e G. 

Here we are specially interested in Lie groups whose elements are labelled by a number of 
continuous parameters. In mathematical terms this means that a Lie group is a manifold Ai together 
with an operation Ai x Ai — > Ai that we will call multiplication that satisfies the associativity 
property g\ ■ (g 2 ■ #3) = (gi ■ g 2 ) • #3 together with the existence of unity gl = lg = g,for every 
g G Ai and inverse gg^ 1 = g~ l g = 1. 

The simplest example of a Lie group is SO(2), the group of rotations in the plane. Each 
element R(9) is labelled by the rotation angle 9, with the multiplication acting as R(9i)R(9 2 ) = 
R{9\+9 2 ). Because the angle 9 is defined only modulo 2n, the manifold of SO(2) is a circumference 
S 1 . 

One of the interesting properties of Lie groups is that in a neighborhood of the identity element 
they can be expressed in terms of a set of generators T a (a = 1, . . . , dim G) as 

00 (—i) n 

D(g) = exp[-ta a T a ] = ^ K -^-a ai . . . a a T a > . . . T a ", (A.46) 

n=0 

where a a E C are a set of coordinates of .M in a neighborhood of 1. Because of the general Baker- 
Campbell-Haussdorf formula, the multiplication of two group elements is encoded in the value of 
the commutator of two generators, that in general has the form 

[T a ,T b }=if abc T c , (A.47) 

where f abc 6 C are called the structure constants. The set of generators with the commutator 
operation form the Lie algebra associated with the Lie group. Hence, given a representation of the 
Lie algebra of generators we can construct a representation of the group by exponentiation (at least 
locally near the identity). 

We illustrate these concept with some particular examples. For SU(2) each group element 
is labelled by three real number cti, i = 1,2,3. We have two basic representations: one is the 
fundamental representation (or spin \) defined by 

Di(aii) = e~^ aia \ (A.48) 
2 

with a 1 the Pauli matrices. The second one is the adjoint (or spin 1) representation which can be 
written as 

D l (a l ) = e~ iaiJi , (A.49) 

where 

/ \ / -1 \ / 1 \ 

J 1 = I 1 , J 2 = , J 3 = -1 . (A.50) 

\ -1 / \ 1 / \ / 
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Actually, J 1 (i = 1,2,3) generate rotations around the x, y and z axis respectively. Representations 
of spin jGN + i can be also constructed with dimension 



dimD J (g)=2j + l. 



(A.51) 



As a second example we consider SU(3). This group has two basic three-dimensional repre- 
sentations denoted by 3 and 3 which in QCD are associated with the transformation of quarks and 
antiquarks under the color gauge symmetry SU(3). The elements of these representations can be 
written as 



D s (a a ) = e^ a x \ Dz(a a ) = e-2 
where A a are the eight hermitian Gell-Mann matrices 



(a = l,...,8), 



(A.52) 



Ai = 



A, 



A 7 = 




A, 



A, 



A, 




A, = 



A fi 




(A.53) 



Hence the generators of the representations 3 and 3 are given by 

T a (3) = l -\ a , T a (3) -- 



V 



(A.54) 



Irreducible representations can be classified in three groups: real, complex and pseudoreal. 

- Real representations: a representation is said to be real if there is a symmetric matrix S which 
acts as intertwiner between the generators and their complex conjugates 



-ST a S 



a q— 1 



S T = S. 



(A.55) 



This is for example the case of the adjoint representation of SU(2) generated by the matrices 
(A.50) 

- Pseudoreal representations: are the ones for which an antisymmetric matrix S exists with the 
property 



T a = -ST a S-\ S T = -S. 



(A.56) 



As an example we can mention the spin-| representation of SU(2) generated by \a % . 
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- Complex representations: finally, a representation is complex if the generators and their com- 
plex conjugate are not related by a similarity transformation. This is for instance the case of 
the two three-dimensional representations 3 and 3 of SU(3). 

There are a number of invariants that can be constructed associated with an irreducible repre- 
sentation R of a Lie group G and that can be used to label such a representation. If T R are the gen- 
erators in a certain representation R of the Lie algebra, it is easy to see that the matrix T R T R 
commutes with every generator T R . Therefore, because of Schur's lemma, it has to be proportional 
to the identity 25 . This defines the Casimir invariant C 2 (R) as 

dimG 

Y,T a R T a R = C 2 {R)l. (A.57) 

a=l 

A second invariant T 2 (R) associated with a representation R can also be defined by the identity 

Tr T R T R = T 2 (R)5 ab . (A.58) 

Actually, taking the trace in Eq. (A.57) and combining the result with (A.58) we find that both 
invariants are related by the identity 

C 2 (R) dim R = T 2 (R) dimG, (A.59) 

with dim R the dimension of the representation R. 

These two invariants appear frequently in quantum field theory calculations with nonabelian 
gauge fields. For example T 2 (R) comes about as the coefficient of the one-loop calculation of the 
beta-function for a Yang-Mills theory with gauge group G. In the case of SU(N), for the fundamental 
representation, we find the values 

N 2 - 1 1 
C 2 (fund) = -^ r , T 2 (fund) = -, (A.60) 

whereas for the adjoint representation the results are 

C 2 (adj) = N, T 2 (adj) = N. (A.61) 

A third invariant A(R) is specially important in the calculation of anomalies. As discussed 
in section (7), the chiral anomaly in gauge theories is proportional to the group-theoretical factor 
Tr [T^{T R , T R }~\ . This leads us to define A(R) as 

Tr [T a R {Tl T R }~\ = A(R), d abc (A.62) 

where d abc is symmetric in its three indices and does not depend on the representation. Therefore, 
the cancellation of anomalies in a gauge theory with fermions transformed in the representation R 
of the gauge group is guaranteed if the corresponding invariant A(R) vanishes. 

25 Schur's lemma states that a representation of a group is irreducible if and only if all matrices commuting with every 
element of the representation are proportional to the identity. 
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It is not difficult to prove that A(R) = if the representation R is either real or pseudoreal. 
Indeed, if this is the case, then there is a matrix S (symmetric or antisymmetric) that intertwins the 
generators T R and their complex conjugates T R = —ST^S^ 1 . Then, using the hermiticity of the 
generators we can write 



Tr [T R {T R , T r }] — Tr [T°{T b R , T C R }] T = Tr T R {T R ,T R } 



(A.63) 



Now, using (A. 55) or (A. 56) we have 



Tr 



T R {T R , T C R }] = -Tr [ST^S-^STiS- 1 , ST^S' 1 }] = -Tr [T a R {T b R , T C R }} , (A.64) 



which proves that Tr [T R {T R , T r }] and therefore A(R) = whenever the representation is real or 
pseudoreal. Since the gauge anomaly in four dimensions is proportional to A(R) this means that 
anomalies appear only when the fermions transform in a complex representation of the gauge group. 
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